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Abstract 



This dissertation presents investigations on dynamics of discrete-time quantum 
walk and some of its applications. Quantum walks has been exploited as an useful 
tool for quantum algorithms in quantum computing. Beyond quantum computa- 
tional purposes, it has been used to explain and control the dynamics in various 
physical systems. In order to use the quantum walk to its fullest potential, it is 
important to know and optimize the properties purely due to quantum dynam- 
ics and in presence of noise. Various studies of its dynamics in the absence and 
presence of noise have been reported. We propose new approaches to optimize the 
dynamics, discuss symmetries and effect of noise on the quantum walk. Making 
use of its properties, we propose the use of quantum walk as an efficient new tool 
for various applications in physical systems and quantum information processing. 

In the first and second part of this dissertation, we discuss evolution process of 
the quantum walks, propose and demonstrate the optimization of discrete-time 
quantum walk using quantum coin operation from SU(2) group and discuss some 
of its properties. We investigate symmetry operations and environmental effects on 
dynamics of the walk on a line and an n— cycle highlighting the interplay between 
noise and topology. 

Using the properties and behavior of quantum walk discussed in part two, in 
part three we propose the application of quantum walk to realize quantum phase 
transition in optical lattice, that is to efficiently control and redistribute ultracold 
atoms in optical lattice. We also discuss the implementation scheme. Another 
application we consider is creation of spatial entanglement using quantum walk 
on a quantum many body system. 
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Preface 



This dissertation which is divided into four parts discusses the dynamics of discrete- 
time quantum walk and some of its applications. Part I, Introduction consists of 
Chapter 1 which introduces the reader to an overview of the subject followed by 
brief discussion of results presented in the latter chapters. Part II, Quantum walks 
and its dynamics consists of Chapters 2 and 3. In Chapter 2, we review the two 
standard forms of quantum walks (continuous- and discrete-time quantum walk) 
and study the dynamics of discrete-time quantum walk in detail. In Chapter 3, we 
study symmetries and effect of noise on the dynamics of the discrete-time quan- 
tum walk. Part III, Applications consists of Chapters 4 and 5. In Chapter 4, 
we demonstrate the use of discrete-time quantum walk to control the dynamics 
of ultracold in optical lattice and in Chapter 5, we create spatial entanglement 
using quantum walk on distinguishable multiparticle system. Part IV, Conclusion 
consists of Chapter 6 which concludes with summery and future perspective. 
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Chapter 1 
Introduction 



1.1 Introduction 



Theoretical studies and experimental evidences from the early and mid twenti- 
eth century led us to consider the physical world is governed by laws of quantum 
mechanics (see any standard text book on quantum mechanics. For example, Feyn- 
man (1970) [1] and R. Shenkar (1994) [2]). This, in early 1980's led Yuri Manin 
(1980) [3] and Feynman (1982) [4] to independently observe and suggest that the 
physical world can be ideally simulated using quantum computers. A decade later, 
theoretical studies and the experimental possibility to construct quantum comput- 
ers based on the laws of quantum mechanics such as superposition and interference 
between quantum amplitudes became one of the active areas of research [5, 6]. 
Deutsch algorithm (1985) [7] and Deutsch-Jozsa algorithm (1992) [8] were among 
the first algorithms to show that quantum computers are capable of solving cer- 
tain computational problems much more efficiently than classical deterministic 
computers. Simon's algorithm (1994) was one of the first example to show that a 
quantum algorithm can solve the problem of computing an unknown function with 
a polynomial number of queries, where a classical algorithm required an exponen- 
tial number of queries [9, 10]. In 1994, Shor devised a quantum algorithm to factor 
arbitrary integers exponentially faster than the best known classical counterpart 
[11, 12]. In 1996, Grover devised an algorithm which can in principle search an 
unsorted database quadratically faster than any known classical algorithm [13]. 
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Shor's algorithm unleashed active research in the area of quantum information 
and quantum computation across a broad range of disciplines: quantum physics, 
computers science, mathematics and engineering. This research has unveiled many 
new effects that are strikingly different from their classical counterparts. Design 
and analysis of quantum algorithms to solve various problems more efficiently than 
classical algorithms has become a vibrant research area (see Mosca (2008) [14] and 
Childs and van Dam (2008) [15] for recent reviews). 

Among classical algorithms, many are based on classical random walk. Markov 
chain simulation, which has emerged as a powerful algorithmic tool [16] is one such 
example. Like classical random walk, the quantum version of it has also become 
an important constituent of quantum algorithms and computation. 

The quantum walk as it is known today is a generalization of the classical random 
walk developed by using the aspects of quantum mechanics such as superposition 
and interference of quantum amplitudes. In the classical random walk the particle 
moves in the position space with a certain probability, whereas the quantum walk, 
which involves a superposition of states, moves by exploring multiple possible paths 
simultaneously with the amplitudes corresponding to different paths interfering. 
This makes the variance of the quantum walk on a line to grow quadratically with 
the number of steps (time), compared to the linear growth for the classical random 
walk [17, 18]. A probabilistic result is obtained in quantum walk upon measure- 
ment. Figure 1.1 shows the comparative difference in the probability distribution 
of classical random walk and quantum walk. 

Like its classical counterpart, which has found applications in almost all branches 
of sciences [19, 20], quantum walk has also proven to be a useful tool for quantum 
algorithms, to study, control and explain the dynamics in various physical systems. 
Below, we will briefly summarize the history of quantum walk and the progress 
made in employing it for quantum algorithms and the study of the dynamics of 
physical systems. 

Though quantum random walk^ was first introduced by Aharonov et al. (1993) 
[21], the idea of exploring multiple paths simultaneously using path integrals dates 

^Unlike its classical counterpart, the evolution of the quantum version is unitary, reversible 
and has no randomness associated with it during the evolution. Therefore, keeping away the 
term 'random', quantum walk has been the preferred usage. 
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-Quantum walk 
-Classical random walk 




-100 -80 -60 -40 -20 20 40 60 80 100 
Particle position 



Figure 1.1: Probability distribution of classical random walk and its quantum 
counterpart, quantum walk. Difference in the variance can been seen in the 
distribution. The distribution is after 100 steps of walk. 



back to works by Riazanov (1958) [22] and Feynman and Hibbs (1965) [23]. Meyer 
(1996) studied the dynamics of quantum particle in a quantum lattice gas formu- 
lation. The dynamics involved the interference of the left and the right moving 
components of the amplitude at given nodes [24, 25]. Building on these ideas, the 
concept of quantum walk was developed and now is studied in two standard forms: 
continuous-time quantum walk and discrete-time quantum walk. Continuous-time 
quantum walk was introduced by Farhi and Gutmann (1998) [26] and discrete-time 
quantum walk was introduced by Watrous (2001) [27]. A specific form of discrete- 
time quantum walk know as Hadamard walk was introduced by Ambainis et al. 
(2001) [17]. An early review by Kempe (2003) [28] discusses the two variants in 
detail with some algorithmic applications. 

Both continuous- and discrete-time quantum walk have been widely used in algo- 
rithms for a variety of problems. Several quantum search algorithms have been 
proposed (see for example, Shenvi et al. (2003) [29]; Childs and Goldstone (2004) 
[30, 31]; Ambainis et al. (2005) [32]; Aaronson and Ambainis (2005) [33]; Magniez 
et al. (2007) [34]). The continuous-time quantum walk has been used by Childs 
et al. (2003) [35] to demonstrate the exponential speedup over classical compu- 
tation for a hitting time problem on a glued tree. Ambainis (2004) [36] review 
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discusses some of these algorithms in detaiL Ambainis (2007) [37] apphed quantum 
walk to give an optimal quantum algorithm for the element distinctness problem. 
It has been applied to various other problems in query model like triangle finding 
by Magniez et al. (2005) [38], checking matrix multiplication by Buhrman and 
Spalek (2006) [39], testing group commutativity by Maganiez and Nayak (2007) 
[40], for evaluating balanced binary game trees by Farhi et al. (2007) [41], Boolean 
formulas by Ambainis et al. (2007) [42] and Reichardt and Spalek [43] and to 
obtain general adversary bound to characterize the quantum query complexity by 
Reichardt [44]. 

Beyond applications in quantum algorithms, quantum walk is emerging as a po- 
tential tool to understand various phenomenon in physical systems and has been 
employed to demonstrate coherent control over quantum many body systems. Oka 
et al. (2005) [45] mapped the breakdown of electric-field driven system to quan- 
tum walk, Engel et al. (2007) [46] and Mohseni et al. (2008) [47] explained 
wavelike energy transfer within photosynthetic systems, Somma et al. (2008) [48] 
described quantum simulation of classical annealing processes, and Chandrashekar 
and Laflamme (2008) [49] demonstrated coherent quantum control over redistri- 
bution of atoms in optical lattice. 

Some experimental progress on the implementation of quantum walk has been re- 
ported. Du et al. (2003) [50] and Ryan et al. (2005) [51] implemented continuous- 
time quantum walk on 2- qubit and discrete-time quantum walk on a 3- qubit 
Nuclear Magnetic Resonance (NMR) system respectively. Grossman et al. (2004) 
[52] with sodium Bose-Einstein condensates and recently, Perets et al. (2008) [53] 
using propagating photons in waveguide lattices have implemented the discrete- 
time quantum walk. Other implementation schemes have also been proposed: 
Travaglione and Milburn (2002) [54] in an ion trap, Diir et al. (2002) [55] and 
Eckert et al. (2005) [56] on neutral cold atoms in an optical traps, Chandrashekar 
(2006) [57] using Bose-Einstein condensates. Ma et al. (2006) [58] using quan- 
tum accelerator mode, Manouchehri and Wang (2008) [59] in an array of quantum 
dots and Sanders et al. (2003) [60] introduced quantum quincunx to physically 
demonstrate quantum walk in cavity quantum electrodynamics capabilities. 
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It is important to understand the dynamics and behavior of quantum walk, both 
in the presence and absence of environmental effects to effectively use it for devel- 
oping algorithms, to simulate any natural physical process or to experimentally 
implement in any physical system. Considerable study has been done in this di- 
rection. Since this thesis involves results of some of these investigations, we will 
first list the main work that constitutes this thesis and briefly elaborate on them 
in the following sections. 

The main constituents of this thesis are : 

• Structure and dynamic properties of discrete-time quantum walk. Construc- 
tion of a generalized model using quantum coin operation from SU(2) group 
to demonstrate control over the dynamics of the walk [61, 62]. 

• Using a coin-embedded shift operator, a generic quantum walk model which 
can be used to conveniently retrieve discrete- and continuous-time quantum 
walk under different conditions is proposed. The coin degree of freedom is 
retained while the standard variants are retrieved and the generic model sim- 
plifies the physical resources required to implement quantum walk compared 
to the resource required for the discrete-time quantum walk [63]. 

• Investigations on symmetries in discrete-time quantum walk on a line and its 
behavior in presence of environmental effects, breakdown of these symmetries 
for a walk on an ra— cycle and restoration with noise. Relevance of these 
investigations in physical systems [64, 65]. 

• Demonstration of coherent control over the redistribution of atoms in op- 
tical lattice, use of noisy channels- bit flip, phase damping, and amplitude 
damping to act as an additional toolbox to control the dynamics of atoms 
and an experimental implementation scheme [49, 57]. 

• Spatial entanglement in many body system using quantum walk and its 
physical relevance [66]. 
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1.2 Dynamics of quantum walk 

In the continuous-time quantum walk, one can directly define the walk on the 
position Hilbert space Tip [26]. In the discrete-time quantum walk, in addition to 
Hp it is necessary to introduce a quantum coin operation, a coin Hilbert space 
Tic to define the direction in which the particle has to move [17]. Due to the 
coin degree of freedom, the discrete-time variant is shown to be more powerful 
than the other in some contexts (see Ambainis et al. (2005) [32]). To match 
the performance of a spatial search using the discrete-time quantum walk, the 
coin degree of freedom has been introduced in the continuous-time quantum walk 
model [31]. In this thesis, we mainly consider the discrete-time quantum walk for 
our study. 

The dynamics of discrete-time quantum walk using the Hadamard operation 



as quantum coin operation was analyzed and defined as Hadamard walk [17]. Most 
of the work thereafter on discrete-time quantum walk considered using Hadamard 
operation as quantum coin operation. Nayak and Vishwanath (2001) [18] analyzed 
using the single-variable parameter U(2) operator 



as a quantum coin^. The variance cr^ was shown to be dependent on the parameter 
6. The discrete-time quantum walk on a particle initially in a symmetric superpo- 
sition state (|0) -|- at position space (I'lpo)), l^ms) = (|0) + ® I'lpo) 
using Be as the quantum coin was shown to return the symmetric probability 
distribution of the walk in the position space. Asymmetry in the initial state of 
the particle introduced an asymmetric probability distribution. This led to their 
conclusion that obtaining a symmetric distribution depends largely on the initial 
state of the particle [17, 18, 67]. 

^Meyer has used the U(2) operator as a scattering operator in lattice quantum gas automata 







model [24] 



8 



Chapter 1. Introduction 



In Chapter 2, we will first define the continuous-time and discrete-time quantum 
walk, equivalence of their form to the Schrodinger and Dirac equations respectively. 
In Section 2.3.3, we present the discrete-time quantum walk using operation from 
the SU(2) group with three Caley-Klein parameters ^, 6 and ( (Euler angles) as 
the quantum coin 



We show that the variance is dependent on the parameter 6 and takes the form 
0"^ ~ (1 — sin(6'))t^, where t is number of steps of the quantum walk. For a particle 
with symmetric superposition as the initial state parameters ^ and ( introduce 
asymmetry in the probability distribution and their effect on the variance is very 
small. For a particle with an asymmetric superposition as the initial state, the 
parameters ^ and ( can be configured to obtain a symmetric probability distribu- 
tion. We discuss the variation of measurement entropy in position space with the 
coin parameters and optimization of the quantum walk for the maximum variance, 
improving mixing time in an n-cycle and controlling the probability distribution 
using coin parameters [61]. 

Recurrence in the dynamics of physical systems is an important phenomenon. 
For a classical conservative system, whether discrete or continuous in time, the 
Poincare recurrence theorem states that any phase-space configuration of a system 
enclosed in a finite volume will be repeated as accurately as one wishes after a finite 
interval of time (with no restriction on the interval) [68, 69]. A similar recurrence 
theorem is shown to hold in quantum theory as well [70, 71]. In a system with a 
discrete energy eigenvalue spectrum {En}; if ^I/(to) is its state vector at the time 
to and e is any positive number, there exists a finite time T such that. 



Robinett's (2004) [72] review article discusses recurrence in quantum systems. In 
Section 2.4, we show the discrete-time quantum walk which evolves with inter- 
ference of quantum amplitudes fails to recur completely in position space but 
fractional recurrence is seen [62]. 




(1.3) 
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Chapter 3 focuses on the symmetries and effects of noise on the dynamics of 
discrete-time quantum walk. The transition of walk from the quantum to the 
classical regime has been studied earlier by introducing decoherence to Hadamard 
walk [73-75] (see review by Kendon (2006) [76]). In a scheme proposed by Chan- 
drashekar [57] to implement the discrete-time quantum walk on Bose-Einstein 
condensates, stimulated Raman kick was used as shift operator. The stimulated 
Raman kick, in addition to shifting different states of atoms to different position 
space flips (bit-flip) the internal state of the atoms. However in spite of bit-flip 
on state of atoms along with shift in position space, the probability distribution 
remained invariant in position space. Motivated by this, the shift operator was 
augmented by other form of operations and the operation for which the distribu- 
tion remains invariant have been studied and are called symmetries of quantum 
walk. These symmetries have been studied using generalized SU(2) coin operation 
Bo£fi as quantum coin. We further generalize the observations of these symmetries 
in the presence of environmental effects, modeled by various noise channels such 
as bit flip, phase-flip and generalized amplitude damping channels [64]. Interest- 
ingly, we find that the symmetry operations are sensitive to the walk topology. 
For example, a symmetry that holds for quantum walk on an one-dimensional line 
does not hold, in general, for a quantum walk on an n— cycle but leads to other 
interesting behavior. Noise on an n— cycle tends to restore these symmetry, both 
by "classicalizing" the walk and also desensitizing the symmetry operation as a 
topology probe for the quantum walk [65]. From Chapter 2, we can conclude that 
the symmetries and addition of small amount of experimentally engineered noise 
or environmental effect can be used as an additional tool to control the dynamics 
and hence the probability distribution of the quantum walk. These observations 
can have important implications for a better insight into quantum walk, and for 
simplifying certain implementations and are also of relevance in the study of dy- 
namics in environment-assisted quantum systems. 
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1.3 Applications 



By employing the properties of quantum walks due to purely quantum dynamics 
and small amount of engineered noise presented in Chapters 2 and 3, we demon- 
strate coherent quantum control over atoms in optical lattice in Chapter 4. Tradi- 
tionally, theoretical studies of the dynamics of atoms in an optical lattice are done 
using mean field approaches [77] and quantum Monte Carlo methods [78-80]. Use 
of quantum walk can serve as an alternate method. In particular, we will consider 
the quantum phase transition from the Mott insulator (MI) to the superfluid (SF) 
state [81, 82] and vice versa. The simulation of the quantum phase transition using 
the quantum walk occurs quadratically faster in one dimension (ID) compared to 
varying the optical lattice depth and letting the atom-atom interaction follow the 
classical random walk behavior. In our study, we use the noise models discussed 
in Chapter 3 to act as an enhanced toolbox to control the redistribution of atoms 
in an optical lattice [49]. We will also discuss the experimental implementation of 
quantum walk on ultra cold atoms in an optical lattice [57]. 

Entanglement in many body system has not only been a computational resource, 
it has also been used as a signature of quantum phase transition [83-85]. With 
increase in dimension of the Hilbert space, the number of invariants and the mea- 
sures of entanglement grows exponentially so that scaling becomes impractical. In 
2002, Meyer and Wallach proposed a global entanglement measure - a polynomial 
measure of multiparticle entanglement to address the scalability problem [86]. In 
Chapter 5 we investigate the evolution of spatial entanglement, particle- number 
entanglement between regions of space in a many particle system subjected to 
quantum walk evolution using Meyer- Wallach measure. We use the coin degree of 
freedom to demonstrate the dependency of entanglement with the coin parameter, 
number of particles in the system and number of steps of quantum walk. 
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Quantum walks 



2.1 Introduction 

Quantum walks are the quantum analog of the classical random walks [21-23] 
developed using the aspects of quantum mechanics such as interference and su- 
perposition. Like their classical counterpart, the quantum walks are also widely 
studied in two forms: continuous-time quantum walk and discrete-time quantum 
walk. In the continuous-time quantum walk, one can directly define the walk on 
the position space [26], whereas in the discrete-time quantum walk, it is necessary 
to introduce a quantum coin operation to define the direction in which the par- 
ticle has to move [17]. The results from the continuous-time quantum walk and 
the discrete-time quantum walk are often similar, but due to the coin degree of 
freedom, the discrete-time variant has been shown to be more powerful than the 
other in some context [32]. To match the performance of the discrete-time quan- 
tum walk, the coin degree of freedom can be introduced in the continuous-time 
quantum walk [31]. We will only consider the discrete-time quantum walk for our 
study. 

The main focus of this chapter is to present in detail the dynamics of the discrete- 
time quantum walk and demonstrate methods to control it. This chapter is or- 
ganized as follows. For the completeness, we will first review the continuous-time 
quantum walk in Section 2.2. In Section 2.3, we will review the discrete-time 
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quantum walk before demonstrating the optimization and randomization of the 
dynamics using quantum coin operations from the SU(2) group. In Section 2.4, we 
discuss the fractional recurrence nature of the quantum walk and conclude with a 
summary in Section 2.5. 



2.2 Continuous-time quantum walk 

To define the continuous-time quantum walk, it is easier to first define the continuous- 
time classical random walk and quantize it by introducing quantum amplitudes in 
place of classical probabilities. 

The continuous-time classical random walk takes place entirely in the position 
space. To illustrate, let us define continuous-time classical random walk on the 
position space Tip spanned by a vertex set of a graph G with edges set E, 
G = (y,E). A step of the random walk can be described by a adjacency matrix 
A which transform the probability distribution over V, i.e., 

, 1 (j, k) G E 
A,,k = (2.1) 
{j,k)^E 

for every pair j, k & V. The other important matrix associated with the graph G 
is the generator matrix H given by 



djl j = k 



H 



-7 



{j,k)eE^ (2.2) 



otherwise 



where dj is the degree of the vertex j and 7 is the probability of transition between 
neighboring nodes per unit time. 

If Pj{t) denotes the probability of being at vertex j at time t then the transition 
on graph G is defined as the solution of differential equation 

±P^{t) = -J2iij,kPk{t). (2.3) 
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The solution of the differential equation is given by 

P{t) = e-"*P(0). (2.4) 

By replacing the probabilities Pj by quantum amplitudes aj(t) = {j\ip{t)) where |j) 
is spanned by the orthogonal basis of the position Hilbert space Tip and introducing 
a factor of i we obtain 

tj^ajit) = J2iij,kak{t). (2.5) 

fcev 

We can see that (2.5) is the Schrodinger equation 

^||^) = H|^). (2.6) 

Since generator matrix is an Hermitian operator, the normalization is preserved 
during the dynamics. The solution of the differential equation can be written in 
the form 

\m) = e-''''\m)- (2.7) 

Therefore, the continuous-time quantum walk is of the form of Schrodinger equa- 
tion, a non-relativistic quantum evolution. 

To implement the continuous-time quantum walk on a line, the position Hilbert 
space Hp can be written as a state spanned by the basis states where j G Z. 
The Hamiltonian H is defined such that, 

HIV;,) = -7|^,_i) + 27|V^,) - 7|V^,+i) (2.8) 

and is made to evolve with time t by applying the transformation 

U{t) = exp(-im). (2.9) 

The Hamiltonian H of the process acts as the generator matrix which will trans- 
form the probability amplitude at the rate of 7 to the neighboring sites, where 7 
is time- independent constant. 
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2.3 Discrete-time quantum walk 

We will first define the structure of the discrete-time classical random walk. The 
discrete-time classical random walk takes place on the position Hilbert space Tip 
with instruction from the coin operation. A coin flip defines the direction in which 
the particle moves and a subsequent position shift operation moves the particle in 
position space. For a walk on a line, a two sided coin with head and tail defines 
the movements to the left and right respectively. 

The discrete-time quantum walk also has a very similar structure to that of its 
classical counterpart. The coin flip is replaced by the quantum coin operation 
which defines the superposition of direction in which the particle moves simul- 
taneously. The quantum coin operation followed by the unitary shift operation 
is iterated without resorting to intermediate measurement to implement a large 
number of steps. During the walk on a line, interference between the left and the 
right propagating amplitude results in the quadratic growth of variance with the 
number of steps. 

The discrete-time quantum walk on a line is defined on a Hilbert space 



where T-Lc is the coin Hilbert space and Tip is the position Hilbert space. For a 
discrete-time quantum walk in one dimension, Tic is spanned by the basis state 
(internal state) of the particle |0) and |1) and Hp is spanned by the basis state of 
the position where j G Z. To implement the discrete-time quantum walk on 
a particle at origin in state 



the quantum coin toss operation B & U{2) which in general can be written as 



(2.10) 



= (cos(5)|0) + e*''sin(5)|l)) ® \%), 



(2.11) 



5, 




is applied. Where (Jx,(Jy and are the Pauli spin operators. Parameters of the 
coin operations (,a,P,'-f can be varied to get different superposition state of the 
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particle. That is, quantum coin operation -B(,q-,/3,7 is used to evolve the particle 
to superposition of its basis states such that it can serve as an instruction to 
simultaneously evolve the particle to the left and right of its initial position. The 
quantum coin operation is followed by the conditional unitary shift operation S 
given by 

S = ^-mioHDimpi ^ (|o)(o| ® e-^^') + (|1)(1| ® e^^') (2.13) 

where P is the momentum operator, / is the step length and |0) and |1) are the 
basis states of the particle. Therefore the operator 5* which delocalizes the wave 
packet in different basis states |0) and |1) over the position (j — 1) and (j + 1) 
when when step length / = 1 can also be written as 

S = |0)(0| ® Yl + ® E (2-14) 

The states in the new position is again evolved into the superposition of its basis 
state and the process of quantum coin toss operation -B^^a,/?,^ followed by the 
conditional unitary shift operation S, 

= 5(Ec,„,;3,7 ® 1) (2.15) 

is iterated without resorting to intermediate measurement, to realize a large num- 
ber of steps of the discrete-time quantum walk. The four variable parameters of 
the quantum coin, C,a,P and 7 (2.12) can be varied to change and control the 
probability amplitude distribution in the position space. 

Most widely studied form of the discrete-time quantum walk is the Hadamard 
walk, using Hadamard operation H as quantum coin operation and the role of 
coin operation and initial state to control the probability amplitude distribution 
has been discussed in earlier studies [17, 87]. We will discuss the Hadamard walk 
in detail in Section 2.3.2. In Section 2.3.3, we demonstrate that a three parameter 
SU(2) quantum coin operation is sufficient to describe the most general form of 
the discrete-time quantum walk. Before that, we will further analyze the structure 
of the discrete-time quantum walk. 
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2.3.1 Discrete-time quantum walk and Klein-Gordon equa- 
tion 

The standard symmetric discrete-time classical random walk leads to 

P{j, t+l) = l [P{j -l,t) + P{j + l,t)], (2. 16) 

when unit time is required each step of classical random walk, P{j, t) denotes 
the probability of finding the particle at position j at discrete time t. Subtracting 
P{j, t) from both sides of (2.16) leads to the difference equation which corresponds 
to differential equation 

The above equation is irreversible because the coin is effectively thrown away after 
each toss. It is also non-relativistic in the sense that it is not symmetric in time 
and space, and leads to a dispersion relation that is essentially non-relativistic 
[88]. In the continuum limit, (2.17) leads to the standard classical diffusion. On 
the contrary, in the discrete-time quantum walk the information of the state of the 
quantum coin in the previous step is retained and carried over to the next step. 
This makes the quantum walk reversible. 

To illustrate this, we consider the wavefunction describing the position of a particle 
and analyze how it evolves with time t. Let t be the time required to implement t 
steps of quantum walk. The two component vector of amplitudes of the particle, 
being at position j, at time t, with left (L) and right (R) moving component is 
given by 

V *fl(j,*) / 

The dynamics for driven by single parameter quantum coin operation is : 

^ -^sin(^) cos{9) ' 
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and followed by the conditional shift operator S, [S{B ® 1)] in terms of left (L) 
moving and right (R) moving component at a given position j and time t + 1 is 
given by 



(2.20) 



^LU,t + l) \ / cos{9)a -zsin(^)at \ / ^^(j,^) 
y V -2sin(^)at cosi9)a J \^Rij,t) 

where action of operator a and on ^^(j, t) is given by 

a^{j,t) = ^{j + l,t) 

at^(j, t) = ^(j-l,t). (2.21) 

Therefore, 

^i(j,t + l) = cos(^)^L(j + l,t)-^sin(^)^R(j-l,t) (2.22a) 

^^(j,t + l) = cos(^)^^(j-l,t)-2sin(^)^i(j + l,t). (2.22b) 

We thus find that the coin degree of freedom is carried over during the dynamics 
of the discrete-time quantum walk making it reversible. 

Further, the components and \E^/j are decoupled to obtain 

^R{j,t + 1) + ^R{j,t- 1) = cos(^)[^«(j + l,t) + ^r{j - l,t)]. (2.23) 

Subtracting 2'^R{j,t) from both sides, we obtain a difference equation which cor- 
responds to differential equation 



^RU,t) = 2(1 -cos(^))^^(j,t), (2.24) 



and similar expression can be obtained for \E'L(j, t), see Appendix A for interme- 
diate steps. This shows that each component follows a Klein-Gordon equation of 
the form 



V^-i^)0-/.V = O (2.25) 



showing the essentially relativistic character of the discrete-time quantum walk. 
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Equations (2.22a) and (2.22b) can also be written as a compact first-order dif- 
ference equations, which in the continuum hmit would be analogous to the one- 
dimensional Dirac equation for a spinor (of spin 1/2). 

Setting the time-step and spatial-step to 1, we obtain from (2.24), the equivalent 
of light speed c and mass m in the discrete-time quantum walk dynamics : 



c 

mc 



/i = — = ^/2\^^^^(e)^]. (2.26) 



Considering ^ = 1, we can write 



m 



2[sec(^) - 1] _ ^2.27) 



cos(^) 



Thus, from Section 2.2 and the above analysis we conclude that a continuous-time 
quantum walk takes the form of Schrodinger equation and a discrete-time quantum 
walk takes the form of the relativistic quantum equation of a free spin-0 particle. 



2.3.2 Hadamard walk 

The simplest version of the discrete-time quantum walk is the walk using Hadamard 
operation as quantum coin operation and is known as the Hadamard walk [17]. 
A particle at origin in one of the basis state |0) or |1) of He (internal state of 
the particle) is evolved into the superposition of Tic with equal probability, by 
applying the Hadamard operation 



(2.28) 




such that 



{H ® 1)(|0) ® l^o)) = (|0) + |1)) 8) l^o) (2.29a) 



V2 
1 

7! 



{H ® 1)(|1) ® l^o)) = (|0) - |1)) ® li'o). (2.29b) 
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The operation H is then followed by the conditional shift operation 5* (2.14), W = 
S{H®1). The process is iterated without resorting to intermediate measurements 
to evolve the particle in superposition of position space and realize a large number 
of steps of the Hadamard walk. After the first two iterations of W, the left and 
right evolving components of the amplitude begin to interfere, deviating from the 
classical evolution, resulting in a quadratic speedup in the growth of the variance. 
The probability amplitude distribution arising from the iterated application of W 
is significantly different from the probability distribution of the classical random 
walk. Figure 2.1. The particle initially in the state |0) drifts to the left and the 
particle with an initial state |1) drifts to the right. This asymmetry arises from 
the fact that the Hadamard operation treats the two states |0) and |1) differently, 
phase difference of —1 in case of state |1). This phase difference, depending on 
the initial state of the particle contributes to the constructive interference on one 
side and to the destructive interference on the other side of the position space. 

0.14 
0.12 
0.1 

i'o.08 

ra 

2 0.06 
0.04 
0.02 

Particle position 

Figure 2.1: Spread of probability distribution of the Hadamard walk with the 
initial state |0) (dashed line) and initial state |1) (solid line) on the position 
space along with distribution of the classical random walk (dashed-circle) . The 

distribution is for 100 steps. 




To obtain left-right symmetry in the probability distribution. Figure 2.3(b), one 
needs to start the walk with the particle initially in a symmetric superposition 

state 

|^,„,) = -^(|0)+z|l))®|^o). (2.30) 
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For Hadamard walk on the hne it is shown that after t steps, the probabihty 
distribution is spread over the interval and reduce quickly outside this 

region [17]. The moments have been calculated for asymptotically large number 
of steps t and the variance is shown to vary as 



Hadamard walk has been extensively studied and the symmetric distribution de- 
pends on the initial state of the particle (see [17, 18, 67, 89]). 

In the following section we will present a quantum walk using an operation from 
SU(2) group as a generalized quantum coin. The parameters in the SU(2) group 
will give control over the dynamics and amplitude distribution of the quantum 
walk. This eliminates the role of initial state of the particle to obtain a symmetric 
amplitude distribution. 

2.3.3 Optimization using SU(2) coin 

In the introduction to this section we defined a discrete-time quantum walk using 
B GU(2) (2.12) as quantum coin operation. To control and optimize the dynamics 
of the discrete-time quantum walk we consider the three parameter operator from 
the group SU(2) rather than the complete U(2) operator. Apart from the global 
phase which does not affect the amplitude distribution, the effect of U(2) operator 
on the quantum walk can be completely reproduced using SU(2) operator. 

We consider the group SU(2) which has for its elements, matrices of the form 



as the generalized quantum coin operation, a's and /3's are real values. For our 
present purpose the parametrization of SU(2) in terms of Euler angles turns out 




(2.31) 




(2.32) 
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to be the convenient one, 




e*^ cos(6') 
e~*^ sin(6') 



e*^ sin(6') 
e"*^ cos(6') 



) 



(2.33) 



We can get the Hadamard operator up to global phase by choosing C = C, = vr/2, 
6* = 7r/4 or by choosing ( = ^ = 0, 9 = 7t/4 and multiplying the operator by 
from left. 

For the analysis of generalized discrete-time quantum walk we consider the initial 
state of a particle to be in the symmetric superposition of basis states |\l/ms) (2.30). 
Implementing = S{B^^g^i^ ® 1) on l^^ms) evolves the particle to, 



The position probability distribution in (2.34) after first step corresponding to the 
left and right positions are |[1 ±sin(2^) sin(^ — (^)]. These probability distribution 
would be equal and lead to a left-right symmetry in position, if and only if ^ = 
(. That is, the parameters ^ 7^ C introduce asymmetry in the position space 
probability distribution from the first step itself. The effect of different values 
of coin parameters after 100 steps of walk is shown in Figure 2.2 plotted using 
numerically obtained values. We thus find that the generalized operator i?^,6»,c as 
a quantum coin can bias the probability distribution of the quantum walk in spite 
of the symmetry of initial state of the particle. This is not true in case of the 
Hadamard walk, where obtaining symmetric distribution depends largely on the 
initial state of the particle [17, 18, 67, 89]. This can be verified by substituting 
= ( = 0,6 = ir/A and multiplying the components of state |1) in (2.34) by —1. 

The walk beyond the first step involves interference and hence, it is appropriate 
to analyze the evolution after t steps using -B^,6»,c coin operator. The analysis 
shows how non-vanishing ^ and ( introduces bias. Therefore, the state after t steps 
can be written as 




n{9)) |0) ® 
(^))|l)®|^+i)]. 



(2.34) 



(2.35) 



m=—t 
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Figure 2.2: Spread of probability distribution for different values of 9, 
C using quantum coin operator -B^.e.c;- Parameter ^ shifts the distribution to 
the left: (a) = (f,f,0) and (c) =' (if,f,0). Parameter C shifts it to the 
right: (b) = (0, f , f ) and (d) = (0, f,ff). The initial state of the particle 
l^ins) = ~75(|0) + ^ iV'o) and the distribution is for 100 steps. 



and proceeds according to the iterative relations, 

Am,t = e'^ cos{e)Am+i,t-i + e'^ sm{e)Bra+i,t^i (2.36a) 
Bm,t = e-'« cos(^)i3„,_i,t_i - e-^^ sin(^)^^_i,4_i. (2.36b) 

After (t+1)^^ step, Am,t and Bm,t shifts to position (m— 1) and (m+1) respectively. 
A little algebra reveals that the solutions to Am^n and B^^n (2.36) can be decoupled 
and shown to satisfy the following equations : 

Am,t = cos{e) {e'^Am+i,t-i + e-^^^™_i,f_i) - A,n,t^2 (2.37a) 
B^,t = cos(^) (e*«S^+i,i_i + e-*«S^_i,t_i) - K^,t_2. (2.37b) 

The decoupled expression (2.37) shows dependence of amplitude on the values 
from previous two steps. By repeatedly substituting the dependence relation for 
amplitudes on the right hand side of (2.37) until amplitude term at t = appears, 
reveals the dependency of Am,t and Bm,t on amplitudes of the initial state and coin 
operation (2.34). We know that to obtain a spatial symmetry of the probability 
distribution from a particle initially in symmetric superposition state, the walk 
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should be invariant under an exchange of |0) ^ |1), and hence should evolve Am,t 
and Bm^t alike (as in the Hadamard walk [90]). From the above analysis we see 
that Am,t and Bm,t are symmetric to each other and evolve alike for all value of 6 
only when ^ = C- Positive ( contributes to constructive interference towards right 
and destructive interference to the left, and vice versa in case of C,, as shown in 
Figure 2.2. The inverse effect can be noticed when the and ( are negative. 

One can show that the coin operator 



Setting ^ 7^ C introduce asymmetry, biasing the walk. Therefore, a quantum coin 
operation 5^,0,^ when ^ ^ ( can be called as a biased quantum coin operation on 
a particle initially in symmetric superposition state^. In Figure 2.2 we show the 
biasing effect for (^, 9, Q = (0°, 60°, 75°) and for (75°, 60°, 0°). 

From (2.37) we know that only when ^ = C the quantum walk on a particle initially 
in symmetric superposition state evolve symmetrically in position space. From 
Section 2.3.2 we also know that the Hadamard walk on particle initially in basis 
state |0) biases the distribution to the left and in state |1) biases the distribution 
to the right. Therefore, for a particle initially in symmetric superposition state, by 
varying the parameter ^ and ( the results obtained for walk starting with one of the 
basis |0) or |1) (or any other nonsymmetric superposition) state can be reproduced. 
Similarly, a particle initially in one of its basis states |0) or |1) can also be evolved 
to obtain symmetric probability distribution by making an appropriate choice of 
the parameter ^ and (. For example, quantum walk on a particle initially in state 
|0) with coin parameters 6 being any value, C, = and ( = 7r/2 and for a particle 
initially in state ,^ = 7r/2 and C = return a symmetric probability distribution 
in the position space (symmetric distribution similar for Figure 2.3). Therefore 
quantum coin operation B^^g ^^ and initial state of the particle complement each 
other. The variance can be varied using the parameter 6. Compared to the 
Hadamard walk or a walk using single variable parameter unitary coin [18], a 

"'^It should be noted that for a particle initially in one of its basis state |0) or the coin 
operation i3o,e,o itself biases the walk, state |0) to the left and |1) to the right. When the initial 
state is a symmetric superposition of the basis state, the biasing on state |0) is compensated by 
biasing on state |1) leaving the effective walk distribution symmetric. 
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walk using a SU(2) operator as a quantum coin gives better access to control the 
dynamics of discrete-time quantum walk irrespective of the initial state of particle. 

We note that using of three parameter operator U(2) ^ SU(2) of the form 

_ / e'^cos{9) e'^sm{9) 
~ y e-^« sin(^) -e-'^ cos{e) 

in place of -Bg,^,^ will also produce the same effect on amplitude distribution during 
the quantum walk evolution. That is, 

S{B^,e,C ® l)|^*n.) = S{Bls^^ ® l)\^^ns) (2.40) 

when l^^ins) is (2.30). The operator B'^^^ reproduces the Hadamard operation for 
^ = ( = and 6 = n/A (whereas -85,0, c needs an additional Pauli z operation, 
cr^). For the convenience of reproducibility of properties of Hadamard walk, in 
this thesis we will switch between these two operators since it does not alter any 
of the results. We justify this claim as the existence of symmetries in quantum 
walk in Chapter 3. 
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Figure 2.3: Spread of probability distribution for different values of 9 using 
the quantum coin operator Bq q q. The distribution is wider for (a) = (0, y^,0) 
than for (b) = (0, ^, 0) and (c) = (0, 0), showing the decrease in spread with 
increase in 9. The initial state of the particle |^ins) = (|0) + i|l)) \ipo) and 

the distribution is for 100 steps. 




(a) = (0°, 15°, 0°) 

(b) = (0°, 45°, 0°) 

- (c) = (0°, 75°, 0°) 
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Figure 2.4: A comparison of variation of a with 6 for different number of steps 
of walk using operator -Bo,e,o using numerical integration. 
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Figure 2.5: Value of C for a given 9 with increase in number of steps. C 
remains constant with increase in number of steps and decreases with increase 

in e. 



Approximating the variance: 

The effect of parameter 6 in coin operation -B5,6»,c on the distribution when ^ = 
C = is shown in Figure 2.3, obtained by numericahy evolving the density matrix. 
The numerical results in Figure 2.4 shows decrease in variance with increase in 
the value of 6 for different number of steps of quantum walk. 
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The change in variance for different values of 6 is attributed to change in the value 
of Ce, a constant for a given 6, Figure 2.5. The dependence can be written in the 
form 

= Cet". (2.41) 

Therefore, starting from Hadamard walk (6' = |; ^ = C = 0), the variance can be 
increased {6 < |) or decreased (6 > |) using B^^e,c,- Figure 2.6 shows the depen- 
dence of C on 9. To analyze the effect of 9 on variance further, it is instructive 
to first consider the extreme values of parameters in B^qx^. If ^ = ^ = = 0, 
-Bo,o,o = o"2, the Pauli z operation, then W^^^ = S and the two superposition 
states, |0) and move away from each other without any diffusion and inter- 
ference having high o"^ = t^. On the other hand, if 6* = |, then i?o,f ,o = (^x, the 
Pauli X operation, and the two states |0) and |1) cross each other going back and 
forth, thereby remaining close to position j = and hence giving very low o"^ ~ 0. 
These two extreme cases are not dominated by interference effect, but they define 
the limits of the variance. Intermediate values of 9 between these two extremes 
show intermediate drifts and quantum interference and after t steps, the proba- 
bility distribution is spread over the interval (— t cos(6'), t cos(6')) [18]. Probability 
distribution is almost beyond |tcos(^)| and for all practical purpose it can be 
considered to be equal to 0. This is also verified by analyzing the distribution 
obtained using numerical integration technique. 

Numerically obtained data of variance for different 9 overlaps with (1 — sin(6'))t^. 
That is, the best fit function of 9 for the numerical data can be written in the 
form 

= Cet^ ^ [1 - sm{9)]t\ (2.42) 

To arrive at an expression for variance as a function of 9, a function of 9 for 
the probability distribution of quantum walk can be used. A function of 9 for 
probability distribution can be approximated by fitting a function that envelop 
the probability distribution obtained for different values of 9 from the numerical 
integration technique. Figure 2.7 ^. 

^Arriving at 9 dependent function from first principles, that is from the analytics of the 
dynamics of quantum walk is an ideal method. For this thesis, I have not considered that 
method due to constrains from the complexity of the dynamics of quantum walk. 
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Figure 2.6: Plot of the variation of Ce when rj = |^ — C| =0° from numerical 
data and the function (1 — sin(^)) to which it fits. The effect of maximum 
biasing, rj = 90° on Ce is also shown and it is very small. 



One of the 6 dependent function that envelop the probabiUty distribution is, 



where, K{e) = ^ cos(6')[l + cos2(26')][l + sin(6')]3. The above probability distribu- 
tion P{j) as a function of 6 was obtained by trying out different function of 6 until 
a reasonable fit for the numerically obtained distribution for quantum walk was 
obtained. Figure 2.7 shows the probability distribution obtained by using (2.43). 

The probability distribution for any 6 is spread over the position j in the interval 
(— t cos(^), t cos(6')). Therefore, for any 6 the position j can be parametrized by 
introducing a function of 0, /(0). That is. 



■^Since K{9) is an approximate function, other forms of K{6) which closely fit the envelop can 
also be considered. 




(2.43) 



j ^ /(0) =tcos{9) sin(0) 



(2.44) 
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Figure 2.7: The approximate probability distribution that envelop the nu- 
merically obtained probability distribution is obtained using (2.43) for different 
value of 9. The distribution is for 100 steps. 



where range from — ^ to |. For a walk with coin -Bo,6i.O; the mean of the distri- 
bution is zero and hence variance can be analytically obtained by evaluating 



tcos(e) 

a' = I P{j)fdj ^ I P{j)fdj = I " P{f{<P)){f{<j,)rf{<j,)d<j,. (2.45) 

J J-tcos(e) 

Substituting appropriate values and simplifying we get 



a 



(^ + ^"^(2^))e^(^)(-'--W^O(tcos(^)sin(0))^ 

X {t cos{9) cos((/))) # = ^2(1 - sm{e)), (2.46) 



That is, 

= Cet" ^{l- sm{e))t^. (2.47) 

The function matches the result Ce = {1 — sin(6')) obtained from the numerical 
data, as shown in Figure 2.6. Intermediate steps are presented in the Appendix 
B. 

We note that biasing the walk by setting ^ 7^ in -8^,6/,^ does not alter the width of 
the distribution in the position space but the probabihty decreases as a function 
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Figure 2.8: Variation of entropy of measurement H(j) with 9 for different 
number of steps of quantum walk. The decrease in H{j) is not drastic tih 9 is 



close to or I . 



of cos(?7) on one side and increases as a function of sin(?7) on the other side, where 
7] = \^ — C\- The mean value j of the distribution, which is zero for -Bo,e,0) attains 
finite value with non- vanishing rj, that contributes for an additional term in (2.45), 

/t cos (6) 
P{j){j-jydj. (2.48) 
'tcos(e) 

The contribution of r/ is a very small decrease in the variance of biased quantum 
walk, see Figure 2.6. 



2.3.3.1 Entropy of measurement 

As an alternative measure of fluctuation in position, we consider the Shannon 
entropy of the quantum walk. Shannon entropy of the walk is obtained using 
position probability distribution Pj obtained by tracing over the coin basis, 

i7(j) = -5^P,logP,. (2.49) 

j 

The quantum walk with a Hadamard coin toss, -Bo,f ,0; has maximum uncertainty 
associated with the probability distribution and hence measurement entropy is 
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maximum. For ^ = C = and low 6, operator i?o,e,o is almost a Pauli z operational, 
leading to localization of walker at ±t. At 6 close to |, with ^ = = 0, 
approaches the Pauli x operation ax, leading to localization close to the origin, 
and again, low entropy. However, as 9 approaches |, the splitting of amplitude 
in the position space increases towards the maximum. The resulting enhanced 
diffusion is reflected in the relatively large entropy at ^, as seen in Figure 2.8. 
Figure 2.8 is the measurement entropy with variation of 6 in the coin -Bo,6i,o ioT 
different number of steps of quantum walk. The decrease in entropy from the 
maximum by changing 6 on either side of ^ is not drastic until 6 is close to or 
|. Therefore for all practical purposes, small entropy can be compensated for by 
the relatively large Cq, and hence a^. The effect of ^ and ( on the measurement 
entropy is of very small magnitude. These parameters do not affect the spread of 
distribution and variation in height reduces the entropy by a very small fraction. 
In other cases, such as mixing of quantum walk on an n— cycle briefly discuss in 
Section 2.3.3.2, it is ideal to use a lower value of 6. 

2.3.3.2 Mixing time on an n-cycle 

The simplest finite Cayley graph is an n-cycle with n vertices in closed path 
[91]. Due to closed structure of the ra— cycle graph, with n positions in Tip, the 
conditional shift operation S (2.14) for a discrete-time quantum walk on a line 
takes the form 

n— 1 n— 1 

= |0)(0| ® Yl n){^Pj\ + |1)(1| ® 5Z '''^J'+l n)(^,|. (2.50) 

j=0 j=0 

The quantum state after t steps of discrete-time quantum walk on particle initially 
in state |\l/jn) (2.11) is written as 

j=0 

where | \l/j^t) is the state at position j after t steps of quantum walk. 

Mixing time: Mixing time is the time it takes for the probability distribution 
on a graph (position space) to converge to a stationary distribution. The classical 
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Figure 2.9: A comparison of mixing time M of the time-averaged probabihty 
distribution of a quantum walk on an n-cycle for different value of 9 using coin 
operation Bq^q q, where n, the number of position, is 101. Mixing is faster for 
lower value of 9. The distribution is for 200 cycles. 

random walk approaches a stationary distribution independent of its initial state 
on a finite graph. Unitary (i.e., non- noisy) quantum walk, does not converge to any 
stationary distribution. But by defining a time-averaged probability distribution, 

mT) = ^X]^(^'^)' (2.52) 

r=0 

obtained by uniformly picking a random time r between and (T — 1), evolving 
for t time steps and measure to find a particle at a given vertex (position j), 
a convergence in the probability distribution can be seen even in the quantum 
case. It has been shown that for the quantum walk on an n-cycle, the mixing 
time is bounded above by M = O(nlogn), almost quadratically faster than the 
classical case which is O(n^) [91]. From previous section, we know that a quantum 
walk can be optimized for maximum variance and wide spread in position space, 
between (— tcos(6'), tcos(6')) after t steps. For a walk on an n-cycle, choosing 6 
slightly above would give the maximum spread over n vertices in the cycle for 
n/2 steps of the quantum walk. Maximum spread during each cycle distributes 
the probability over the cycle faster and this would optimize the mixing time. For 
optimal mixing time, it turns out to be ideal to fix ^ = in -Bg,e,C; since biasing 
impairs a proper mixing. Figure 2.9 is the time averaged probability distribution 
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of a quantum walk on an n-cycle graph after time t = nlogn, where n is 101. It 
can be seen that the variation of probabihty distribution over the position space 
is least for 6 = 15° compared to 6 = 45° and 6 = 75°. 

2.3.4 Randomizing coin operations from SU(2) group 

In Section 2.3.2 and 2.3.3 we discussed quantum walk using Hadamard operation 
and SU(2) operation as quantum coins. Though we used a generalized SU(2) coin 
i?^,6»,c (2-33), identical coin parameters were used during each step of the walk and 
we also saw the dependence of variance on the coin parameter 6. 
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Figure 2.10: Quantum walk using randomly picked coin operation for each 
step of walk from SU(2) group with randomly picked £,,9,( ^ {0; "^/^j for -B^,^,^ 
and -Bo,6»,o respectively. We see that the variance of the distribution is very much 
close to the variance of the classical random walk (CRW) distribution. 

Evolution of quantum walk using randomized coin operation can be constructed 
by randomly chosen quantum coin operator for each step from a set of operators 
in SU(2) group. That is, 

S{B^,eu. ® 1) ■ ■ ■ S{B^M ® 1) ■ ■ ■ S{B^M ® (2.53) 

with randomly chosen parameters ^,^,C ^ {0,vr/2} for each step. Though the 
coin parameters are randomly chosen for each step, the evolution is unitary and 
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involves interference of amplitudes, and the effect is seen in probability distribu- 
tion, see Figure 2.10. From the numerical evolution we also see that variance of the 
distribution is much closer to variance of the classical random walk distribution. 
However by restricting the range of the coin parameters that can be used for the 
walk, the probability distribution can be localized or made to spread wide in posi- 
tion space. One simple example we can consider is, by randomly picking different 
6 for each step from a subset of complete range of 6, subset with 6i G {0, 7r/4} and 
6*2 G {7r/4, vr/2} for B^^q^ q and B^ q^^c^ respectively, where other parameters are still 
picked from the complete range, ^, C ^ {0, 7r/2}. Probability distribution obtained 
is shown in Figure 2.11. Walk using 6i G {0, 7r/4} spreads wider in position space, 
whereas the walk using 62 G {7r/4,7r/2} localizes the distribution. 
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Figure 2.11: Quantum walk using two different set of randomly picked param- 
eters G {0,7r/2},6'i G {0,7r/4} and 62 G {7r/4,7r/2} for Bf^^e^^c. and -Bg,e2,C 
respectively. The distribution is for 200 steps, walk using 61 G {0,7r/4} spreads 
wider in position space, whereas the walk using 62 G {7r/4, 7r/2} localizes the 

distribution. 



2.4 Recurrence of quantum walk 

In the dynamics of physical system, from a single free particle to stellar dynamics, 
understanding of recurrence phenomenon have significantly contributed to a better 
understanding of the system dynamics [20]. For a classical conservative system, 
whether discrete or continuous in time, the Poincare recurrence theorem states 
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that any phase-space configuration of a system enclosed in a finite volume will 
be repeated as accurately as one wishes after a finite interval of time (with no 
restriction on the interval) [68, 69]. A similar recurrence theorem is shown to hold 
in quantum theory as well [70, 71]. In a system with a discrete energy eigenvalue 
spectrum {En}', if "^{to) is its state vector at the time to and e is any positive 
number, there exists a finite time T such that, 

|*(T)-^(to)| <e. (2.54) 



For classical random walk we can consider the recurrence probability Po(t), i.e., 
the probability of periodicity of dynamics that the particle returns to origin during 
the time evolution (t steps). It is characterized by the Polya number 

'Pcrw = 1 Y^oo T-) / , \ • (2.55) 

If the Polya number equals one, the classical random walk is recurrent, otherwise 
the walk is transient, i.e., with nonzero probability the particle never returns to 
the origin. For a classical random walk to be transient the series X]t^o-^o(^) ^^st 
converge [92]. Polya proved that the one- and two- dimensional classical random 
walk are recurrent [93], while in higher dimension for each dimension a unique 
Polya number is associated and the classical random walk is transient. 

In standard quantum mechanics, initially localized wave packet in state |\E'(to)) 
which can spread significantly in a closed system can also reform later in the 
form of a quantum revival, i.e., the spreading reverses itself and the wave packet 
relocalizes [72]. The relocalized wave packet can again spread and the periodicity 
in the dynamics can be seen validating the quantum recurrence theorem. The time 
evolution of a state |^E'(t)) over entire position j, or a wavefunction \E'(j, t) at a 
position j after time t is given by a deterministic unitary transformation associated 
with the Hamiltonian. During quantum evolution we deal with amplitudes and 
probability density P(j, t) = |\l/(j, t)P at position j after time t appears only when 
we collapse the wave packet to perform measurement. 
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Unlike standard wave packet spreading in free space and harmonic potential^, the 
quantum walk spreads the wave packet in multiple possible paths with amplitudes 
corresponding to different paths interfering. In this section we show that due 
to particle-position entanglement and interference effect during the evolution of 
the quantum walk, the wave packet delocalizes over the position space as a small 
copies of the initial wave packet. These delocalized copies of fractional wave packet 
fails to satisfy the complete quantum recurrence theorem during quantum walk 
evolution. That is for quantum walk, the recurrence we considered is return of 
the unit probability amplitude at the origin during the dynamics. However, due 
to revival of fractional wave packets, a fractional recurrence can be seen in the 
quantum walk. The probabilistic characterization of the recurrence in quantum 
walk has been done in [17] and the characterization using quantum Polya number 
defined in [94] can be used to show the fractional recurrence nature of quantum 
walk. In Section 2.4.2.2, we also show the exceptional cases of quantum walk that 
can be constructed by suppressing or minimizing the interference effect and get 
closer to complete recurrence. 

2.4.1 Quantum recurrence theorem 

Quantum recurrence theorem in the dynamics of a closed system states that there 
exist a time T when 



where \I/(T) = is the state of the system after time T, \l/(to) = l^o) is the 
initial state of the system and e is any positive number [70]. 

The recurrence of the complete state of the system or exact revival happens when 
all the expectation values of observables Ao of the two states |\1/t) and |\l/o) are 
equal to one another, that is. 



''^During standard wave packet spreading in free space and harmonic potential, a wave packet 
which is initially Gaussian, spreads retaining the Gaussian shape causing increase in the full- 
width at half maxima (FWHM). 




(2.56) 



(2.57) 
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In classical dynamics, the characterization of the nature of recurrence can be con- 
veniently done using probabilistic measures. Measurements on a quantum system 
leaves the state in one of its basis states with certain probability. Therefore, re- 
currence in quantum systems can be analyzed using comparative evolution of the 
two identically prepared quantum system with the initial states |\&o,o) at position 
j = and time t = 0. We will consider two cases of comparative analysis. 

Case 1: Consider two identically prepared particle wave packets which revive 
completely in the position space at time T. One of the two particle wave packets 
at position j = and time t = is first evolved to spread in position space and 
then reverse the spreading till it relocalizes completely at position j = at time 
T. The measurement performed on this particle will collapse the wave packet at 
the relocalized position with an expectation value 

(*o,t|^pI*o,t) = (^o,o|^p|*o,o) = 1, (2.58) 

where Xp is the position operator. After the measurement at time T, the system 
is further evolved for an additional unit time and the corresponding state can be 
given by I'^j^Tm+i)- The subscripts j and (Tm + 1) stand for position and time 
(T + 1) with intermediate measurement being performed at time T. The second 
particle wave packet at position j = and time t = is evolved up to time (T + 1) 
directly without any measurement being performed at time T and the state can be 
written as |\E'j,T+i)- Since both the wave packets completely relocalize at position 
j = after the evolution for time period T, irrespective of the measurement being 
performed, expectation value for both the particles after time (T + 1) would be 
identical, 

{^j,n,+l\Xp\^,,TM+l) = (^',,T+l|Xp|^',-T+l), (2.59) 

with j spanning over all position space. 

Case 2: Consider two identically prepared particle wave packets which does not 
relocalize completely at position j = at time T, i.e., revive fractionally or does 
not revive at all. The measurement will collapse the wave packet and return the 
expectation value 

{^o,t\Xp\^o,t) = S ; 0<5<1. (2.60) 
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In this case the two identically prepared particle evolved to time (T+ 1), one with 
a measurement being performed at time T and an other without any measurement 
being performed. With very high probability the measurement will not return the 
same expectation value, that is 

(^,-Tm+i|^pI^.,Tm+i) ^ {^j,T+l\Xp\^,,T+l, ) (2.61) 

with j spanning over all position space. Non-zero values in both expectation values 
in the inequality can act as a signature of fractional recurrence of the quantum 
state at time T and a zero expectation value on the left hand side shows transient 
dynamics. 

From the analysis of the above two cases we can conclude that, if a system is 
completely recurrent at time T the two states |\Efj_TM+i) |\Efj_T+i), one with 
intermediate measurement at time T and other without any measurement at time 
T will be equal to one another. This will be helpful for us in understanding the 
fractional recurrence nature of the quantum walk. 

2.4.2 Fractional recurrence of quantum walk 
2.4.2.1 On a line 

The state of the particle wave packet after implementing the discrete-time quantum 
walk of t steps on a line with unit time required to implement each step can be 
written as 

l**) = W^|Acl*o,o) = 5^|*,,t). (2.62) 

j 

Where is the state after time t over entire position space, I'^j^t) is the state 
of the delocalized wave packet at each position j in Tip and I^E'o.o) is the coin and 
position state of the wave packet before implementing the quantum walk. In the 
quantum walk process that involves a deterministic unitary evolution, the par- 
ticle wave packet delocalize and spread over the position space forming a mini 
wave packet. During this delocalization process the mini wave packets interfere 
and entangle the position and coin Hilbert space, Tip and Tic- The interference 
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Figure 2.12: The plot of [1 - Po{t)], where Po{t) is the probabihty of particle 
at the origin with t being the number of steps in the discrete-time quantum walk 
evolution. The plot for quantum walk on a line using different coin operation 
parameter 6 = 15°, 45°, and 75° is shown. With an increase in 6, the quan- 
tum Polya number, which characterizes the fractional recurrence nature of the 
quantum walk, also increases. The plot with = 0° is for a walk on an n— cycle 
with a completely suppressed interference effect, where n = 51. For a walk to 
be completely recurrent there should exist a t = T where [1 — Pq{T)] = 0. Any 
non-zero value can be used to characterize the fractional recurrence nature of 
the quantum walk and an absolute 1 for all t would show that the quantum 
walk is completely transient, (b) Close up of the plot. 



and the entanglement between the Tip and "He during the standard quantum walk 
evolution does not permit complete relocalization of the wave packet at initial 
position after any given number of steps t. Therefore the argument leading to 
(2.60) and (2.61) holds to show that the complete recurrence of the quantum state 
does not occur during the evolution of the quantum walk process on a line. That is, 



In a discrete-time quantum walk evolution on a line and on an n— cycle dominated 
by the interference of quantum amplitude ^, there exists no time T when the quan- 
tum state of the system revive completely and repeat the delocalization and revival 
at regular interval of time. 



^By choosing extreme value of 9 (0 or 7r/2) in the quantum coin operation, quantum walk on 
a line can be evolved without constructive or destructive interference taking place. 
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The above statement can also be quantified in the following way: After imple- 
menting a quantum walk, the wavefunction describing the particle at position j 
and time t can be written as a two component vector of amplitudes of particle 
being at position j at time t with left and right propagating component 

l*-)-f!f!!)^ (2.63) 

Lets analyze the dynamics of discrete-time quantum walk on wave packet '^{j,t) 
driven by single parameter quantum coin 

/ cos(^) sin(^) \ 

Be=\ ^ ' ^ M 2.64 

V sin(^) -cos{e) J 

and shift operator S (W = S{Bg ® 1)). In terms of left and right propagating 
component it is given by 

^L{j,t + 1) = cos{e)^R{j -l,t) + sm{e)^ l{j -I, t) 

^R{j,t + 1) = sm{e)^ii{j + l,t)-cos{e)^L{j + l,t). (2.65) 



Then the probability of being at position j and t is 

P(j,t) = |vl>^(j,t)|2 + |vl/^(j,t)|2 (2.66) 

and sum of probability over the entire position space is 

^P(j,t) = l. (2.67) 
j 

After time t with unit time required to implement each step of the quantum walk 
on a line, the wavefunction will be spread between j = —t to +t and can be written 
over position space as 

t t 

J2 ^U^) = E [*^(-^'^) + '^RO,t)] . (2.68) 
j=-t j=-t 
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It should be noted that for even number of steps the amphtude at odd labeled 
positions is and for odd number of steps amplitude at even labeled position is 0. 

For classical random walk, each step of walk is associated with the randomness 
and the probability of the entire particle therefore, however small the probability 
is at the origin, it is attributed to the recurrence of the entire particle. That is, 
< Pcri«(0,T) < 1 at some time T shows the recurrence of complete particle 
at the origin with some finite probability. For the discrete-time quantum walk 
evolution described by (2.65) we found that the information of the coin degree of 
freedom is carried over to the later steps during the dynamics of the walk making it 
reversible (2.37). The randomness and the probability in quantum walk comes into 
consideration only when the wave packet is collapsed to discard the signature of the 
quantum coin degree of freedom from earlier steps. Therefore, for a discrete-time 
quantum walk to be completely recurrent, the condition 



has to be satisfied for some time T. Any P{0,t) < 1 shows that the particle is 
present in superposition of position space (origin and other positions), which we 
will call as the fractional recurrence nature of the quantum walk. From (2.65) and 
(2.66), we can conclude that (2.69) is satisfied only when 6 = 7r/2 and t is even, 
that is when there is no interference of the quantum amplitudes. For ^ = the 
two left and right component move in opposite directions without returning and 
for any 0<6'<7r/2we get P(0,t) < 1 showing the fractional recurrence of the 
quantum walk. 

Alternatively, using quantum Fourier analysis to study the evolution of the discrete- 
time quantum walk on a line, it is shown that the amplitude at the origin and to 
a very good approximation for any position j between two dominating peaks (for 
example see Figure 2.3) decreases by 0{1/\^)[17]. Therefore, 



and their exists no time T = t where the walk is completely recurrent showing the 
fractional recurrence nature of the quantum walk. 



P(0,T) = |^i(0,T)|2 + |^«(0,T)|2 = l 



(2.69) 




(2.70) 
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A probability based characterization of the recurrence nature of the quantum 
walk, quantum Polya number was defined for an ensemble of identically prepared 
quantum walk systems by the expression 

oo 

= 1 - n [1 - ^o(^)] ' (2.71) 
t=l 

where Po{t) is the recurrence probability of the particle. Each identically prepared 
particle is subjected to different number of steps of quantum walk from 1 to oo 
and the probability of the particle at the origin is measured and discarded [94]. 
The probability that the particle is found at the origin in a single series of such 
measurement records is the quantum Polya number. The quantum Polya number 
was calculated for various coined quantum walks and it is shown that in the higher 
dimension it depends both on the initial state and the parameters of the coin 
operator whereas, for the classical random walk the Polya number is uniquely 
determined by its dimensionality [95]. 

To show the quantum walk to be completely recurrent adopting a probability 
based characterization given by (2.71) needs to have at least one of the many 
particles, each evolved to different steps t from 1 to oo, to return Po(t) = 1. 
If < Po{t) < 1 for any t, then only with certain probability Po{t) the wave 
packet collapses at j = 0, that is, prior to measurement the particle existed in 
superposition of position space. If < Po{t) < 1 for all the particles evolved to 
different steps from 1 to oo then the fractional recurrence of the quantum walk is 
characterized by the Vqw 

In Figure 2.12 the plot of [1 — Po{t)] is shown for a discrete-time quantum walk 
on a line where the different coin operation parameter 6 = 15°, 45°, 75°. With 
increase in 6 the quantum Polya number, which can also be called as the fractional 
recurrence number Vqw, also increases. 

2.4.2.2 On an n— cycle 

We have discussed the dynamics of discrete-time quantum walk on an n— cycle and 
the mixing time using coin degree of freedom in Section 2.3.3.2. In this section we 
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Figure 2.13: Probability at the initial position after different number of steps 
of quantum walk on 15— and 24— cycle. This clearly shows that even for steps 
as large as 5000, their is no signature of complete recurrence. The distribution 
is obtained using Hadamard operation, -Bo,7r/4,o &s quantum coin operation. 



will discuss the recurrence nature of quantum walk on an n— cycle. 

As discussed earlier, in the case of classical random walk walk on a line a non- 
zero probability is sufficient to show the recurrence of the entire particle and a 
uniform distribution returns a non-zero probability at the origin. Whereas uniform 
distribution in the quantum case does not reveal the complete recurrence nature 
of the particle as it does in the classical case. 

In Figure 2.13 we show numerically that the probability of finding the particle at 
the initial position j = after different number of steps of quantum walk on 15— 
and 24— cycle. The distribution is obtained for up to as large as 5000 steps. At 
no time t = number of steps, a unit probability value is returned at the initial 
position showing that the wave packet evolved using quantum walk on an n— cycle 
fails to revive completely at the initial position. 

The failure of the wave packet to completely revive and recur at initial position 
can be attributed to the interference effect caused by the mixing of the left and 
right propagating components of the amplitude. By suppressing the interference 
effect during the evolution in a closed path one can get closer to the complete 
relocalization, revival at initial position j = 0. For example, the wave packet 
can be completely relocalized at j = on an n— cycle at t = n and make the 
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Figure 2.14: Probability at the initial position on 8— cycle, (a) The distribu- 
tion is obtained using Hadamard operation i?o,7r/4,0) due to return of amplitudes 
to initial position (constructive interference at the origin) before interference of 
amplitude dominates uniformly over the entire vertices, recurrence is seen, (b) 
The use of quantum coin, So,7r/6,o during the evolution does not lead to com- 
plete constructive interference at origin and hence complete recurrence is not 

seen. 



quantum walk recurrent by choosing an extreme value of coin parameters ^, C) = 
(0°,0°,0°). By choosing 6 = 6, 6 being very small and close to 0° during the 
evolution, the interference effect is minimized and will return a near complete 
recurrence on an n— cycle. In Figure 2.12, the plot of [1 — Po(t)] is shown for 
9 = 0° and n = 50. The interference effect is completely suppressed and the 
quantum walk recurs after every 50 steps. 

From the numerically data for upto 5000 steps of quantum walk we note that for 
small n— , especially when n is even, the left and right propagating amplitude re- 
turn back completely to the initial position (origin) before the mixing and repeated 
interference of the amplitudes takes over at non-initial position in the evolution 
process. Therefore, for a quantum walk on an ra— cycle with n being even up to 8 
a complete revival and recurrence of wave packet at initial position j = is seen. 
Figure 2.14 [67]. For a quantum walk on particle initially in symmetric super- 
position state I'^ins) with coin operation i?o,6»,o and n > 8, due to larger position 
Hilbert space the interference effect at the non-initial position dominates reduc- 
ing the recurrence nature of the dynamics. In Figure 2.15, for quantum walk on 
10— cycle, a small deviation from complete recurrence is shown. 



47 



Chapter 2. Quantum walks 



I 0-8 

"en 
o 

Q. 

0.( 



0.' 

CO 

n 
o 



0.2 



,,l 



-10-cycle 



50 100 150 200 250 
Number of steps 



o 0.8 




1000 2000 3000 4000 
Number of steps 



5000 



(b) 



:|o.99 
o 

CL 

~ 0.98 



0.97 



0.96" 



-10-cycle 



1000 2000 3000 4000 
Number of steps 



5000 



(c) 

Figure 2.15: Probability at the initial position on 10— cycle. The distribution 



is obtained using Hadamard operation B^ 



0,7r/4,0- 



A small deviation from com- 



plete recurrence can be seen, (a) and (b) are probability at initial position for 
quantum walk up to 250 steps and 5000 steps respectively and (c) is the close 
up of the probability and we note that the probability is not exact 1 at any time 

within 5000 steps. 



For example, we will consider a small, odd n = 5. If the positions on a cycle 
are marked as j = 0, 1,2,3,4 after the third step of the quantum walk, the left 
propagating amplitude move from position 2 to 3 and the right propagating ampli- 
tude move from 3 to 2. That is, after the second and third step of quantum walk 
using shift operator of the form S'^ for a walk on an n— cycle (3.1) and Hadamard 
operation H as coin operation on a particle initially in superposition state 
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Number of steps Number of steps 

(a) (b) 

Figure 2.16: Probability at the initial position after different number of steps 
of quantum walk on 5— cycle, (a) is the complete plot whereas (b) is a close up 
of probability values between 0.9 and 1. At no step in the plot, probability is a 
unit value. This numerically shows that the quantum walk on an n— cycle does 

not recur for n=5. 



(2.30) takes the form, 

S'^iH ® l)S'\^,ns) = \{\0)^ + { |1) + i\0) )®|^o) - ® 1^2)} (2.72) 



S'iH ® 1)S%H ® l)51*i„,) = ^{|0) ® 1^2) + (|0) + |1) + 1\0)) ® 1^4) 

- (|1) - t\0) + ® 1^1) +t\l)(g) |^3)}.(2.73) 

The left and right propagating amphtude crossover without suppressing the mix- 
ing, therefore the constructive interference effect continues to exist even at posi- 
tions other than the origin during the evolution. In Figure 2.16 the probability 
of finding the particle at the initial position is shown. Due to the small size of 
the Hilbert space, the probability is seen to be close to unity but the closer look 
reveals that its only a fractional revival. 

We also note that localization effect found in 2D [96], in quantum walk using 
multi quantum coins to diminish the interference effect [97] or in walk using selec- 
tive randomization of coin operations (Section 2.3.4) can result in increasing the 
fractional recurrence number Vqw on a line and higher dimension. If the particle 
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wave packet is evolved in a position Hilbert space Hp with the edges that permits 
the wave packet to escape, the fractional recurrence nature of the quantum walk 
does not allow the quantum walk to be completely transient. Therefore, the frac- 
tional transient nature of the quantum walk is seen to complement the fractional 
recurrence nature. 

From this section we conclude that as long as the wave packet spread in position 
space interfering, forming mini wave packets during the evolution of the quan- 
tum walk process, it fails to satisfy the complete recurrence theorem. However, 
fractional recurrence can be seen. 
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2.5 Summary 



• Expression used to understand the continuous- and discrete-time quantum 
walk which are mathematically identical to Schrodinger and Dirac equations 
respectively are reviewed. Our simple decoupling analysis of the evolution 
shows the similarity between the discrete-time quantum walk evolution ex- 
pression and the Dirac equation. 

• Our construction of discrete-time quantum walk model using three param- 
eter quantum coin operation -B^,6»,c from SU(2) group optimizes the control 
over the quantum walk evolution. For example, parameter 6 to control the 
variance and maximize the variance for a walk on a line, control the mixing 
time on an n— cycle, and parameters ^ and ( to bias and control biasing in 
the walk. 

• We have shown that the quantum walk evolution dominated by the inter- 
ference of quantum amplitudes fails to satisfy complete recurrence theorem. 
However, fractional recurrence characterized by the quantum Polya number 
can be seen. 
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Symmetries and noise on 
quantum walk 

3.1 Introduction 

The probability distribution of discrete-time quantum walk on a line remains in- 
variant in position space when the operations to implement each step of the walk is 
augmented by certain operations. We refer to these discrete operations as symme- 
tries of the quantum walk. In this chapter we will study some of these symmetry 
operations. We further generalize the observations of these symmetries in the pres- 
ence of environmental effects, modeled by various noise channels such as, bit flip, 
phase-flip (decoherence without net dissipation), generalized amplitude damping 
(decoherence with dissipation) and squeezed generalized amplitude damping chan- 
nels [5, 64, 98, 99] on the coin space. We have found it convenient to explain the 
symmetries and effect of noise using quantum trajectories, and the numerical re- 
sults are obtained by numerical integration by evolving the density operator of the 
system. 

We extend these studies to quantum walk on an n— cycle. Interestingly, we find 
that the symmetry operations is sensitive to the walk topology, in the sense that 
the symmetry which holds for quantum walk on an one-dimensional line does not 
hold, in general, for a quantum walk on an n— cycle but leads to other interesting 
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behavior. The difference between the walk on the hne and an n— cycle can be 
attributed to the different ways the interference occurs between the propagating 
wavefunction. Quantum walk on the line involves interference between the for- 
ward and backward propagating wavefunctions whereas, the walk on an n— cycle 
involves interference between forward (backward) propagating waves from both the 
sides of the loop along with the interference between the forward and backward 
propagating wavefunctions. Noise on an ra— cycle tends to restore these symmetry 
both by classicalizing the walk and also desensitizing the symmetry operation as 
a topology probe for the quantum walk. 

These observations can have important implications for a better insight into, and 
for simplifying certain implementations of, quantum walks and are also of relevance 
to studies in quantum optics and condensed matter systems. Later in this chapter 
and in Chapter 4 we show that the application of these ideas can help simplify 
certain experimental implementations of quantum walk and can be used as an 
additional degree of freedom in applications of quantum walk in physical systems. 

This chapter is organized as follows. In Section 3.2, we discuss symmetries and 
noise on quantum walk on a line with Section 3.2.1 focusing on bit flip and phase 
flip symmetries and Section 3.2.2 focuses on the environmental effects. In Section 
3.3 we extend the studies to walk on an n— cycle and discuss the breakdown in 
symmetry. In Section 3.4 we discuss the experimental implications and conclude 
with summary in Section 3.5. 



3.2 Symmetry and noise operations on a line 
3.2.1 Bit flip and phase flip symmetries 

As defined in Section 2.3 each step of the discrete-time quantum walk consists of 
the quantum coin operation -Bg,^,^ (2.33) followed by the shift operation 

S = |0)(0|®5^|^,_i)(V',| + |l)(l|®^|^,+i)(^,| 

= |0)(0| ®a+ |1)(1| (3.1) 
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Here a and a? are unitary operators that are notationally reminiscent of annihi- 
lation and creation operations, respectively. Lets consider the application of the 
modified conditional shift operator of the form 



instead of S (3.1). Where X = 'i?> the Pauh x operator. Since S' = {X ® 1)5', 
i.e., it is equivalent to an application of bit flip following S*, conditioned on the 
internal state being |0) (|1)). That is, at any given position j, the particle will 
move to the left (right) and changes its internal state to |1) (|0)). Thus, 



5'(|0) ® 1^,)) = |1) ® and S'{\1) ® |^,)) = |0) ® (3.4) 



A relevant observation in this context is that there are physical systems where the 
implementation of 5" is easier than that of S. We will discuss one such system 
in Section 4.5.1 [57]. In that case, applying a compensatory bit flip on the inter- 
nal state, after each application of S", reduces the modified quantum walk to the 
usual scheme. In all, this would require {t — 1) compensatory bit flip operation in 
addition for a t step quantum walk, which adds to the complexity of the experi- 
mental realization. However, this additional complexity can be eliminated. For a 
quantum walk using B^^0x_ (2.33) as coin operation, applying a bit flip in each step 
can be shown to be equivalent to a spatial inversion of the position probability 
distribution. A quick way to see why bit flips are harmless is to note that they are 
also equivalent to relabeling the edges of the graph on which the quantum walk 
takes place, so that each end of each edge has the same label [100] ^ We may in 
this sense call a bit flip together with spatial inversion a symmetry of the quantum 
walk on a line. To be specific, when any unitary operation augmented during 
each step of quantum walk which may leave the position probability distribution 
unaffected, then it is called a quantum walk symmetry. 

^It is worth noting that in [55], bit flips arc employed to improve the praetieal implementation 
of a quantum walk on atoms in an optical lattice 



= (X®I)exp(-z(|0)(0|- |1)(1|)®P/) 



(3.2) 



S' = |l)(0|®5^|^,_i)(V;,| + |0)(l|®5^|V',+i)(^, 




1)(0| ®a+ |0)(1| 



(3.3) 
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Experimentally, the symmetries are useful in identifying variants of a given quan- 
tum walk protocol that are equivalent to it. This motivates us to look for other 
(discrete) symmetries of the quantum walk, which we study below. We begin with 
Theorem 1, where we note four discrete symmetries, associated with the matrices 
(f = 1, 2, 3, 4) (3.5), of the quantum walk. Thereafter two of these symmetries, 
B^^^ and B^'^\ are identified with operations that are relevant from the perspective 
of physical implementation. It is an interesting open question with relevance to 
practical implementation of quantum walks, whether other such symmetries of the 
quantum walk exist. 



Theorem 1. If B = fi^.e,^^ in (2.39) is replaced by any of 5^, B^^\ B^^\ or B^*\ 



the resulting position probability distribution of the quantum walk remains invari- 
ant. 

Proof. With the notation B = {&g,r} and B^^^ = {bql}, we find 



b'-^J = b,,re"^^ ■ = 5,,,e*'^<^ ; bf) = b.^rc''^^ ; b^^l = bg,re''^ , (3.6) 



where the matrix indices q,r take values and 1, z = +v^— 1, and the overbar 

denotes a NOT operation (0 ^ 1). The state vector obtained, after t steps, using 

^ Three parameter quantum coin operation (2.39) is used instead of SU(2) operation so that 
B reduces to Hadamard operation operation for ^ = = 0, = 7r/4 



given by. 




e*'' sin(6') 
_ei(<^-0 cos(^) 

eiinO sin(^) 
-e^^-^-^) cos(^) 





(3.5) 
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B and 5^^) as the coin rotation operations, are, respectively, given by 

l^i) = {SBY\a,l3)= b,,,q,_,---bg,,,,bq,,a\qu^ + 2Q -t), 

qi,q2,--- ,qt 

= Yl &.i,a(e^'^)'*-^'"-^'^+1?*,/3 + 2g-t), (3.7) 

qi,<i2,--- ,qt 

where Q = qi + ■ ■ ■ + qt- Consider one of the element |a, b) in the computational- 
and-position basis. Now, 

(a,6|^i) = e^''^a,6|vl>2), (3.8) 

where rj = qt_i + ■ ■ ■ + gi + a, which is fixed for given a and 6, and determined by 
b = f3 + 2Q — t and qt = a. As a result, 

|(a,fe|vl/i)r + |(a,fe|vl/i)|2 = |(a,6|vl/2)r + |(a,6|vl/2)r- (3-9) 

A similar proof of invariance of the position distribution can be demonstrated to 
hold when B is replaced by one of the other E^'s (f = 2,3,4). On account of 
the linearity of quantum mechanics, the invariance of the walk statistics under 
exchange of the B^^^^s and B holds even when the initial state \a, (3) is replaced 
by a general superposition or a mixed state. ■ 

Interchanging B and the B^^^^s may be considered as a discrete symmetry operation 
G : B ^ B* (where B* denotes any of the B^^^^s in (3.5)), that leaves the positional 
probability distribution invariant. We express this by the statement that 

W ~ GW, (3.10) 

where G refers to the application of G at each step of the walk, and W refers to the 
walk operation of evolving the initial state through t steps and then measuring in 
the position basis. Knowledge of this symmetry can help simplify practical quan- 
tum walks. Below we identify two of these quantum walk symmetry operations 
B -v^ B^^^ and B ^ B^'^\ associated with physical operations of interest. 
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We first consider the phase shift operation 



$(^)^|0)(0|+e^*|l)(l| 



(3.11) 



as a symmetry operation of a quantum walk. In our model, the quantum operation 
for each step is augmented by the insertion of $(0) just after the operation S{B^ 
1). At each step, the walk evolves according to 



1 

e''^ 

1 





1 




d + 



® d 




1 



® d 




1 

t 



® d'^ 



e*^ cos(6') e*'» sin(^) 



sin(^) 



cos(^) 



cos(^) 



sin(^) 



(3.12) 



This is equivalent to replacing B by B^^\ which, according to Theorem 1, leaves 
the walk distribution invariant. Thus the operation $(0), applied at each step, is 
a symmetry of the quantum walk. 

As a special case, the phase flip operation Z = a^, a Pauli z operation applied 
at each step, obtained by setting </> = vr, is a symmetry of the quantum walk. 
Representing the inclusion of operations $ or Z at each step of the walk by $ or 
Z, respectively, we express this symmetry by the statements : 



w ~ zw. 



(3.13a) 
(3.13b) 



Unlike the phase flip operation, bit flip is not a symmetry of the quantum walk on 
a line. However, the combined application of bit flip along with angular reflection 
R {9 ^ 71/2 — 9, i.e., sin6' -H- cos 9, and C, ^ —() and parity P (a -H- a^) turns out 
to be a symmetry operation. These three operations commute with each other. 
By the inclusion of PRX, the walker evolves by {[PRX]SBy. At each step, the 
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walker evolves according to 



Wp = PR 




e*^cos(6') e*'' sin(6') 



sin(6') — e *^ cos(6') 



cos{e) 



(3.14) 



This is equivalent to replacing B by 5*^^^ with = vr, which, according to Theorem 
1, should leave the walk distribution invariant. Thus the operation PRX applied at 
each step, is a symmetry of the quantum walk. It will be convenient henceforth to 
choose ^ = C = 0, so that R will simply correspond to the replacement 6 tx /2 — 6. 

Representing the inclusion of operations P, R or X at each step of the walk by P, 
R or X, respectively, we express this symmetry by the statements: 



W 
XW 



PRXW, 
PBW. 



(3.15a) 
(3.15b) 



The above expression (3.15a) was proved immediately above and (3.15b) follows 
from (3.15a), since the operations P, R and X mutually commute, and X"^ = I. 
It expresses the fact that applying the X operation at each step is equivalent to 
replacing a quantum walk by its angle-reflected, spatially inverted counterpart. 
The observation made at the beginning of this section pertains to the special case 
of ^ = 45°. By a similar technique the following symmetries may be proved. 



XW - XZW - ZXW. 



(3.16) 



The first equivalence easily follows from (3.13). 



3.2.2 Environmental effects 

A quantum walk implemented on a physical system is is inevitably affected by 
noise due to the environment. We consider three physically relevant models of 
noise: a phase flip channel (which is equivalent to a phase damping or purely de- 
phasing channel), a bit flip channel and a generalized amplitude damping channel 
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Figure 3.1: The effect of envhonmental decoherence on the position proba- 
bility distribution of a quantum walk subjected to a noisy channel. Coin bias 
is of the form (2.39) with 6 = 60°,^ = C = 0. The noise is modeled as a phase 
flip (solid line) and bit flip (dashed line) channel, characterized by (3.24) and 
(3.30), respectively, at various noise levels p : (a) p = 0.005 ; (h) p = 0.05; 
(c) p = 0.1; (d) p = 0.5, which corresponds to a fully classical random walk. 
Comparing Figure (d) with Figure 3.2(d), we note that the distribution in the 
case of maximal bit flip noise is the same. The distribution is for 100 steps. 



(T > 0). In all cases, our numerical implementation of these channels evolves 
the density matrix employing the Kraus operator representation for them. How- 
ever to explain symmetry effects, it is convenient to use the quantum trajectories 
approach, discussed below. 
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Figure 3.2: The effect of envhonmental decoherence on the position proba- 
biUty distribution of a quantum walk subjected to a noisy channel. Coin bias 
is of the form (2.39) with 6 = 30°,^ = ( = 0. The noise is modeled as a phase 
flip (solid line) and bit flip (dashed line) channel, characterized by (3.24) and 
(3.30), respectively, at various noise levels p: (a) p = 0.005 (b) p = 0.05; (c) 
p = 0.1; (d) p = 0.5, which corresponds to a fully classical random walk. The 

distribution is for 100 steps. 



3.2.2.1 Phase damping and bit flip channels 

In studying the status of the walk symmetries in the presence of noise, it is advan- 
tageous to employ the quantum trajectories approach [101]. This simplifies the 
description of an open quantum system in terms of a stochastically evolving pure 
state, which allows us to adapt the symmetry results for the pure states, given in 
the preceding section, to mixed states. 
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Figure 3.3: Variation of standard deviation with noise level, for both phase 
noise and bit flip noise, (a) and (b) is for 9 = 30°; (c) is for Hadamard walk 
(6 = 45°); (d) and (e) is for 6 = 60° in the quantum coin operation i?o,e,o- 
In the classical limit of p = 0.5, the standard deviation converges to a fixed 
value for bit flip noise, irrespective of 6, but different for phase flip noise. The 
convergence happens because, at maximum bit flip noise {p = 0.5), the mea- 
surement outcome in the computational basis is completely randomized. On 
the other hand, the non-convergence in the case of phase flip noise is due to 
the fact that the asymptotic mixed state obtained via a phase damping channel 
depends on the value of the 6. The standard deviation was calculated for 100 



a 'quantum trajectory'. (More precisely, a trajectory refers to the sequence of 
states produced by these operations, for which the above serves as a convenient 
representation.) If all the B^s with x taking values from 1 to t are the same, 
then X is the usual 'homogeneous' quantum walk W. In general, the B^s may be 
different f/(2) operators, B^^\ More generally, each step of the walk may include 
generalized measurements whose outcomes are known (Section 3.2.2.2). If each 
walk step in X is subjected to a fixed symmetry operation G, the result is a new 
quantum trajectory 



steps of quantum walk. 



We call the sequence of walk step operations. 




(3.17) 
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We have the following generalization of Theorem 1 to inhomogeneous quantum 
walks on a line. 

Theorem 2. Given any quantum walk trajectory X = (SBt) ■ ■ ■ (SBi), the sym- 
metry Q holds, i.e., X ^ GX. If the operation $ (3.11) is restricted to Z, that is 
operation Z at each step of the walk, then the symmetries hold even when some 
of the S"s are replaced by S'^'s. 

Proof. In the proof of Theorem 1, we note that if, in each step of the walk, we 
alter the rotation B by the transformation G, the proof still goes through. That 
is, l^i) = {SBt){SBt^i)---{SB,)\a,^) and j^a) = {SBt){SBt,) ■ ■ ■ {SBt)\a, P) 
produce the same position distribution. 

Suppose that in some of the walk steps, S is replaced with S'^. In place of (3.7), 
we have 

l^i) = {SBt)---{S^B,)---{SB,)\a,f3) 

K,'it-i---K2,qAi,'^\it^(^ + ^{Qi-Q2)-{ti-t2)), (3.19a) 

qi,q2,--- ,qt 

\^,) = {SBi'^)...{S^Bi'y)...{SB{'^)\a,(3) 

= J2 ■ ■■K,,,AiA-^r-'^"''""^''\Qt,^ + 2{Q, - Q2) - (ti - t2)), 

qimr-- ,qt 

(3.19b) 

where Qi = J2k Ik ^^e ti steps k where operator S is used, and Q2 = J2i 1i 
for the t2 steps / where operator S*^ is used. Here Q = Qi + Q2 and t = ti + ^2- 
Observe that the exponent of (—1) is effectively evaluated in modulo-2 arithmetic. 
We can thus replace Q by Qi — Q2 in the exponent. Following the argument in 
Theorem 1, we find that 

(a,6|^i) = e*®(a,6|^2), (3.20) 
where Q = Qi — Q2 — a + a. ■ 

As a corollary, the symmetries Z and PRX hold good because they reduce to 
special cases of G. A question of practical interest is whether PRX and Z are 
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symmetries of a noisy quantum walk. Suppose we are given a noise process TV in 
the Kraus representation : 

p Af{p) = ExpEl Yl ^1^- = I- (3-21) 

x=0 X 

With the inclusion of noise, each step of the quantum walk becomes augmented 
to {USBk), where 11 is a random variable that takes Kraus operator values E^. 
Thus, A/" corresponds to a mixture of upto trajectories or 'unravellings' 

Xi = {U{k)SBt)---{U{h)SB,), (3.22) 

each occurring with some probability pi, where YliPi ~ ^- and PRX are 

symmetries of an unraveling Xi, then the operations Xi and 

BXi = {U{lt)DUBt) ■ ■ ■ {U{h)DUBi) , (3.23) 

where D denotes Z or PRX, must yield the same position probability distribution. 
In the case of bit-flip and phase-flip channels, there is a representation in which 
the i?a;'s are proportional to unitary operators. 

Theorem 3. If trajectories Xi are individually symmetric under operation G, then 
so is any noisy quantum walk represented by a collection {Xi^pi}. 

Proof. The state of the system obtained via A/" is a linear combination (the 
average) of states obtained via the ^^..'s. Thus, the invariance of the X^^s under 
G implies the invariance of the former. ■ 

This result, together with those from the preceding Section, can now be easily 
shown to imply that the symmetry D is preserved in the case of phase-flip and 
bit-flip channels. 

Decoherence via a purely dephasing channel, without any loss of energy, can be 
modeled as a phase flip channel [5]: 

S{p) = il-p)p + pZpZ. (3.24) 
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An example of a physical process that realizes (3.24) is a two-level system inter- 
acting with its bath via a quantum non-demolition (QND) interaction given by 
the Hamiltonian 

H = + H^, + 

k k k 

Here H5 stands for the Hamiltonian of the system (S), for our noise channels the 
system is only the particle (coin degree of freedom). H/j and tlsR are reservoir 
(R) and system-reservoir (SR) interaction, respectively. Operators h and W are 
the annihilation and creation operators. The last term on the RHS of (3.25) is a 
renormalization inducing 'counter term'. Since [H5,H5/j] = 0, (3.25) is of QND 
type. 

Following [102] (apart from a change in notation which switches |0) i — > |1)), 
taking into account the effect of the environment modeled as a thermal bath, the 
reduced dynamics of the system can be obtained, which can be described using 
Bloch vectors as follows. Its action on an initial state 

1(1H-(M0)» (.40)) 
\ {a+(0)) \(l-(aM))l 

is given in the interaction picture by 

> I 1(1 +{'^.(0))) (<T_(0)>e-<M'7(.)\ 
= I (a,(0)>e-(M^«. l(l-(..(0)»j- 

Here (j+ and a_ are the standard raising and lowering operators given by 

a+ = |l)(0| = l(a, + ^a,) ; a_ = |0)(1| = ^(a, - ^a,) (3.28) 

where ax-, (Jy and az are the Pauli spin operator in x, y and z directions, respec- 
tively. The initial state (3.26) may be mixed. (The derivation of the superoperator 
S in terms of environmental parameters for the pure state case, given explicitly 
in [102], is directly generalized to the case of an arbitrary mixture of pure states. 
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since the environmental parameters are assumed to be independent of the system's 
state.) 

Comparing (3.27) with (3.24) allows us to relate the noise level p in terms of 
physical parameters. In particular, 

p=^(l-exp[-(M'7W])- (3.29) 

When 7(t) ^0 (either because the coupling with the environment is very weak or 
the interaction time is short or the temperature is low), p ^ 0, tending towards 
the noiseless case. On the other hand, under strong coupling, 7(t) is arbitrarily 
large, and p — )■ 1/2, the maximally noisy limit. The result of implementing channel 
(3.24) is to drive the position probability distribution towards a classical Gaussian 
pattern [100]. The effect of increasing phase noise in the presence on quantum 
walk is depicted in Figure 3.1, for the case of ^ = 60°, and in Figure 3.2, for the 
case of 6* = 30°. The onset of classicality is observed in the Gaussianization of the 
probability distribution. This is reflected also in the fall of standard deviation, as 
shown in Figure 3.3. 

Decoherence can also be introduced by another noise model, the bit flip channel 

[5], 

S{p) = {l-p)p + pXpX. (3.30) 

As with the phase damping channel, the bit flip channel also drives the probability 
distribution towards a classical, Gaussian pattern, with increasing noise [100]. The 
effect of increasing bit flip noise in the presence of biased walk is depicted in Figures 
3.1 and 3.2. Here again, the onset of classicality is observed in the Gaussianization 
of the probability distribution, as well as in the fall of standard deviation, as shown 
in Figure 3.3. 

A difference in the classical limit of these two noise processes, as observed in Figure 
3.3, is that whereas the standard deviation (in fact, the distribution) is unique in 
the case of the bit flip channel irrespective of 6, in the case of phase flip noise, the 
classical limit distribution is 9 dependent. This is because phase flip noise leads, 
in the Bloch sphere picture, to a coplanar evolution of states towards the o"^ axis. 
Thus all initial pure states corresponding to a fixed 6 evolve asymptotically to the 



66 



Chapter 3. Symmetries and noise 



5 



0, 





---15° 

- - 30° 

— 45° 
° 60° 
o 75° 




o o. o . 



0.1 0.2 0.3 0.4 0.5 

Bit flip noise level 



Figure 3.4: Variation of the ratio of standard deviation without any symmetry 
operation to the bit flip symmetry operation with increasing bit flip noise level. 
The standard deviation was calculated for 100 steps of walk. 
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Figure 3.5: Variation of the ratio of standard deviation without any symmetry 
operation to the phase flip symmetry operation with increasing bit flip (phase 
flip) noise level. The standard deviation was calculated for 100 steps of walk. 



same mixed state [5, 102]. This also explains the contrasting behavior of bit flip 
and phase flip noise with respect to 9, as seen by comparing Figures 3.1 and 3.2. 



Representing the walk distribution by its standard deviation a, we may describe 
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Figure 3.6: Variation of the ratio of standard deviation without any symmetry 
operation to the bit flip symmetry operation with increasing phase flip noise 
level. The standard deviation was calculated for 100 steps of walk. 



symmetry by the ratio of a without the symmetry operation to a with the sym- 
metry operation. Figure 3.4 depicts the symmetry operation X for various bit flip 
noise levels. The convergence of the curves representing various 6's is a consequence 
of the complete randomization of the measured bit outcome in the computational 
basis. This implies that although X is not a symmetry of non-Hadamard walk 
(6* 7^ 7r/4), it does become one in the fully classical limit. On the other hand, the 
symmetries PRX remain unaffected by noise. We note that, since the quantum 
walk here is evolved from the symmetric state ^ (|0) + and the bit flip and 
phase flip noise are not partial to the state |0) or this is equivalent to setting 
P to 1, which explains the fact that the distributions in Figures 3.1 and 3.2 are 
spatially symmetric. Thus, RX by itself becomes a symmetry operation, which is 
manifested in the fact that in Figure 3.4 the values of the curve for complementary 
angles are the inverse of each other. Figure 3.5 shows that for either of the two 
noises, Z is a walk symmetry. 

Figure 3.6 depicts the symmetry of the RX operation at all phase flip noise levels, 
as evident from that fact that the values of the curve for complementary angles 
are the inverse of each other. From Figures 3.4, 3.5 and 3.6, we note that for 
the Hadamard walk, all three symmetry operations Z,X and R are individually 
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preserved. This is expected because here P = 1 as stated earher, R = 1 by 
definition, so that the symmetry of PRX imphes X = 1. 

With the non-Hadamard coin operation and an initial arbitrary state, the full 
symmetries Z and PRX would be required, as proved by the following theorem. 

Theorem 4. The operations PRX and Z are symmetries for the phase-fiip and 
bit-fiip channels. 

Proof. We may look upon the phase flip channel (3.24) as a probabilistic mixture 
(in the discretized walk model) of 2* quantum trajectories with 11 G {I, Z}. By 
virtue of Theorem 3, it suffices to show that any given unraveling is invariant 
under Z and PRX. Consider an unraveling 

= ■ ■ ■ {I SB) (ZSB) {I SB) ■■■ = ■■■ (SB) (ZSB) (SB) ■■■ . (3.31) 

This is the same as: ■ ■ ■ {SB){SB'){SB) ■ ■ ■ , where B' = ZB, noting that Z 
commutes with S. Now, 

ZX^ = ■■■ {IZSB){ZZSB){IZSB) ...= ..■ {ZSB){ZSB'){ZSB) ■ ■ ■ , (3.32) 

which, by Theorem 2, is equivalent to X^. Also, by Theorem 2 

PRXX^ = ■ ■ ■ {PRX SB) {Z PRX SB) {PRX SB)--- 

= ---{PRX SB){PR ZX SB){PRX SB)--- 

= ---{PRX SB){PR ZXZ ZSB){PRX SB)--- 

= ---{PRX SB){PR{-X) SB'){PRX SB)-- - , (3.33) 

which, by Theorem 2, is equivalent to X^, since an overall phase factor of ±1 is 
irrelevant. Thus, the phase-fiip channel is symmetric with respect to the operations 
Z and PRX. 

Regarding the bit-fiip channel (3.30): as in the above case, consider an unraveling 

X^ = --- {SB) {XSB) {SB) -- - . (3.34) 
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This is the same as: ■ ■ ■ {SB){S^B"){SB) ■ ■ ■ , where, as may be seen by direct 
calculation, B" = XB. Now, 

ZX^ = ■■■{Z SB){X Z SB){Z SB)--- 

= ---{Z SB){X ZX XSB){Z SB)--- 

= ---{Z SB){{-Z) S^B"){Z SB)-- - , (3.35) 

which, by Theorem 2, is equivalent to X^, since an overall phase factor of ±1 is 
irrelevant. Further, 

PRXX^ = ---{PRX SB){X PRX SB){PRX SB)--- 
= ---{PRX SB){PRX X SB){PRX SB)--- 
= --- {PRX SB){PRX S^B"){PRX SB)-- - , (3.36) 

by Theorem 2, is also equivalent to X'^. ■ 



3.2.2.2 Generalized amplitude damping channel 

Here we study the behavior of quantum walk subjected to a generalized ampli- 
tude damping (with temperature T > 0), which would reduce at T = to the 
amplitude damping channel. As an example of a physical process that realizes the 
generalized amplitude damping channel, we consider a two-level system interact- 
ing with a reservoir of harmonic oscillators, with the system-reservoir interaction 
being dissipative and of the weak Born-Markov type [103, 104] leading to a stan- 
dard Lindblad equation, which in the interaction picture has the following form 
[105] 

d ^ 

= Yl (S^'^P'^I - ^I^fcP' - P'T^iT^k) , (3.37) 

k=l 

where 

7^l = (7o(A^th + l)/2)l/'7^ ; 7^2 = (7oiVth/2)^/'7^t (3.38) 

and 

iVth = {exp{hw/kBT) - l)-\ (3.39) 
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Figure 3.7: Amplitude damping channel acting on a Hadamard walk at tem- 
perature T = 0. The distribution corresponding to intermediate values of p 
clearly show the breakdown of the RX symmetry. However, the extended sym- 
metry, PRX (where P stands for parity operation (spatial inversion)) holds 
good. This is observed at all temperatures, (a) Probability distribution of find- 
ing the particle on which amplitude damping channel is acting. This shows that 
even at T = 0, for sufficiently high coupling, the distribution turns classical, (b) 
Amplitude damping channel with a bit flip symmetry. 



is the Planck distribution giving the number of thermal photons at the frequency 
w, and 7o is the system-environment coupling constant. Here 



7^ = (T_ cosh(r) + e'V+ sinh(r), (3.40) 



and the quantities r and 5 are the environmental squeezing parameters and 

cr± = ^ (cr^ ± iuy) . (3.41) 

For the generalized amplitude damping channel, we set r = 5 = 0. If T = 0, so 
that A'th = 0, then 7^2 vanishes, and a single Lindblad operator suffices. The 
generalized amplitude damping channel is characterized by the following Kraus 
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Figure 3.8: PRX symmetry seen to hold in quantum walk (with 9 = 60° and 
9 = 30°) subjected to amplitude damping (T = 0). The two cases are spatial 
inversions of each other. This holds for a generalized amplitude damping at any 
temperature, (a) walk with 9 = 60° ; (b) 9 = 30° and bit flip. 




Figure 3.9: Onset of classicality is seen to be accentuated in a Hadamard walk 
subjected to generalized amplitude damping followed by bit flip with increas- 
ing temperatures (without bit flip the damping would have been to the right). 
Figure 3.7(a) depicts the T = case (x = 1 in (3.43)). (a) Finite tempera- 
ture corresponding to x = 0.75 (b) T = oo, corresponding to x = 0.5. It may 
be noted that even at T = oo, for sufficiently small coupling the distribution 

remains non-classical. 
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Figure 3.10: Variation of standard deviation with amphtude damping noise 
level for various value of 9, 15°, 30°, 45°, 60° and 75°. Note that the standard 
deviation for complementary angles converge to the same value. 



operators, 



1 

v^T^pW 



E2 = VT^ 



l-p{t) 
1 



E3 = VT^ 











where 



= 1 _ e-7o(2^«,+i)t. 
When T = 0, X = 1, and for T ^ cx), x = 1/2. 



X - 2 



1 + 



2iVth + 1 



The density operator at a future time can be obtained as [105] 



A2 

;i-^i) 



where 



(3.42) 



(3.43) 



(3.44) 



A, = {a-M) = e-^°(^^-+^)*(a.(0)) - (l - e-^"^^^-^^)*) , (3.45) 



^2 = e-^(27V,, + l)t^^^^Q)^^_ 



(3.46) 
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Figures 3.7, 3.8 and 3.9 depict the onset of classicahty with increasing couphng 
strength (related to p) and temperature (coming from x)- Figure 3.7, which shows 
the effect of an amphtude damping channel on a Hadamard walk at zero temper- 
ature, illustrates the breakdown of RX symmetry even though the initial state is 
^(|0) + This is because, in contrast to the phase-flip and bit-flip channels, 
the generalized amplitude damping is not symmetric towards the states |0) and 
However, the extended symmetry PRX is preserved both for Hadamard as 
well as quantum walks with 6 ^ n /2, as seen from Figures 3.7 and 3.8, respectively. 

From Figures 3.7(a) and 3.9(a,b), the onset of classicality with increasing tem- 
perature is clearly seen. Figure 3.10 presents the standard deviation for quantum 
walks on a line with various biases, subjected to amplitude damping noise. The 
standard deviation for complementary angles {6 ^ 71/2 — 9) is seen to converge 
to the same value in the fully classical limit. This may be understood as follows. 
First, we note that since PRX is a symmetry of the quantum walk, and the ef- 
fect of P does not show up in the standard deviation plots, RX by itself is an 
apparent symmetry. Further, in the classical limit the measurement outcome be- 
ing a unique asymptotic state for the (generalized) amplitude damping channel, 
effectively X ~ 1, which makes R a symmetry operation. 

The following theorem generalizes Theorem 2 to an open system subjected to a 
generalized amplitude damping channel. 

Theorem 5. The operations Z and PRX are symmetries for the generalized am- 
plitude damping channel. 

Proof. By virtue of Theorem 3, it suffices to show that any given unraveling is 
invariant under Z and PRX. Consider an unraveling 



■ ■ ■ {EoSB){E,SB){E2SB){E:,SB) ■ ■ ■ 
■■■(5i?(o))(5ti?(i))(5i?(2))(5ti?(3))---, 



(3.47a) 
(3.47b) 
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where the non-unitary matrices are given by B^^x) = E^B. Now, 

ZX^ = ■■■{EoZSB){EiZSB){E2ZSB){E3ZSB)--- 

= ■■■ {ZEoSB){{-Z)EiSB){ZE2SB){{-Z)E3SB) ■ ■ ■ 

= ■■■ {ZSB(,)){{-Z)S^B(r,){ZSB^,^){{-Z)S^B^,)) ■ ■ ■ 

= ■■■ {SB[l]){-S^B[l])iSBll])i-S^Bll]) .... (3.48) 

Ignoring the overall ±1 factor in (3.48), and comparing it with (3.47a), and noting 
that that the derivation of the proof of Theorem 2 did not require the matrices Bj 
to be unitary, we find along similar lines that ZX^ is equivalent to X^. 

The following may be directly verified 

PRXX^ = ... (EoPRXSB) (EiPRXSB) {E2PRXSB) {E3PRXSB) . . .(3.49a) 
= ...{EoSB^^^){EiSB'-^^){E2SB^^^){E3SB^^^)... (3.49b) 
= • • . {SB^;^){S^b'^]){SB^^){S^B^;^) . . . , (3.49c) 

which, by Theorem 2, is equivalent to X^. For proof of (3.49b), see the proof of 
Theorem 1. Equation (3.49c) is obtained analogously to (3.47b), except that the 
matrix B^'^^ is used instead oi B. ■ 

This may be expressed by the statement 

AfW - MZW, (3.50a) 
MW - A/'PRXl?, (3.50b) 

which generalizes (3.10). These results show that the symmetries persist for de- 
phasing (phase flip), bit flip and (generalized) amplitude damping channels. 

3.3 Symmetry and noise operations on an n— cycle 

In previous section we considered symmetries for a quantum walk on a line, and 
the influence of noise on them. In this section we will extend the ideas to quantum 
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walks on an n— cycle. We will show the breakdown of symmetry when the quantum 
walk is implemented on an n— cycle, the effect of noise and its influence on the 
restoration of symmetry. 



3.3.1 Breakdown in symmetry 

In contrast to the case of quantum walk on a line, none of the four discrete sym- 
metries of Theorem 1 hold in general for unitary quantum walk on an n— cycle or 
closed path. Thus, if B (2.33) is replaced by any of 5^, B'-^\ B^^\ or 5^, Efiven 
by (3.5), the spatial probability distribution is not guaranteed to be the same. 

Theorem 6. The operation G : B ^ B* is in general not a symmetry of the 
quantum walk on an n— cycle. 

Proof. For the cyclic case, in place of (3.7), we now have 

= {SBY\a,l3)= Kt,<it-i---K2,qAi,o.\qt,/3 + 2Q-t {modn)), 

(3.51a) 

\^>,) = (^i?W)*|a,/3) 

= E ^.*,.*-i---^.2,.i&.i,a(e^'^r~^+-+'^^+1gt,/3 + 2g-t (modn)), 

91,92, •■• ,9t 

(3.51b) 

where n is the number of sites in an n— cycle. For an arbitrary state |a, h) in the 
computational-and-position basis, we have 

(a,6|^l) = ^ &a,gt_i • ■ ■&g2,gi^gi,a (3.52a) 

(a, h\^,) = ■ ■ ■ bq2,qAUe'^r-'^---^''^'' 

Y ba,q.^^---K,<iAUe'^n, (3.52b) 
iji,92,--- ,9t-ieS 

where Q is the set of binary (t — l)-tuples qi,q2, - ■ ■ At-i such that 



Q = gi + 92 H h qt^i + a 

76 
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satisfies 



h = l3 + 2Q-t (mod n). 



(3.54) 



We find that 



e = a - a + Q and Q = {h + t - fi) /2 + mR, 



(3.55) 



where m = 0,1,2, ■■■ , \ t/n\. Thus, the terms in the superposition (3.52a) are 
not in general identical with those in (3.52b), apart from a common factor, unless 
= 0, 27r, 47r, ■ ■ ■ . A similar argument can be used to show that the terms in 
(a, are not in general the same as those in (a, b\'^2)- Given the independence 

of from the coefficients fegz.gi^gi.a' is not necessary that 



The equality holds in general (for arbitrary unitary matrix B and time t) if and 
only if = 0, 2it, Air, ■ ■ ■ . Repeating the argument for B^'^\ B^^^ and B^'^\ we find 
that all the four discrete symmetries of Theorem 1 break down in general. ■ 

We note that for phase flip symmetry, where e^'^ = —1, the superposition terms in 
(3.52b) may differ from the corresponding terms in (3.52a) only with respect to 
sign. Given that all the terms like bg^^q^, etc. are built from a small set of 

trigonometric functions of the three parameters 6, ( and ^, certain values of t may 
render the right hand sides of (3.52a) and (3.52b) equal. However in general, this 
equality will not hold for arbitrary t. 

An instance of breakdown of phase flip symmetry in the unitary quantum walk on 
a cycle is demonstrated in the example of Figure 3.11. The profile of the position 
probability distribution varies depending on the number of sites and the evolution 
time. 



\{a,b\^^)\' + \{a,b\^i)\' = \{a,b\^,)\' + \{a,b\'^2)\' ■ 



(3.56) 
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Figure 3.11: An instance of breakdown of phase flip symmetry in a unitary 
quantum walk on a cycle, where the two extreme points (located at positions 
±50) on the plot are spatially adjacent. The number of sites is 101 and t = 5000 
(in units of discrete time-steps), with coin parameter 6 = 30°. (a) The solid 
curve represents the positional probability distribution without any symmetry 
operation applied, while the dashed curve represents that with a phase flip 
operation applied at each walk step, (b) The same as the above, but with time- 
averaging applied over every 50 steps, in order to more clearly bring out the 

breakdown in symmetry. 



3.3.2 Breakdown using a generalized phase gate 

We will introduce a generalized phase gate, an element from the parameter group 



G{/3) = 




(3.57) 



to act on the "He- We find that the operation W — > GW leaves the probability 
distribution p{j, t) of the particle on the line invariant; hence the walk is symmetric 
under the operation 

G{(3) -.Iq) ^e'^^lq) (3.58) 

for |g) in the computational basis (eigenstates of the Pauli operator cr^) and 
g = 0, 1. The physical significance of G is that it helps identify a family of quantum 
walks that are equivalent from the viewpoint of physical implementation, which 
can sometimes allow a significant practical simplification [64]. For example, sup- 
pose the application of the conditional shift is accompanied by a phase gate. The 
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walk symmetry implies that this gate need not be corrected for, thereby resulting 
in a saving of experimental resources. The inclusion of a phase gate on the coin 
operator is equivalent to a phase gate at each lattice site in the sense of quantum 
lattice gas automata (QLGA) [24, 25], with the physical meaning of a constant po- 
tential. The evolution rules for single-particle QLGA can be classified into gauge 
equivalent classes, there being a difference between the class of rules for periodic 
(n— cycle) and non-periodic 1-dimensional lattice and this feature can be used to 
distinguish between these two spatial topologies [106]. 

It turns out that in the case of quantum walk on an n-cycle, this symmetry breaks 
down. To see this, we note that the t-fold application of the operation GSB 
on a particle with initial state \^in) on the line and on an n— cycle produces, 
respectively, the states 

qi,q2,--- ,qt 

(3.59a) 

(GSBYl^,.^) = 

e''^'^Bg,^,^_,---B^,^g,{B^,oa + Bg,,b)\qt)\2Qt-t mod n), 

gi,(j2,--- ,qt 

(3.59b) 

where Qt = q'i + q'2 + ' ■ ■ + All terms in superposition (3.59a) contributing to the 
probability to detect the walker at a given position j = 2Qt—t have the same phase 
factor, e**^*^, which is fixed by Qt = {j + t)/2 (where, it may be noted, j and t are 
both even or both odd). Thus, this factor does not affect the probability to detect 
the walker at j, whence the symmetry. In the case of quantum walk on an n-cycle 
the breakdown of the symmetry, see Figure 3.12, can be attributed to the topology 
of the cycle, which introduces a periodicity in the walker position (determined by 
a congruence relation with modulus given by the number of sites), but not in the 
phase of the superposition terms. As a result, fixing j fixes Qt mod n = {j +t)/2 
mod n, but not Qt itself, so that the phase terms in the superposition (3.59b) do 
not factor out globally. Thus if /3 is non- vanishing, then in general the symmetry 
G is absent in the cyclic case. We quantify the breakdown in symmetry by means 
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Figure 3.12: Position probability distribution for a Hadamard walk, 
i?(0°, 45°, 0°) on (a) a line, and (b) an n-cycle, with (n = 51) and initial state 
(l/-v/2)(|0) + i|l)), for the unitary case (70 = 0) and r = 50. Each figure 
presents the distribution with and without being subjected to the phase opera- 
tion G(40°). In (a), there is perfect symmetry, since both distributions coincide. 
In (b) the two plots do not overlap, indicating the breakdown of the symmetry. 




Particle position 



Figure 3.13: Restoration of phase flip symmetry in a noisy quantum walk 
on a cycle, where the two extreme points on the plot are spatially adjacent. 
The number of sites is 101 and t = 5000 (in units of discrete time-steps), with 
coin parameter 6 = 30°. The figure depicts the position probability distribution 
with or without phase flip symmetry operation applied at each step, with phase 
damping noise level p = 0.02 (3.24). After sufficiently long time, the quantum 
walk reaches the uniform distribution, typical of classical random walk. 
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of the Kolmogorov distance (or, trace distance [5]), given by 



d(t) = (l/2)^|P(j,t)-P'(j,t)|, 



(3.60) 



3 



between the particle position distributions obtained without and with the symme- 
try operation, given by P{j,t) and P'{j,t), respectively. The breakdown in sym- 
metry for a noiseless cyclic quantum walk is depicted by the bold curve in Figure 
3.14 as a function of the number of turns r (where t = ts, with n = 2s + 1). 

3.3.3 Effect of noise and symmetry restoration 

Remarkably, this symmetry is restored above a threshold value of noise. The 
pattern in Figure 3.13 corresponds to phase noise with p = 0.02 applied to a 
quantum walk, either with or without a phase flip symmetry operation. We note 
that the introduction of noise tends to classicalize the random walk, hence causing 
it to asymptotically reach a uniform distribution [76, 107]. The above-mentioned 
symmetry restoration happens well before the uniformity sets in. The initial lack 
of symmetry gradually transit to full symmetry as the noise level is increased. 
Thus, the role of symmetry operations and noise is quite different in the case of 
quantum walk on an n— cycle as compared with that on a line. 

We now describe the n-cycle quantum walk on the particle, a two level system, 
when subjected to noise. The situation is modeled as an interaction with a thermal 
bath, characterized by phase damping or a generalized amplitude damping channel, 
the latter process being represented by the Kraus operators (3.42). The density 
operator of the coin evolves according to 



The full evolution of the walker, described by density operator p(t), is given by 




(3.61) 



X 




(3.62) 



3 



where the E^s are understood to act only in the coin space. 
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Figure 3.14: Kolmogorov distance (i(r) against the number of turns (r) of the 
cychc quantum walk in the noiseless and noisy case with n = 51. For the unitary 
case (A = 0, 70 ; bold line) the walk becomes increasingly asymmetric as the 
number of turns is increased, until about 7-10 turns, after which it fluctuates 
around d ~ 0.15. The plots represent generalized amplitude damping noise at 
different noise levels at temperature T = 3.5 (in units where h= ks = 1). The 
walk is evolved with the initial state of the particle parameters 5 = 30° , rj = 40° 
in (2.11), with 5(20°, 10°, 30°) and G(10°). 



The curves in Figure 3.14 plot d{T) as a function of turns in the case of unitary and 
noisy quantum walk (parametrized by 70), and demonstrate the gradual restora- 
tion of symmetry with time on account of the noise. Although the figure employs 
generalized amplitude damping noise, qualitatively the same behavior can be seen 
for a phase damping noise. Here, a general feature is that (i(r) is non-zero when 
r < 2, being equivalent to (noisy) quantum walk on a line. Thereafter, (i(r) at first 
increases with increasing turns, being dominated by unitary evolution, and even- 
tually falls down, being dominated by noise. It is observed that for sufficiently low 
noise levels, the time at which this turnover in slope happens remains constant, 
for given 6. This is depicted in Figure 3.14 for the case of a generalized amplitude 
damping channel corresponding to a fixed temperature and varying 70. However, 
we note that for strong enough noise, the turnover happens earlier. 

Typical noisy probability distributions are depicted in Figure 3. 15(a), (b) at an 
instant where the symmetry has been almost fully restored while the walk is well 
within the quantum regime. Figure 3.15(c) represents a classicalized distribu- 
tion, indicated by the regular envelope (that will eventually turn into a uniform 
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Figure 3.15: Position probability distribution for a Hadamard walk 
S(0,45°,0) on an 77,-cycle (n = 51) with initial state (l/\/2)(|0) + i|l)) when 
subjected to generalized amplitude damping noise with A = 0.1 and finite T 
(= 6.0) for T = 11. For (a) and (b) 70 = 0.025. (b) is the distribution for 
the quantum walk augmented by operation G(20°) and (a) is without any aug- 
mented operation; note the overlap of (b) on (a), walk remains quantum yet 
with symmetry completely restored ; (c) classicalized pattern (indicated by the 
regular envelope) obtained with larger noise level corresponding to 70 = 0.1. 

distribution). 

We define coherence C as the sum of the off-diagonal terms of states in Tic ® Tip, 
where He and Tip are the coin and position Hilbert space of the quantum walker, 
respectively. If the state of the quantum walker is 



the sum of the absolute values of all off-diagonal terms of p in the computational- 
position basis. The coherence function is defined as the quantity C{m), where 
m G {1, 2, ■ ■ ■ , M}, obtained by partitioning C into M intervals of size s/M, such 




(3.63) 



ab;jk 



where |a), \b) G Tic and \k) G Tip, then 




(3.64) 



83 



Chapter 3. Symmetries and noise 



that for the mth interval 



(m - l)(s/M) < \j -k\ < m{s/M). 



(3.65) 



Physically, C{m) is a measure of coherence between two points on a (in general, 
noisy) quantum walker, as a function of their mutual separation. Let Co(m) 
represent the coherence function of the corresponding noiseless walk. At any turn 
r, we define the normalized coherence function by 



Since noise tends to destroy superpositions, and the breakdown in symmetry is 
essentially a phenomenon of superposition of the forward and backward waves, 
noise tends to restore symmetry, as seen from Figure 3.15. This is brought out by 
Figure 3.16 for two possible values of G. In the figure, in spite of its considerable 
spikiness, the bold curve, representing c(m = M), shows an overall fall. A similar 
trend as depicted in this figure, has been numerically checked for various other 
values of G. This raises the question whether symmetry restoration of the cyclic 
quantum walk can be considered as a good indicator of classicalization. Here we 
note that from Figure 3.15(a),(b) the probability distribution pattern is seen to 
be clearly quantum, even though symmetry has been almost fully restored. This 
is suggestive of the notion that that symmetry tends to be restored even in the 
regime where the walk still possesses some quantum features. 

The reason D{t) is not a faithful indicator of classicalization of the walk has to 
do with the effect of noise on the sensitivity of the symmetry operation G{(5) 
to the topology of the path. Since this operation senses the closure of the path 
through the superposition of the forward and backward waves, the suppression of 
superposition through noise will also have the effect of desensitizing the operation 
to the closure of the path, thereby moving the noisy cyclic quantum walk towards 
a noisy quantum walk on a line from the perspective of this operation, before 




(3.66) 



and, analogously, normalized Kolmogorov distance by 




(3.67) 
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Figure 3.16: Normalized coherence function c(M) (bold) and the normalized 
Kolmogorov distance D(t) (dashed) for an n— cyclic quantum walk as a function 
of turns T. An overall reduction of both c{M) and D{t) with time is seen for the 
generalized amplitude damping noise characterized by temperature T = 3.5 and 
A = 0.1. The Hadamard walk is evolved with the initial state parameters (in 
degrees) 6 = 45°, rj = 0° (2.11), and G{f3). The line-with-circle plot represents 
D{t) with G(20°) and dashed line represents for G(40°). The c(M) (sohd line) 
remains roughly the same for both /3. For a clear depiction of the notion of 
symmetry restoration in the walk even in the quantum regime. Figure (3.15). 



further classicalization transforms it into a cychc classical random walk. And 
as shown in Section 3.2, all the above symmetries are respected by a (noisy) 
quantum walk on a line, both in the case of phase damping noise and generalized 
amplitude damping^. This brings out the point that decoherence (T2 process) is 
the principal mechanism responsible for the restoration of symmetries. It also 
highlights the interplay between topology and noise in a quantum walk on an 
n— cycle. A similar interplay may be expected also in the case of quantum walk 
with other nontrivial topologies. 

Further extensions would be quantum walks on a more general graph [76, 107] or 
in higher dimensions d > 2 [109]. In the former, the ID walk is generalized to an 
n-cycles and to hypercubes, including the effect of phase noise in the coin space, 
and decoherence in position space. In the latter, the Hadamard transformation is 

■^In NMR nomenclature, phase damping is called a T2 process, and generalized amplitude 
damping, which is a Ti ,T2 process [108]. 
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generalized to a non-entangling tensor product of Hadamards, or to an entangling 
discrete Fourier transform or the Grover operator. They bring in many novel 
features absent in the quantum walk on a line. 

The interplay between geometry and decoherence has been noted before in the 
case of delocalized bath modes [110], as against localized bath modes [110-112]. 
This is of relevance as the noise processes considered here [64, 105] are described 
by the interaction of the system with delocalized bath modes. 

3.4 Experimental implications 

Experimental realization of quantum walk using any of the proposed schemes is 
not free from noise due to environmental conditions and instrumental interference. 
In particular, noise can be a major issue in the scaling up of the number of steps 
in already realized quantum walk systems. Understanding the symmetries of the 
noisy and noiseless quantum walk could greatly help in the improvement of im- 
plementation technique (see for example Section 4.5.1) and in further exploration 
of other possible systems where quantum walk can be realized on a large scale. 

The experimental study of the decoherence and decay of quantum states of a 
trapped atomic ion's harmonic motion subjected to engineered reservoirs, both of 
the phase damping and amplitude damping kind, have been reported in [98, 99]. 
The phase reservoir is simulated by random variation of the trap frequency u 
without changing its energy (non-dissipative), while the amplitude reservoir is 
simulated by random electric field along the axis of the trap (dissipative). Coupling 
the reservoirs reported in [98, 99] to the scheme presenting the combination of 
pulses required to implement a quantum walk on a line and on a cycle in an ion 
trap [54] provides a convenient set up to demonstrate the symmetry-noise interplay. 
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3.4.1 NMR quantum- information processor 

Continuous time [50] and discrete time [51] quantum walk have been successfully 
implemented in a nuclear magnetic resonance (NMR) quantum-information pro- 
cessor. In NMR spectroscopy of the given system (molecule) the extent of its 
isolation from the environment is determined in terms of its phase coherence time 
T2 and its energy relaxation time Ti. If the pulse sequence is applied to the NMR 
quantum- information processor within the time T < T2,Ti, the environmental ef- 
fects on the system are less significant. The pulse sequence exceeding the time 
T2 can be considered to be affected by the dephasing channel and the pulse se- 
quence exceeding the time Ti can be considered to be affected by the amplitude 
damping channel. In experiments of time scale greater than the time T2 or Ti, a 
refocusing pulse sequence is applied to compensate for the environmental effects. 
In the implementation of discrete-time quantum walk, the pulse sequence was im- 
plemented within the time Ti and T2 and by introducing dephasing, the transition 
from quantum walk to the classical random walk was shown [51]. 

The environmental effect (noise) on quantum walk symmetries presented in this 
chapter can be verified in the NMR system by scaling up the number of steps of 
quantum walk realized. By applying a controlled amount of the refocusing pulse 
sequence, the effect of different levels of noise can be experimentally verified. 

3.4.2 Ultracold atoms 

There have been other proposals for physical realization of quantum walk using 
Bose-Einstein condensate (BEC) [57] where the unitary shift operator induces a 
bit flip. A stimulated Raman kick is used as a unitary shift operator to translate 
the Bose-Einstein condensate in the Schrodinger cat state to a superposition in 
position space. Two selected levels of the atoms in the Bose-Einstein condensate 
are coupled to the two modes of counter-propagating laser beams. The stimulated 
Raman kick, in imparting a translation in position space, also flips the internal 
state of the Bose-Einstein condensate. An rf pulse (vr pulse) is suggested as a 
compensatory mechanism to flip the internal state of the condensate back to its 
initial value after every unitary shift operator. From the PRX symmetry pointed 
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in this chapter, it follows that there is no need for the compensatory operation for 
an unbiased quantum walk for a particle initially in the state :^(|0) + The 
availability of walk symmetries could also be useful for exploring other possible 
physical implementations which induce such symmetry operations along with the 
translation. 

In the most widely studied version of quantum walk, a quantum coin operation 
is used after every displacement operation. Continuous external operations on a 
particle confined in a trap reduces the confinement time of the particle. Reducing 
the number of external operations will benefit the scaling up of the number of steps 
of the quantum walk. One can configure a system where a separate quantum coin 
operation is eliminated by using a single coin-embedded conditional shift operator 
to shift the particle in superposition of position space retaining the superposition 
of the coin state [63]. Further, systems of this kind are expected to be affected by 
amplitude damping as one of the basis states might be more stable than the other 
one in the trap. The present study of effect of noise could help to optimize and 
use quantum walk with in the system limitations. 

3.4.3 Other condensed matter systems 

Breaking of the symmetry due to the change in walk topology causes long-range 
correlations to develop, in analogous to the hydrodynamics of ordered systems such 
as spin waves in: ferromagnets, antiferromagnets (where it is the spin wave of stag- 
gered magnetization), second sound in He^, nematic liquid crystals [113]. Here the 
correlations may be identified with symmetry-broken terms (whose measurement 
probability depends on /3 in the walk augmented by G{f3)) in the superposition 
of the quantum walker. One finds that correlations are set up rapidly over large 
distances with increase in the winding of the walker, until symmetry is broken 
throughout the cycle. However, as noted above, the randomization produced by 
noise causes the reappearance of symmetries. The symmetry breaking and the 
symmetry restoring agents are thus different, the former given by the topological 
transition from a line to an n-cycle, the latter being the noise-induced randomiza- 
tion. 
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Coherence is also widely used to understand quantum phase transitions, the tran- 
sition from superfluid to Mott insulator state in an optical lattice being one specific 
example [81, 82]. In Chapter 4 the quantum phase transition using quantum walk 
in a one dimensional optical lattice is discussed [49]. Using various lattice tech- 
niques, desired geometries to trap and manipulate atoms can be created. In most 
physical situations one deals with closed geometries. The characteristics of the 
ra— cycle walk, in particular the re-appearance of the symmetry (implying a family 
of implementationally equivalent noisy cyclic quantum walks) while still in the 
quantum regime, presented here could be of direct relevance to such situations. 

The ubiquity of the ideas developed in this chapter can be seen from the fact that 
the quantum dynamics of a particle on a ring (cycle) subject to decoherence along 
with dissipation finds its place in the physics of quantum dots. The effective action 
of a quantum dot accounting for the joint effect of charging and coupling to an 
environment [114, 115] mirrors the behavior of the quantum dynamics of a particle 
on a ring (cycle) subject to a dissipative damping mechanism describing the dis- 
sipation of the energy stored in dynamic voltage fluctuations into the microscopic 
degrees of freedom of the quasi-particle continuum. In the absence of dissipation, 
the action describes the ballistic motion of a quantum particle on a ring. The 
ring topology reflects the 27r-periodicity of the quantum phase, which is in turn 
related to the quantization of charge, thereby highlighting the point that the main 
source of charge quantization phenomena, in the approach developed in [114, 115], 
is the periodicity, of the relevant variable, due to the ring topology. With the 
increase in the effect of dissipation, the particle begins to forget its ring topology 
(full traversal of the ring become increasingly unlikely), leading to a suppression 
of charge quantization phenomena. This behavior is similar to that seen here for 
the case of quantum walk on a cycle, where with an increase in the effect of the 
environment, i.e., with increasing noise, the walker in unable to perceive the cyclic 
structure of the walk space. That the topology-noise interplay studied here has an 
impact on a concrete condensed matter system, viz. the crossover from strong to 
weak charge quantization in a dissipative quantum dot, highlights the generality 
and scope of these ideas. 
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3.5 Summary 



• Our work considers variants of quantum walks on a line which are equivalent 
in the sense that the final positional probability distribution remains the 
same in each variant. In particular, we consider variants obtained by the 
experimentally relevant operations of Z or X applied at each quantum walk 
step, with the symmetry operations given by Z and PRX. This could be 
experimentally advantageous since practical constraints may mean that one 
of the variants is preferred over the rest. What is especially interesting 
is that these symmetries are preserved even in the presence of noise, in 
particular, those characterized by the phase fiip, bit fiip and generalized 
amplitude damping channels. This is important because it means that the 
equivalence of these variants is not affected by the presence of noise, which 
would be inevitable in actual experiments. The symmetry of the phase 
operation under phase noise is intuitive, considering that this noise has a 
Kraus representation consisting of operations that are symmetries of the 
noiseless quantum walk. However, for the PRX symmetry under phase noise, 
and for any symmetry under other noisy channels (especially in the case of 
generalized amplitude damping channel), the connection was not obvious 
before the analysis was completed. 

• An interesting fact that comes out while extending the studies to walk on an 
n— cycle is that the symmetry breaks down in general but is restored above a 
certain noise level. This symmetry-topology-noise interplay presented would 
be of relevance to quantum information processing systems, and have wider 
implications to the condensed matter systems. 

• Our results in this chapter to study noise model are supported by several 
numerical examples obtained by evolving the density operator in the Kraus 
representation. However, analytical proofs of the effect of noise on symme- 
tries are obtained using the quantum trajectories approach, which we find 
convenient for this situation. 
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Quantum phase transition using 
quantum walks 

4.1 Introduction 

In Chapter 2 we saw that the dynamics of the quantum walk can be controlled 
by varying the quantum coin parameters. From Chapter 3 we concluded that the 
symmetries and addition of small amount of experimentally engineered noise or 
environmental effect can be used as an additional tool to control the dynamics 
and hence the probability distribution of the quantum walk. 

These enhanced properties of the quantum walk purely due to quantum dynamics 
and small amount of engineered noise [64, 73, 74] can be used to enlarge the toolbox 
for controlling the dynamics in a physical system, to demonstrate coherent quan- 
tum control over, for example, atoms, photons or spin chain systems. Therefore, 
the use of the quantum walk to study the dynamics of particles in ID magnetic 
systems [84], atoms in optical lattice [81], and photonic Mott insulators [116] or 
to observe complex quantum phase transitions [117] and quantum annealing [48] 
would be of both theoretical and experimental relevance. In this chapter we will 
consider the use of a quantum walk to control and study the dynamics of ultracold 
bosonic atoms in an optical lattice. 
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Traditionally, theoretical studies of the dynamics of atoms in an optical lattice 
are done using mean field approaches [77] and quantum Monte Carlo methods 
[78-80]. The quantum correlation induced between atoms and position space by 
the quantum walk can serve as an alternate method for theoretical studies. In 
particular, we will consider the quantum phase transition from the Mott insulator 
(MI) - a regime where no phase coherence is prevalent - to the superfiuid (SF) - a 
regime with long-range phase coherence [81, 82] and vice versa by redistributing 
the density profile of atoms using quantum walk. The simulation of the quan- 
tum phase transition using the quantum walk occurs quadratically faster in one 
dimension (ID) compared to, varying the optical lattice depth and letting the 
atom-atom interaction follow the classical random walk behavior. Compared to 
ballistically moving the atoms, which just relocates the atoms to the new lattice 
site, quantum walk on atoms spread the wavepacket of each atom over the lattice 
sites (for example, SF regime). We will consider the discrete-time quantum walk 
with quantum coin operation -Bo,6»,o small amounts of three physically relevant 
models of noise that can be experimentally induced: a bit-flip channel, a phase-flip 
channel and an amplitude-damping channel discussed in Chapter 3 [64], to act as 
an enhanced toolbox to control the redistribution of atoms in an optical lattice. 
We will also discuss the scheme for the experimental implementation of quantum 
walk on ultracold atoms in an optical lattice. 

This chapter is organized as follows. In Section 4.2 we review phase transition in 
optical lattice using Bose-Hubbard model. The implementation of the quantum 
walk on atoms in ID MI and SF regimes is discussed in Section 4.3, where we 
analyze the dynamics and present the density profile obtained by using some of 
the properties of the quantum walk. In Section 4.4 we discuss the quantum walk 
with a noisy channel as a tool to control the redistribution of atoms. In Section 
4.5 we propose a scheme for an experimental implementation of quantum walk on 
atoms (Bose-Einstein condensates) in optical lattice and conclude with summary 
in Section 4.6. 



94 



Chapter 4. Quantum phase transition 



4.2 Quantum phase transition in optical lattice 

Bosons do not obey the exclusion principle, as temperature T — t- they all begin to 
get into the lowest energy eigenstate available, with energy e = and the effects of 
quantum degeneracy begin to emerge. That is, a large fraction of an atomic wave 
packets of bosonic atoms when cooled below a critical temperature Tc, begins 
to condense into the lowest quantum state called as Bose-Einstein condensate 
(BEC). A BEC at low enough temperature is a superfiuid described by a wave 
function that exhibits long-range coherence [118]. When the BEC is transferred 
to the lattice potential, the atoms move from one lattice site to the next by tunnel 
coupling. Bose- Hubbard model is the standard model used to discuss the dynamics 
of the boson and quantum phase transition in optical lattice. The simple possible 
approximation for the wavefunction of the many body system is a product of single 
particle state, usually known as Gross-Pitaevskii or mean-field approximation. 
Below we will briefly discuss the Bose-Hubbard model [119]. 



4.2.1 Bose-Hubbard model 

The elementary degree of freedom in this model the name implies spinless 

bosons, this takes the place of the spin-1/2 fermionic electrons in the original 
Hubbard model introduced as a description of the motion of electrons in transition 
metals. The Hamiltonian of the Bose-Hubbard model is 

(jk) j j 

J is the hopping element which allows hopping of bosons between the nearest 
neighbor pair of sites represented by (jk). bj and their Hermitian conjugate bl 
are the boson operators, which annihilates and creates bosons on the sites j and 
/c of a regular lattice in D dimensions respectively. The two operators obey the 
commutation relation 

[b,,bi] = 6,, (4.2) 
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while the two creation or annihilation operators always commute. 




(4.3) 



is the boson number operator which counts the number of bosons on each site, /ij 
is the chemical potential of the bosons, denotes the energy offset due to external 
harmonic confinement of the atoms in the jth lattice site, and U is the repulsive 
interaction between two atoms in a single lattice site. U > 0, represents the sim- 
plest possible on-site repulsive interaction between the bosons (on-site repulsion). 
Hb is invariant under global U{1) = 0(2) phase transformation under which 



We can also notice that the J term couples neighboring sites in a manner that 
prefers a state that breaks the global symmetry. However these terms compete 
with the U term in H^, which are local and prefer states that are invariant under 
symmetry transformation. Therefore we can expect a quantum phase transition 
in Hb as a function of J/U between a state in which the U{1) symmetry (4.4) is 
unbroken to one in which it is broken. As the consequence of U{1) symmetry we 
have the conservation on the total number of bosons 



and this is easily verified by noting that the Ni, commutes with H. We can notice 
that the chemical potential /i in is coupled to the conserved total number of 
bosons Nf,. In contrast, it can be noted that remains invariant under (4.4) for 
any value of /i and fi term does not break any symmetries. This gives no choice but 
to examine for all /i and can be done by using a mean- field theory approach. 

Using the mean- field theory approach, the properties of are modeled by the 
best possible sum, Hmf, of single Hamiltonians 



(4.4) 




(4.5) 




(4.6) 
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where the complex number is a variational parameter and a field to represent 
the infiuence of the neighboring sites; this field has to be self consistently deter- 
mined. This breaks the U{1) symmetry and does not conserve the total number 
of particles. As the mean-field Hamiltonian is the same on every site, the ground 
state does not spontaneously break a translational symmetry of the lattice. 

To determine the optimal value of the parameter "^b, the ground state wavefunci- 
ton of Hmf for an arbitrary \E'_b is determined. This wavefunction will simply 
be a product of single-site wavefunctions. Next the expectation value of in 
the wavefunction is evaluated by adding and subtracting Hmf by H^. Then the 
mean-field value of the ground state energy is written in the form 



where M is the number of the lattice sites, Emf{^b) is the ground state energy 
of Hjv/F, Z is the number of nearest neighbors around each lattice point (the 
coordination number), and the expectation values are evaluated in the ground 
state of HjvfF- Final step is to minimize (4.7) over variations in ^ b- By taking 
the derivative of (4.7) with respect to ^ b-, it can be show that at the optimum the 
value oi^B 



this relation, however does not hold at a general point in parameter space. 

When J = 0, the sites are decoupled, and the mean field theory is exact, \E'5 = 0, 
and we simply have to minimize the on-site interaction energy. The on-site Hamil- 
tonian involves only the operator n, and the solution involves finding the boson 
occupation number (which are integer- valued eigenvalues of fi) that minimizes H^. 
We get the ground state wavefunction 



Eq _ Emf{^b) 
M ~ M 



(4.7) 



^B = ZJ{h)- 



(4.8) 



m = no ifJ^/U)), 



(4.9) 
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where the integer- valued function no (77) is given by 



no 



0, 


for 


fi/U < 0, 


1, 


for 


< ii/U < 1, 


2, 


for 


1 < fi/U < 2, 



(4.10) 



n, for n — 1 < ji/U < n 

Each site has exactly the same number of bosons, which jumps discontinuously 
whenever fi/U goes through positive integer. When fi/U is exactly equal to a 
positive integer, there are two degenerate state on each site (with boson numbers 
differing by 1) and so the entire system has a degeneracy of 2^^. This large 
degeneracy implies a macroscopic entropy which can be lifted once we turn on 
a nonzero J. Figure 4.1 is the phase diagram of the ground state of the Bose- 
Hubbard model (4.1). As shown in the Figure 4.1, the regions with b = 
survive in lobes around J = 0, (4.9) characterized by a given integer value 
of 77-0 (yu/f/). Only at degenerate point with fi/U = integer does a nonzero J 
immediately lead to a state with \E'b 7^ 0. First we will consider the additional 

98 



Chapter 4. Quantum phase transition 



property of the lobes with = 0. The expectation value of the number of each 
site is given by 

{byb,)=nMU). (4.11) 

Existence of an energy gap and the fact that A'^f, commutes with Hb are the 
important ingredients of the above expression. First recall that for J = 0, provided 
fi/U was not exactly equal to positive integer, there was a unique ground state, 
and there was a nonzero energy gap. As a result, when we turn on a small nonzero 
J , the ground state will move adiabatically without undergoing any level crossings 
with other state. Now the J = state is exact eigenstate of Nh with eigenvalue 
MnQ{fi/U), and the perturbation arising from a nonzero J commutes with Ni,. 
Consequently the ground state will remain an eigenstate of N^, with precisely 
the same eigenvalue, MnQ{fi/U), even for small nonzero J. These regions with a 
quantized value of the density and an energy gap to all excitations are known as 
Mott insulators. The Mott insulators are also known as incompressible because 
their density does not change under changes of the chemical potential fi or other 
parameter in H^: 

The boundary of the Mott insulating phase is a second order quantum phase 
transition, i.e., a nonzero turns on continuously. 

When we turn to the phase with 7^ 0. The mean-field parameter "^b varies 
continuously as the parameters are varied. As a result all thermodynamic variables 
also change, and the density can be varied smoothly across any real positive value. 
So this is a compressible state in which 

^ # (4.13) 

The presence of 7^ implies that the U{1) symmetry is broken, and there is a 
nonzero stiffness to twists in the orientation of the order parameter. This state is 
a superfluid and that the stiffness is just a superfluid density. 

From the above description we can conclude that when J dominates the Hamilto- 
nian, the ground-state energy is minimized if the single-particle wave functions of 
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all atoms are spread out over the entire M-lattice site. If /ij = const, (homoge- 
neous system) then the many-body ground state is called the superfluid state and 
is given by 



M ^ ^ 



\Mu=o<x\^^j \V), (4.14) 

where |\^) is a vacuum state. In this state the probability distribution for the 
local occupation Uj of atoms on a single lattice site is Poissonian. The state is 
well described by a macroscopic wave function with long-range phase coherence 
throughout the lattice. With increase in the ratio U/ J, the system reaches a 
quantum critical point, the fluctuations in atom number of a Poisson distribution 
become energetically very costly, and the ground state of the system will instead 
undergo a quantum phase transition from the superfluid state to the Mott insulator 
state, a product of local Fock states of n atoms in each lattice site is given by 
[81, 82], 

M 

\i^Mi)j=o oc l[{b]r\V). (4.15) 

4.3 Quantum walk on atoms in ID Mott insula- 
tor and superfluid regime 



Lets first consider the localized atomic wave functions in the Mott insulator regime 
with one atom in each of the M lattice sites and implement the quantum walk 
on it. Atoms are first initialized into a symmetric superposition of any of the two 
internal trappable state, hyperfine levels, |0) and |1)^ at position j, 



M 




(4.16) 



The "He of each atom is spanned by the two hyperfine levels and Tip is spanned by 
the lattice site. The total system is then in the Hilbert space Tim = (n^'Hc^) ®'Hp. 
The unitary shift operation 5* (2.14) written in the form (3.1), on the above system 

7r/2 radio frequency pulse can evolve the atoms into an equal superposition of the two 
internal trappable states 
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will evolve each atom into superposition of the neighbor lattice site, establishing 
the quantum correlation between the states of the atom and the neighboring lattice 
site, 



M 
2 



^'|*A//)j=o oc [[ [ J. (4.17) 



For a system initially in state given by (4.16), t steps of quantum walk is im- 
plemented by iterating the process of S followed by the quantum coin operation 
B^^g^t^ (2.39), W = {B(^fi^c ® 1)'S', upto t times. During the iteration, the atom and 
position correlations overlap resulting in. 



M_ 

2 



{wy\'^Mi)j=o^ n {P,-M®\3-t) 



i=-M. 

■J 2 



/3,_(i+i)|0) ® \j - (t + 1)) + + (3,+t\0) ® \j + t) 

+7,_i|l) ® \j -t) + + 7,+i)|l) ® \j + t)). (4.18) 



This can be written as 



M 
2 



{Wy\<ilMi)j=o « n ( E ['^-10) + ® 1^) ' (4.19) 

/3a; and 7^ are the probability amplitudes of state |0) and |1) at lattice site x, which 
range from (j — t) to (j + t). For a quantum walk of t steps on a particle initially 
at position j = using -Bo,6»,o the quantum coin, the probability distribution 
spread over the interval (— t cos(6'), t cos(6')) in position space and ceases quickly 
outside this region [18, 61]. Therefore, after t steps the density profile of M atoms 
initially between ±^ will correlate with the position space and spread over the 
lattice site ±("1^ + tcos(6')). Figure 4.2 is the redistribution of atomic density of 
40 atoms initially in the Mott insulator state when subjected to a quantum walk 
of different number of steps with Hadamard operator H = -80,45° .0 as the quantum 
coin. In the lattice region range where all the 40 atoms are correlated can be seen 
superfluid region. 

The redistributed atoms would be in either of their internal states and this can 
be retained to study the phase transition of the two-state bosonic atoms in an 
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Figure 4.2: Density profile of the evolution of 40 atoms starting from Mott 
insulator state in correlation with the position space when subjected to a quan- 
tum walk of different numbers of steps with the Hadamard operator -Bo,45°,o 
as the quantum coin. The distribution spreads with increase in the number of 

steps. 



optical lattice [120] or all atoms can be transferred to one of the internal states. 
A technique based on adiabatic passage using crafted laser pulses can be used for 
a nearly complete transfer of population between two states [121], and a popular 
example of one such technique is stimulated Raman adiabatic passage (STIRAP) 
[122]. 

Along with the correlation between the atoms and the position space, (4.18) also 
reveals the overlap of the probability amplitude of different atoms in the position 
space. That is, during each step of the quantum walk, the amplitude of the 
states of each atom overlaps with the amplitude of the states of the atoms in the 
neighboring lattice site. After the iteration of the t steps > the overlap of 

all atoms could be seen at the central region of the lattice. 

When the number of steps of the quantum walk is equal to the number of lattice 
site (t = M), the overlap of the fraction of amplitude of all M atoms exists 
within the lattice sites -|-:?££2£M making the region identical to the superfluid 
state. Beyond -|-ill£££M the number of atoms in correlation with that position 
space decreases and hence the number of overlapping atoms also decrease. To 
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make all M atoms to spread between lattice site with the usual method of 
lowering the optical potential depth following the classical random walk protocol 
takes steps. Therefore, we can conclude that using quantum walk a long range 
correlation can be induced quadratically faster than any other technique using the 
classical random walk protocol. 

1.5 
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Figure 4.3: Density profile of the evolution of 40 atoms starting from super- 
fluid state in correlation with the position space when subjected to a quantum 
walk of different numbers of steps with the Hadamard operator -Bo,45°,o as the 
quantum coin. The distribution spread with increase in the number of steps 
and is almost uniform between the lattice site ±20 when t = 25. At this stage 
the optical potential depth can be increased to cancel the correlation and obtain 

the Mott insulator state. 

Similarly, Figure 4.3 is the redistribution of atomic density of 40 atoms initially 
in the superfluid states when subjected to the quantum walks of different number 
of steps using Hadamard operator -80,45° ,o the quantum coin. To reach Mott 
insulator state the distribution should be uniform over the lattice site without cor- 
relation. A distribution that is approximately uniform within the region ±^ can 
be obtained using the quantum walk. In Figure 4.3, when t = 25 the distribution 
is almost uniform between the lattice site ±20 retaining the correlation with the 
position space. Once the distribution is uniform, the optical potential depth can 
be increased to cancel the correlation and obtain the Mott insulator state. The 
uniformity of the distribution can also be improved by introducing a noise channel, 
as we show later. 
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The variance and the probabihty distribution can be controlled using the param- 
eters 6, C, and ( in -B^,^,^- Figures 4.4 and 4.5 shows the density distribution 
obtained by implementing the quantum walk on atoms in Mott insulator and su- 
perfluid state {t = M = 40) respectively with different values of 6 in the coin 
operator -Bo,0,o- The spread is wider for 6 = 30° and decreases with increase in 6. 
In Figure 4.5 the distribution is almost uniform for 6 = 60° between the lattice 
site ±20. At this value the optical potential depth can be increased to cancel the 
correlation and obtain the Mott insulator state. 




Lattice site (j) 



Figure 4.4: Distribution of atoms initially in Mott insulator state when sub- 
jected to quantum walk (t = M = 40) using different values for 9 in the operator 
Bqq Q. The spread is wider for 9 = 30° and decreases with increasing 9. 



The Hamiltonian of the system in general can be described by the ID Bose- 
Hubbard model for two-state atoms, 

{j,k) j,a='[,l 

+ Yl (^^%-i)t^.t + dRb\^+iJ^,) , (4.20) 
j 
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Figure 4.5: Distribution of atoms initially in superfluid state when subjected 
to a quantum walk {t = M = 40) using different values for 9 in the coin operator 
Bo,e.o- The distribution is almost uniform for = 60° between the lattice site 
±20. At this value the optical potential depth can be increased to cancel the 
correlation and obtain the Mott insulator state. 



here f and l represent the terms for atoms in state |0) and state |1) respectively. 
U is the interaction between atoms in state |0) and V^i^i) is the interaction 
between atoms in same state, di and are the left and right displacement 
terms in the Hamiltonian operation. The Hamiltonian is evolved with coin toss 
operation in a regular interval of time r, the time required to move the atom to the 
neighboring site. The optical lattice can be dynamically manipulated to evolve 
atoms in a superposition of lattice site without giving time for the atom-atom 
interaction, and the optical potential depth of the system can be configured just 
above the level where there is no direct tunneling. Then the dynamics of atoms 
in lattice will be dominated by the last term of the Hamiltonian in (4.20) ignoring 
the atom-atom interaction. Therefore scaling up the scheme of using the quantum 
walk on systems with large number of atoms in each lattice site or to infinitely large 
numbers of lattice site is also straight forward; whereas, the atom-atom interaction 
play a prominent role during the usual method of varying the potential depth. 
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4.4 Quantum walk with a noisy channel as a 
toolbox 

We saw the effect of quantum walk with different number of steps and the coin 
parameters 6 on atoms in optical lattice in the previous section. To demonstrate 
the effect of the quantum walk with a noisy channel on atoms in optical lattice 
we consider a bit-flip channel, a phase-flip channel and an amplitude-damping 
channel. The bit-flip channel flips the state of the particle from |0) (|1)) to |1) 
(|0)) (Pauli X operation) and a phase- flip channel flips the phase of |1) to — 11) 
(Pauli Z operation). We use the notation p for the noise level where, < p < 1. 
Therefore, the bit- (phase-)flip channel flips the state (phase) with probability p 
during each step of the quantum walk. An amplitude-damping channel leaves state 
|0) unchanged but looses the amplitude of state |1) with probability p resulting 
in an asymmetric distribution [5]. From Chapter 3 we know that the maximum 
decoherence effect using bit- and phase- flip channel is at p = 0.5 due to symmetries 
induced by these two channels during the quantum walk; whereas, the amplitude- 
damping channel does not obey any symmetry and hence the maximum effect is 
for p = 1 [64]. Our numerical implementation of these channels evolves the density 
matrix employing the Kraus operator representation. 

Figure 4.6 is a comparison of the distribution obtained without noise channels to 
the distributions with phase-flip {p = 0.02 and p = 0.1) and amplitude-damping 
channel {p = 0.2) on atoms initially in the Mott insulator state. Similarly, the 
redistribution of atoms in superfluid state when subjected to quantum walk with- 
out and with noise channels is shown in Figure 4.7. In Figure 4.7, with the phase 
flip noise of p = 0.02, a uniform distribution is obtained. Then the optical poten- 
tial depth can be increased to cancel the correlation and obtain the Mott insulator 
state. Other plots shows the different redistribution of atoms that can be obtained 
using different noise channels. 

The effect of bit-flip on the distribution is close to the one obtained using phase- 
flip channel and hence the redistribution in the Figures 4.6 and 4.7 are shown only 
for phase-flip and amplitude-damping channels. Increasing noise level affects the 
variance of the quantum walk [64], which proportionally affect the atom-position 
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Figure 4.6: Atoms in Mott insulator state after implementing the quantum 
walk (t = M = 40,6 = 45°) with noise channel. With increased phase damping 
from pp = 0.02 to Pp = 0.1 the distribution (nj) at the central region gets closure 
to Gaussian. Amplitude damping of state |1) followed by a bit flip at each step 
introduces asymmetry to the distribution, pa = 0.2 (without a bit flip the shift 

would have been to the left). 



correlation and the atom-atom overlap region. Therefore, higher noise level classi- 
calize the dynamics and hence, it is important to restrict the noise level to a very 
low value (p ^ 0.1) to use noisy channel as a useful tool. 

Some of the distribution presented in Sections 4.3 and 4.4 using quantum walks 
are visibly similar to some of the distribution presented using different technique; 
in [123] using time evolution density matrix renormalization group (t-DMRG) 
and in [124, 125] using quantum Monte Carlo simulation. Figure 4.8 and Figure 
4.9. These technique present the change in density profile distribution with time, 
whereas in this chapter we have discussed the change with number of steps of the 
quantum walk. 

In this section, we have shown the use of the quantum walk to study the dy- 
namics of atoms in an optical lattice and expedite the process of quantum phase 
transition. We have also used the quantum walk with experimentally realizable 
noisy channels to show the additional control one can have over the evolution and 
atomic density redistribution. Theoretically, the evolution of the density profile 
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Figure 4.7: Atoms in superfluid state after implementing a quantum walk 
(t = M = 40,9 = 45°) with noise channel. With noiseless quantum walk, 
the atoms spread, getting close to a uniform distribution. For phase damping 
Pp = 0.02 the distribution is uniform between ±20 and gets closer to Gaussian 
at Pp = 0.1. Amplitude damping of state |1) followed by a bit flip introduces 
asymmetry to the distribution, pa = 0.2 (without a bit flip the shift would have 

been to the left). 
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Figure 4.8: Density profile of the evolution of atoms initially in Ml state from 
other references, (a) Evolution of 100 atoms with time for atoms initially in MI 
state [125]. (b) Evolution of 20 atoms with time (the dynamics is driven by the 
potential depth and tunneling of atoms between the lattice) [123]. These profile 
have visibly similar distribution compared to Figure 4.2, Figure 4.4 and part of 

Figure 4.6 in this chapter. 
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Figure 4.9: Density profile of the evolution of 40 atoms initially in SF state 
from [124]. The initial SF profile in om' plots and in [124] are different and 
hence the similarity are not as close as we have for atoms in MI state). The 
profile in Figure 4.3 after 25 steps of quantum walk is similar to the Profile at 
TG limit (Tonks-Girardeau limit - leading to MI state) in the above figure. 



with a quantum walk can be used in place of quantum Monte Carlo simulation to 
study the correlation and redistribution of atoms in optical lattices. We expect 
the quantum walk to play a wider role in simulating and expediting the dynamics 
in various physical systems. 



4.5 Implementation 



Experiments on ultra cold atoms has been one of the active areas of research 
in the last two decades. Optical lattices ranging from a simple periodic, square 
and cubic, to a more exotic ones, such as hexagonal and Kagome lattices using 
superlattice technique (see for example, [126-132]) have been created to trap and 
manipulate ultra cold atoms. Manipulation of cold atoms in time-varying optical 
lattice has also been reported [133, 134]. This has provided flexibility in designing 
and studying quantum phases and quantum phase transitions in coherent and 
strongly correlated ultracold atoms. 
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Various schemes have been proposed to implement quantum walk on neutral atoms 
in an optical lattice [55, 56]. In [135], the controlled coherent transport and split- 
ting of atomic wave packets in spin dependent optical lattice has been experimen- 
tally demonstrated using rubidium atoms. A BEC of up to = 3 x 10^ atoms is 
initially created in a harmonic magnetic trap. A three dimensional optical lattice 
is superimposed on the BEC and the intensity is raised in order to drive the system 
into a Mott insulating phase [82]. 

Two of the three orthogonal standing wave light fields is operated at one wave- 
length, Xy^z = 840 nm and the third along the horizontal direction is tuned to the 
wavelength A^^ = 785 nm between the fine structure splitting of the rubidium Dl 
and D2 transitions. Along this axis a quarter wave plate and an electro-optical 
modulator (EOM) is placed to allow the dynamic rotation of the polarization vec- 
tor of the retro-reflected laser beam through an angle 6 by applying an appropriate 
voltage to the EOM. After reaching the Mott insulating phase, the harmonic mag- 
netic field is completely turned off but a homogeneous magnetic field along the x 
direction is maintained to preserve the spin polarization of the atoms. The light 
field in the y and z direction is adiabatically turned off to reduce the interaction 
energy, which strongly depends on the confinement of the atoms at a single lattice 
site. 

A standing wave configuration in the x direction is used to transport the atoms. By 
changing the linear polarization vector enclosing angle 6, the separation between 
the two potentials is controlled. By rotating the polarization angle 6 by vr, with 
the atom in a superposition of internal states, the spatial wave packets of the atom 
in the |0) and the |1) state are transported in opposite directions. The final state 
after such a movement is then given by l/v^(|0,x — 1) + zexp(z/3p)|l, x -|- 1)). 
The phase /3p between the separated wave-packets depends on the accumulated 
kinetic and potential energy phases in the transport process and in general will be 
nonzero. The coherence between the two wave-packets is revealed by absorption 
imaging of the momentum distribution. A 7r/2 microwave pulse is applied before 
absorption imaging to erase the which-way information encoded in the hyperfine 
states. 
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However, to increase the separation between the two wave-packets further, one 
could increase the polarization angle 6 to integer multiples of 180°. To overcome 
the limitation of the maximum voltage that can be applied to the EOM, a vr pulse 
after the polarization is applied, thereby swapping the roles of the two hyperfine 
states. The single particle phase Pp remains constant throughout the atomic cloud 
and is reproducible. After the absorption imaging a Gaussian envelope of the 
interference pattern is obtained. 

One can build upon the above technique to implement a quantum walk. The above 
setup can be modified by dividing the separations (splitting) into small steps and 
introducing a rotation (vr/2 pulse for Hadamard rotation) after each separation 
without intermediate imaging. The absorption imaging of the distribution of the 
atomic cloud after t steps would give the interference pattern similar to the prob- 
ability distribution of the quantum walk. 

This effect of the addition of phase during the quantum walk process can be 
easily understood from the phase damping channel and arbitrary phase rotation 
discussed in Chapter 2. The addition of vr pulse to overcome the limit of EOM is 
a bit flip operation in the quantum walk. 

In the following section we propose a scheme to implement a discrete time quantum 
walk on ultra cold atoms in BEC state. BEC being a macroscopic wave packet, 
the scheme we propose can implement a quantum walk retaining the macroscopic 
features of the wave packet [57] or the scheme can be accordingly modified to 
implement the walk at an individual atom level. For the transition from superfiuid 
to the Mott insulator and vice versa using quantum walk, we need to consider the 
quantum walk on individual atom. 

To implement a macroscopic quantum walk, atoms in BEC state is first evolved to 
macroscopic superposition state and a stimulated Raman kick, two selected levels 
of the atom are coupled to the two modes of counterpropagating laser beams to 
coherently impart a translation of atoms in the position space. After each transla- 
tion, the wave packet is again evolved into the macroscopic superposition state and 
the process is iterated to implement large number of steps of macroscopic quan- 
tum walk. With a certain modification to this scheme, that is, by evolving atoms 
to the superposition of the states at individual atom level and implementing the 
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shift operation before the interatomic interaction takes over to form a macroscopic 
superposition, the quantum walk at individual atom level can be realized. 

4.5.1 Quantum walk using Bose-Einstein condensate 

As discussed in Section 2.3, two degrees of freedom, the coin Hilbert space He and 
the position Hilbert space Hp are required to implement the discrete-time quantum 
walk. The state of the BEC formed from the atoms in one of the hyperfine states |0) 
or |1) can be represented as \Obec) or \1bec) (macroscopic states of condensed 
atoms). The BEC formed is then transferred to an optical dipole trap with long 
Rayleigh range 2/?^. With the appropriate choice of power and beam waist, uq^ of 
the trapping beam, the BEC can remain trapped at any point within the distance 
±z from the focal point in the axial direction of the beam (See Appendix C for 
calculations) . 

The position of the BEC formed in one of the Hyperfine state \Obec) or \1bec) 
and trapped at the center of the optical trap is described by a wave packet I^E'jo) 
localized around a position jo, i.e., the function {jl'^jo) corresponds to a wave 
packet centered around jo- The BEC formed in one of the Hyperfine state is 
evolved into the macroscopic superposition (Schrodinger cat) state 

l^.n) = {a\OBEc) + b\lBEc)) ® |^,o) (4-21) 

by applying the rotation (coin operation B), where |ap + |6p = 1. We will discuss 
this process in detail in Section 4.5.1.1. The coin Hilbert space Tic can then be 
defined to be spanned by the two internal trappable macroscopic states \Obec) 
and \1bec) of the BEC. The position Hilbert space Hp is spanned by the positions 
within the long Rayleigh range optical trap (trappable range). Once the BEC is in 
the superposition state, a unitary shift operator, stimulated Raman kick, Section 
4.5.1.2, corresponding to one step length / in the form (3.2) is applied; 

S' = {X® I) exp{-i{\OBEc){OBEc\ - \UEc)k^BEc\) ® PI). (4.22) 

^Distance at which the diameter of the laser beam size increases by a factor of \f2. 
■^MinimmTL radius of the beam, at the focal point. 
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P being the momentum operator. Step length / is chosen to be less than the 
spatial width of the condensate wave packet I^E'jo). 

The application of the unitary shift operator (5") on the wave function of (4.21) 
entangles the position and the coin space and implements the quantum walk. 
Stimulated Raman kicks evolves the state to, 

= {a\lBEc) ® e-^" + b\OBEc) ® e*") (I ® |^,o)) , (4-23) 

where the wave packet is centered around jo i ^- Note that the values in the coin 
space have been flipped. This can be corrected by applying a compensating bit 
flip on the BEC, S = {X (g> T)S'. From the symmetries of quantum walk discussed 
in Chapter 3 we know that the distribution is invariant to bit flips at each step and 
therefore, augmenting S' with the X operation is unnecessary. To realize a large 
number of steps of the macroscopic quantum walk, the process of shift operation 
followed by the quantum coin operation (iJ^^e^^ (8)1)5" is iterated without resorting 
to intermediate measurement in the long Rayleigh-range optical dipole trap. 

However, the coin operation B^^qi^ can be performed such that, its action will 
evolve the individual atoms into the superposition of |0) and |1) and not the 
superposition of macroscopic state \Obec) and \Ibec), Section 4.5.1.1. This can 
be used to drive the transition between the superfluid and the Mott insulator state. 
To realize transition from the superfluid to the Mott insulator phase, atoms are 
first redistributed by implementing the quantum walk using the rotation operation 
(coin operation) i?5,6i,c on individual atoms in a long-Rayleigh-range optical trap 
and the optical lattice is switched on at the end to confine atoms and cancel 
the correlation. For transition from Mott insulator to the superfluid, the atoms 
initially in the optical lattice are transferred to long-Rayleigh-range dipole trap 
and quantum walk is implemented. 

4.5.1.1 Macroscopic cat state of Bose-Einstein condensate 

Various schemes have been proposed for producing macroscopic superposition or 
Schrodinger cat states in BECs [136-138]. The scheme described in [136], shows 
that the two species (two interacting Bose condensates) can be evolved to the 
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macroscopic superposition of the internal atomic states. In this scheme, atom- 
atom interactions are mediated through atom-atom colhsions and a Josephson-hke 
laser coupling that interchanges internal atomic states in a coherent manner. In 
certain parameter regimes the ground state of the Hamiltonian is a superposition 
of two states involving a particle number imbalance between the two condensates. 
Such a state produced by the normal dynamic evolution of the system represents 
a superposition of two states which are macroscopically (or mesoscopically) distin- 
guishable, and hence can be called a Schrodinger cat state. But, two interacting 
Bose condensates are needed for this scheme. In [137], it is shown that a macro- 
scopic superposition state can be created by a mechanism involving the coherent 
scattering of far-detuned light fields and it neglects the collisional interactions be- 
tween particles. The major drawback of the above schemes is that the time needed 
to evolve to a cat state can be rather long, and thus problems due to decoherence 
would be greatly increased. An other scheme proposed in [138] involves the cre- 
ation if macroscopic superposition by an adiabatic transfer of the ground state 
of the Josephson-coupling Hamiltonian, that is, after the initial state preparation 
the Josephson coupling is turned on for some amount of time and turned off. The 
resulting modified quantum state is a Schrodinger cat state. But, the production 
of such a state involves considerable experimental difficulty. 

The ideal scheme to demonstrate the quantum walk is to confine the BEC that 
has an attractive interaction between atoms in two hyperfine levels |0) and |1) in 
a single optical potential well. The rotation operation is applied on the BEC in 
the potential well to transfer (or rotate) the atoms part of way between states 
|0) and |1) using a resonant rf pulse ^ of duration r and detuning A from the rf 
resonance [139]. The resonance rf pulse couples the atomic hyperfine states |0) 
and |1) with a coupling matrix element huifi/2, where is the Rabi frequency 
and the duration of pulse is much shorter than the self-dynamics of the BEC. The 
amplitude of these states evolves according to the Schrodinger equation, 

"'Microwave pulses are also used. 
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At this stage, each atom evolves into the superposition state a(r)|0) +6(r)|l), with 
a(r) = cos{ujrt/2) and 6(r) = sin(a;ijr/2) for detuning A = 0. The A^-particle 
wave function of the BEC is a product of the single-particle superpositions of |0) 
and |1), that is, it is still a microscopic superposition and is given by 

[a{r)\0) + bir)\l)f = g J 1 a(r)^-"KrniV - n,n), (4.25) 

where \N — n,n) is the state with N — n atoms in state |0) and n atoms in state 
|1). The individual atoms in the superposition state interact among themselves. 
Interatomic interactions, which provide nonlinear terms through binary collision 
as seen from the viewpoint of single-particle dynamics helps in generating highly 
entangled many body states [140, 141]. That is, the Hamiltonian that governs the 
BEC with its attractive inter-atomic interactions, after time ri, (4.25) evolves into 
the macroscopic superposition state in which all atoms are simultaneously in level 
|0) and level |1) [136, 142, 143], 

[a(ri)|0) + h{n)\l)f = a(n)^|iV,0) + 6(n)^|0,iV) (4.26) 

where |A^, 0) = \Qbec) and |0,A^) = \1bec)- Symmetric probability distribution 
aiji) = b{ri) = 4= can be obtained by carefully choosing ur and r. 



Microscopic superposition state of atoms : 

Equation (4.25) gives the microscopic superposition state of A^-atoms. Implement- 
ing the unitary shift operator before the individual atoms in the superposition 
states start interacting will prevent the evolution of atoms to macroscopic super- 
position state. Therefore we can choose between the macroscopic, Schrodinger cat 
state or the superposition of individual atoms, microscopic superposition state. 
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4.5.1.2 Unitary operator - Stimulated Raman kicks 

A unitary shift (controlled-shift) operation S is applied on the BEC to spatially 
entangle the position and the coin space and implement the quantum walk. Various 
schemes have been worked out to give momentum kick to ultracold atoms in a 
trap [143, 144]. A technique was reported in [143] where coherent rf-induced 
transitions were used to change the internal state of the atoms in the magnetic trap 
from a trapped to an untrapped state and thus displacing the untrapped atoms 
from the trap. This method, however, did not allow the direction of the output- 
coupled atoms to be chosen. Later, controlling the direction of the fraction of the 
outgoing BEC using a stimulated Raman process between magnetic sublevels was 
experimentally demonstrated [144]. This technique was used to extract sodium 
atoms lasers from the trapped BEC. 

A stimulated Raman process can also be used to drive the transition between two 
optically trappable states of the atom |0) and |1) using the virtual state |e) as 
an intermediary state, and impart a well-defined momentum to spatially translate 
the atoms in the coherent state (BEC). A unitary shift operation S' thus can be 
applied on the atom in the BEC using a stimulated Raman process. A pair of 
counterpropagating laser beams 1 and 2 (Figure 4.10) with frequency ui and U2 
and wave vectors ki and k2 is applied on the BEC for a 2A^-photon transition time 
(A^ is the number of atoms in the BEC) to implement one unitary shift operation. 
These beams are configured to propagate along the axial direction of the optical 
dipole trap. A stimulated Raman transition occurs when an atom changes its 
state by coherently exchanging photons between the two laser fields, absorption 
of the photon from laser field 1 and stimulated emission into laser field 2 or by 
absorption from field 2 and stimulated emission into field 1. The BEC initially in 
eigenstates of the atom |0) (|1)) can absorb the photon from field 1 (2) and re-emit 
the photon into field 2 (1). This inelastic stimulated Raman scattering process 
imparts well-defined momentum on the coherent atoms, 

P = h{ki - fca) = h6z (4.27) 
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Figure 4.10: Light field configuration for a stimulated Raman transition pro- 
cess to give directional momentum kick. A is the detuning of the laser from its 
transition frequency. \S) signifies (l/\/2)(|0) ± |1)). 



to the left during |0) — ^ |e) — ^ |1) and 

P = h{k2 - h) = -hSz (4.28) 
to the right during |1) ^ |e) ^ |0), where 

\ki - = \k2- ki \ = 6, 

Thus, conditioned to being in coin state |0) (|1)), the atoms in the BEC receive a 
momentum kick (shift) h6z {—h6z). This process of imparting momentum is cahed 
stimulated Raman kick and can be analyzed as a photon absorption and stimulated 
emission between three bare states |0), |1), and |e) driven by two monochromatic 
light fields. The matrix element for photon absorption and stimulated emission 
during each stimulated Raman kick can be written as 

(rifc, - l,nfc2 + 1, l\HjJnk^,nk^,0) (4.29) 

during |0) — |e) ^ |1) and 

(rifci + l,rifc2 - l,0\HiJnk,,nk^, 1) (4.30) 
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during |1) -A |e) |0). Hj^ and Hj^ are the interaction Hamiltonians (electric 
dipole Hamiltonians). The field causing the transitions between bare states |0) and 
|e) is detuned from resonance by Ai and has a constant dipole matrix element Vi 
and phase (pi. The field causing transitions between bare states |e) and |1) is 
detuned from resonance by A2 and has a constant dipole matrix element V2 and 
phase ip2- The time-independent Hamiltonian in the rotating- wave approximation 
for this system can be written in terms of projection operations as 

Hj^ = hA,\e){e\ + ^Ai + A2)|l)(l| - ^[|e)(0| exp(-z^i) 

+ |0)(e|exp(z(pi)] - ^[|e)(l|exp(z<^2) + |l)(e| exp(-zv^i)] (4-31) 

In the above interaction picture, the energy of the bare state |0) is chosen to be 
zero. 

Hj^ = hA2\e){e\ + h{A, + A2)|0)(0| - :^[|e)(l| exp(-zv^2) 

hv 

+ |l)(e| exp(z(^2)] - ^[|e)(0| exp(z(^i) + |0)(e| exp(z(^i)]. (4.32) 

In the above interaction picture the energy of the bare state of the atom |1) is 
chosen to be zero. Vi = ^' , for i = 1,2. fii is the dipole operator associated 
with the states |0),|e) and |e),|l) and Eq is the electric field of the laser beam. 

With the appropriate choice of V, A, and (f, the probability of being in bare state 
|1) (or |0) depending on the starting state) after time t can be maximized to be 
close to one, where t is the time required for one stimulated Raman kick. As 
discussed in Section 4.5.1.1, after every stimulated Raman kick on atoms due to 
the attractive interaction between atoms, they evolve into the BEC. 

The coin state of the atoms in the BEC |0) (|1)) after stimulated Raman kick fiips 
the states to |1) (|0)) and symmetries due to bit fiip discussed in Chapter 3 will 
take care of these fiip without affecting the distribution. 
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4.5.1.3 Physical setup for the implementation 

Magnetic trap technique [145] played a major role in the first formation and early 
experiments on the BEC. But the magnetic trap has a limitation, it can trap and 
manipulate atoms only of certain sublevels. For example, If we consider '^'^Rb atom 
as an example, the \F = 2, m/ = 2) state can be confined in the magnetic trap 
whereas \F = l,mf = 1) cannot be confined using the magnetic trap [146]. Under 
appropriate conditions, the dipole trapping mechanism works independent of the 
particular sub-level of the electronic ground state. Thus, internal ground state can 
thus be fully exploited using an optical dipole trap technique and be widely used 
for various experiments [147]. Today the BEC from bosonic atoms has been very 
consistently formed and manipulated using various configurations of magnetic and 
optical traps, an all optical dipole trap technique has also been developed [148]. 

A BEC first formed in one of their internal states (eigenstates) of atoms |0) or 
|1), using any of the techniques [145, 148] mentioned above, is transferred to a 
far detuned, long Rayleigh range {zpi) optical dipole trap. A sizable number of 
steps of the quantum walk can be implemented within the axial range ±z without 
decoherence in a long Rayleigh-range trap. Appendix C has calculated numerical 
values of potential depth, power required to trap ^'^Rb atoms at distance z = Xn 
from the focal point of the trapping beam (after compensating for gravity) using 
light fields of various frequency and beam waists Uq. In the same way one can work 
out the required power and beam waists to trap ultracold atoms at a distance z 
from the focal point for different species of atoms using laser fields of different 
detuning from the resonance. 

Once the BEC is transferred into the optical dipole trap, a resonant rf pulse (rota- 
tion) of duration r and detuning A is applied to make it evolve into the Schrodinger 
cat state (or a microscopic superposition state depending on the requirement) af- 
ter time Ti. The Schrodinger cat state for the ^^Rb BEC trapped in one of the 
states |0) = |-F = 1, my = 1) or |1) = |F = 2, = 2), can be realized by applying a 
radio frequency (rf) pulse, fast laser pulse, a standard Raman pulse or microwave 
techniques. 

The BEC in the superposition state is then subjected to a pair of counterpropa- 
gating beams to implement the stimulated Raman kick for duration t^, 2N photon 
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Figure 4.11: Physical setup to implement the quantum walk using the BEC. 
In the figure a is a BEC, b is a dipole trap with a long Rayleigh range, and c is 
the counterpropagating laser beam, ui and 002 used to implement the stimulated 
Raman kick, unitary operator. A coin rotation and Raman pulse are applied 
throughout the trap region. The stimulated Raman pulse and the rf pulse are 
applied alternatively to realize the quantum walk. 




Figure 4.12: Multiple microtraps which can be used to confine the atom 
observe the distribution of atoms after implementing a considerable number of 
steps of the quantum walk on the BEC. 

transition time (A^ being the number of photons). After implementing the stimu- 
lated Raman kick, the BEC is left to translate for a duration P/ml, the time taken 
by the BEC to move distance / during or after which the rf vr pulse is applied. 
This process of applying a rf pulse (coin rotation) and a stimulated Raman kick 
(shift), followed by a compensatory rf pulse (bit flip) is iterated throughout the 
trapping range to realize a large number of steps (Figure 4.11). 

Measuring the probability distribution : 

After t steps of the quantum walk the superposition in position space is made 
to collapse by applying multiple microtraps to confine and redistribute atoms in 
the position space. The multiple microtraps with an equal spacing / between each 
potential well can be created using the known techniques [133, 149] (Figure 4.12). 
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The microtraps are switched on, turning off the long Rayleigh range optical trap 
simultaneously after a time interval of T. 

T = ttd+ [t -l)Ti+tP/ml. (4.33) 

T is the time taken by the BEC to travel a t step quantum walk each of distance 
/ in a line. The time is the 2A^-photon transition time required to implement 
one unitary shift operation, r is the time duration to bring the BEC to the su- 
perposition state, P/ml is the time taken by the BEC to move one step distance 
/ with momentum P. Fluorescence measurement is performed on the atoms in 
the microtrap to identify the final position/distribution of atoms in the position 
space. 

If the quantum walk was implemented on a BEC evolving the wave packet to the 
macroscopic superposition during each coin operation, the distribution of atoms 
will be similar to the distribution of single particle quantum walk. We should note 
that by retaining the macroscopic behavior during the walk we get the quantum 
walk distribution in a single measurement at the end. Where as the experiment 
would have to be repeated a large number of times if we intend to get a probability 
distribution of the quantum walk using a single particle. 

If the quantum walk is implemented only by evolving the atoms to the superposi- 
tion at individual atom level during each coin operation we get the distributions 
similar to the ones presented in the earlier sections of this chapter. 

Decoherence and physical limitations : 

The decoherence of the BEC also leads to the decoherence of the quantum walk. 
When the atoms in the trap are not coherent they no longer follow the coherent 
absorption and stimulated emission of light. Some atoms absorb light field 1 and 
emit light field 2 and bring atoms to state |1) and some absorb light field 2 and 
emit light field 1 and bring atoms to state |0), displacing atoms in both directions 
in space giving no signature of displacement in the superposition of position space. 
This will contribute to collision and heating, and finally, atoms escape out of the 
trap. The number of implementable steps of quantum walk without decoherence 
depends mainly on (a) the Rayleigh range, as the BEC moves away from the 
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trap center the width of the wavepacket increases and contributes to the internal 
heating of the atoms. Beyond a certain distance Xn from the dipole trap, center 
atoms in the trap decohere resulting in the collapse of the quantum behavior, (b) 
The stimulated Raman kick and rf pulse used to implement the quantum walk 
also contribute to the internal heating and decoherence of the atoms after a few 
iterations. With the careful selection of beam waist and laser power one can have 
a trap with a long Rayleigh range 2;^ (Appendix C) so that the one- dimensional 
quantum walk can be implemented in a line of close to one centimeter length. 
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4.6 Summary 



• We proposed the use of the quantum walk to redistribute atoms, and studied 
its dynamics in an optical lattice and to expedite the process of quantum 
phase transition. We have demonstrated the coherent control over the atoms 
using the coin degree of freedom during the evolution of the walk. We have 
also used an experimentally realizable noisy channels to show the additional 
control over the evolution and atomic density redistribution. Theoretically, 
the evolution of the density profile with a quantum walk can be used in place 
of quantum Monte Carlo simulation or the time evolution density matrix 
renormalization group (t-DMRG) to study the correlation and redistribution 
of atoms in optical lattices. We expect the quantum walk to play a wider 
role in simulating and expediting the dynamics in various physical systems. 

• We proposed a scheme to implement a discrete-time quantum walk on Bose- 
Einstein condensate (BEC). BEC being a macroscopic wave packet, the 
scheme we propose can implement a quantum walk retaining the macroscopic 
features of the wave packet or the scheme can be accordingly modified to 
implement the walk at an individual atom level. To realize the quantum 
phase transition from superfluid to the Mott insulator state and vice versa, 
implementation of walk on individual atom level will be effective. 



123 



Chapter 5 



Spatial entanglement using 
quantum walk on many body 
system 



5.1 Introduction 



Entanglement in many body systems has been more than just a computational 
resource, it has been used as a signature of quantum phase transition [83-85]. 
Therefore, in the last few years, the study on entanglement in many body system 
has been one of important area of research interfacing between condensed matter 
system and quantum information sciences [150]. 

Measure of entanglement in a pure bipartite system is a function of the eigen- 
values of the reduced density matrix. However, as the number of particles in the 
system increases, the complexity of finding the appropriate entanglement measures 
also increases making scalability of entanglement measure an enormous task. To 
address the scalability problems, a scalable entanglement measures which do not 
diverge with the system size were proposed and global entanglement measure, a 
polynomial measure of multipartite entanglement is one such measure [86]. Using 
this global entanglement measure, namely the Meyer- Wallach measure, we investi- 
gate the evolution of spatial entanglement- particle-number entanglement between 
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regions of space in a many distinguishable particle system subjected to quantum 
walk process. Spatial entanglement has also been explored using different meth- 
ods, for example, in an ideal bosonic gas it has been studied from off-diagonal 
long-range order [151]. For our investigations of spatial entanglement using multi 
particle quantum walk, we consider distinguishable particles whose dynamics can 
be controlled using the quantum coin parameters, initial state of the particles, 
number of particles in the system and the number of steps of quantum walk. In 
particular we consider distinguishable particles in one dimensional open and closed 
chains. The spatial entanglement thus created can be used, for example to create 
entanglement between distant atoms in optical lattice [152] or as a channel for 
state transfer in spin chain systems [153-155]. At this stage, we find calculating 
spatial entanglement for indistinguishable particle case computationally enormous 
time consuming. 

This chapter is organized as follows. In Section 5.2 the entanglement between the 
particle and the position space has been discussed. In the same section we also 
introduce spatial entanglement and discuss spatial entanglement using single and 
many particle quantum walk. In Section 5.3, we present the measure of the spatial 
entanglement on the system using Meyer- Wallach global entanglement measure 
scheme for particles in one dimensional lattice and in a closed chain (ra— cycle). We 
also demonstrate the control over the entanglement using the dynamical properties 
of quantum walk. We conclude with a summary in Section 5.4. 



5.2 Entanglement 

5.2.1 Position-particle entanglement 

Quantum walk entangles the particle Hilbert space Tic and the position Hilbert 
space Hp. As discussed in Chapter 2, He is spanned by the basis state (internal 
state) of the particle |0) and |1) and Tip is spanned by the basis state of the 
position where j G Z. Lets consider a discrete-time quantum walk on a 

particle initially in state given by symmetric superposition state (2.30) at position 
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Figure 5.1: Entanglement of single particle with position space when subjected 
to quantum walk. The initial state of the particle is ;^(|0) + and is evolved 
in position space using different values for 9 in the quantum coin operation 
-Bo,6»,o- The entanglement initially oscillates and approaches asymptotic value 
with increase in number of steps. For smaller values of 9 the entanglement 
is higher and decreases with increase in 9. Initial oscillation is also larger for 

higher 9. 



j with coin operation 

_ / cos(^) sin(^) 

Bq e — \ 

■ ' \^ sm{e) - cos(^) 

After the first step, Wg = S'(i?o,e,o ® 1); the state takes the form 

= 7(|0)®|V^,-i)) + 5(|l)®|^,+i)) (5.2) 

where 7 = (^£2£W+pMj and 5 = (^EBMzpWj. Schmidt rank of is 2 
which implies entanglement in the system. The value of entanglement can be 
further quantified with increase in number of steps by computing the Von Neumann 
entropy of the reduced density matrix of the position subspace. Position-particle 
entanglement in quantum walk on regular graph has been studied in detail in [156]. 

Figure 5.1 is the plot of the entanglement against the number of steps of the 
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Figure 5.2: Entanglement of single particle with position space when subjected 
to quantum walk. The initial state of the particle is given by form (2.11) with 
6 = 40° and rj = 30° and is evolved in position space using different values for 
9 in the quantum coin operation Bq q q. The entanglement initially oscillates 
and approaches asymptotic value with increase in number of steps. For smaller 
values of the entanglement is higher and decreases with increase in 9. Initial 
oscillation is also larger for higher 9. 

quantum walk on particle initially in symmetric superposition state using different 
values for 6 in the operation Wg. Von Neumann entropy of the reduced density 
matrix of the coin is used to quantify the entanglement between the coin and the 
particle position, 

Ecit) = -Y,^M2i\) (5.3) 

j 

where are the eigenvalues of the reduced density matrix of the coin after t steps 
(time). The entanglement initially oscillates and reaches an asymptotic value 
with increasing number of steps. In the asymptotic limit, entanglement value 
decreases with increase in 6 and this dependence can be attributed to spread of the 
amplitude distribution in position space. That is, with increase in 6, constructive 
interference of quantum amplitudes towards the origin gets prominent narrowing 
the distribution in the position space. In Figure 5.2, the process is repeated for 
a particle initially in a non symmetric superposition state = [cos(40°)|0) + 

e^30°sin(40°)|l)] ® l^po). 
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5.2.2 Spatial Entanglement 

Spatial entanglement is the entanglement between the lattice points. This entan- 
glement takes the form of non-local particle number correlations between spatial 
modes. To observe spatial entanglement we need to first associate the lattice with 
state of the particle. 

5.2.2.1 Using single particle quantum walk 

In a single particle quantum walk, each lattice point is associated with a Hilbert 
space spanned by two subspaces. The first is the zero-particle subspace which 
does not involve any coin (particle) states. The other is the one-particle subspace 
spanned by the two possible states of the coin, |0) and To get the spatial 
entanglement we will write the state of the particle in the form of the state of 
the lattice. Following from (5.2) the state of the particles after first two steps of 
quantum walk takes the form : 

1^2) = Wel^i) 

= 7[cos(^)|0)|^,_2) + sin(e)|l)|^,)] 

+5[sin(^)|0)|^,)-cos(^)|l)|V;,+2)]. (5.4) 

In order to get the state of the lattice we can redefine the position state in the 
following way: the occupied position state {tpj) as which means that the j-th 
position is occupied and rest of the lattice is empty. Therefore, we can rewrite 
(5.4) as, 

= 7[cos(^)|0)lV2) + sin(^)|l)|l,)] 

+5[sin(e)|0)|l,)-cos(e)|l)|l,+2)]. (5.5) 

When j— th position is unoccupied (empty) the lattice state is written as \0j). 
Since we are interested in the spatial entanglement we project this state into one 
of the coin state so that we can ignore the entanglement between the coin and the 
position state and consider only the lattice states. Here we will choose the coin 
state to be |0) and take the projection to obtain the state of the lattice in the form 
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l^iat) = |0)(7cos(^)|l,-2) + 5sin(^)|l,)). (5.6) 

Each lattice site j can be considered as a Hilbert space with the basis state 
(occupied state) and |0j) (unoccupied state). Then, the above expression (5.6) in 
the extended Hilbert space of each lattice can be rewritten in terms of occupied 
and unoccupied lattice states as 

|M/;j=7cos(^)lV20,)+5sin(e)|0,_2l,). (5.7) 

We can see that after first two steps of quantum walk the lattice states and 
\lj-2) are entangled. One can check that the lattice states and llj+2) are 
entangled if we choose the coin state to be With increase in number of steps 
the state of the particle spread in position space and the projection over one 
of the coin state reduces the state considered to measure spatial entanglement. 
Therefore, with increase in number of steps the spatial entanglement from single 
particle quantum walk decreases. 

5.2.2.2 Using many particle quantum walk 

We will extend the study of the evolution of spatial entanglement as the quantum 
walk progresses on a many particle system. 

To define a many particle quantum walk we will consider a one dimensional lattice 
with one particle at each position as initial state. Figure 5.3. M independent 
identical particles in M lattice with each particle having its own coin and position 
Hilbert space will have a total Hilbert space n = inc® Hpf^. We consider the 
particles to be distinguishable for the time being. 

The evolution of each particle is independent as they are not interacting, and hence 
the evolution operator is simply Wf^'^. The initial state that we will consider for 
many particle system in one dimension will be 
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Figure 5.3: Many particle state with one non- interacting particle at each 

position space. 



For an M particle system after t steps of quantum walk, the Hilbert space consists 
of the tensor products of single lattice position Hilbert space which are (2t+M + l) 
in number. That is, after t steps of quantum walk, each of M particles spread 
between (j — t) to (j + 1). In principle, each lattice point can be empty without 
any particle or can have a particle in state |0) or That is, at each lattice 
point, each of M distinguishable particles are spanned by the state |0), |1) and 
no particle state. Therefore, after redefining the particle-position Hilbert space in 
the form of state of the lattice, the dimension of each lattice point will be 3*^ and 
the dimension of total Hilbert space is (3^)®*'^. 

Let us first consider the analysis of the first two steps of Hadamard walk {6 = 7r/4 
in (5.1)) on a three particle system with initial state 




We will label the three particles at positions —1, and 1 as A, B and C. Since 
the evolution of these particles are independent, we write down the state after the 
first step as a tensor product of each of the three particle, 

\^?) = wn^l') = m\-2)+5\m]^ 

®[7|0)|-1) + 5|1)| + 1)]3 

®[7|0)|0) + (5|l)| + 2)]c (5.10) 
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where 7 = {1 + i)/2 and S = {1 — i)/2. After two step the tensor product of each 
of the three particle is given by 



7 



7 



|0)|-3) + |l)|-l) 
V2 



d 



|o)|-i)-|i)l + i) 

V2 



J A 



|0)|-2) + |l)|0) \ ^ ^/ |0)|0)-|l)| + 2) 



V2 



V2 



|0)|-1)+J1)I + 1) ^ , ^^|0)l + l)-|l)l + 3) 



V2 



V2 



(5.11) 



By projecting this state onto a 1-D coin state we can get the state of the lattice 
to calculate spatial entanglement. We choose this state to be |0) ® |0) ® . . . |0). 
Then the state of the lattice after projection and normalization is 



l^lat) =7'|A)-3|B)-2|C)_i 

+ 7^5 ( |A)_3|B)-2|C)i + |A)_3|B)o|C)_i + |AC)_i|B)_2 
+ 752 ( |A)_3|B)o|C)i + |AC)_i|B)o + |A)_i|B)_2|C)i) 
+ 53|A)_i|B)o|C)i. 



(5.12) 



From the above expression the spatial entanglement due to quantum walk on 
a three particle system can be calculated. But with increase in the number of 
particles and number of steps of quantum walk, the measure of entanglement gets 
complicated. For which a scalable entanglement measure scheme we discuss in 
next section would be useful. 



5.3 Calculating spatial entanglement in a multi- 
partite system 

In a system with two particles, the state is separable if we can write it as a tensor 
product of the individual particle states, and entangled if not. However, it is 
somewhat more involved to define entanglement for a system with more than two 
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Number of partices 

Figure 5.4: Evolution of spatial entanglement with increase in the number 
of steps of the quantum walk for different number of particles in an open one 
dimensional lattice chain. The entanglement first increases and with the further 
increase in the the number of steps, the number of lattice positions exceeds 
the number of particles in the system resulting in the decrease of the spatial 
entanglement. The distribution is obtained by implementing quantum walk on 
particles in the initial state "^(|0) + i\l)) and Hadamard operation i?o,7r/4,o ^ 

quantum coin operation. 

particles. In general a state is said to be partially entangled if it can be written as 

\ij) = \(f)i)0\(j)2)CB)---\(f)k) (5.13) 

where k < M. In case k = M the state is said to be fully separable and will 
then denote the state of the z— th particle. If on the other hand k = 1 then the 
state will be fully entangled. 

There are quite a few good entanglement measures for the multipartite state [86, 
157-162]. We will be using the Meyer- Wallach (M-W) entanglement measure in 
the lattice as it does not diverge with increasing system size and it is relatively 
easy to calculate [86]. The M-W measure is the entanglement of a single particle 
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Figure 5.5: Value of spatial entanglement for 20 particles on a one dimensional 
lattice after 20 steps of quantum walk using different values of 9 in the quantum 
coin operation i3o,6»,o- (a) and (b) are the distributions for particles initially in 
state "^(|0) + and state |0) (|1)) respectively. In (a), with increase in 
from to 7r/2 the spread of the distribution in lattice position decreases first up 
to ^ = 7r/4 and again increase. Due to asymmetric distributions for particles in 
case of (b) the above effect is not very prominent. For both (a) and (b) when 
9 = vr/2, for every even number of steps of quantum walk, the system returns 
to the initial state where entanglement is 0. Entanglement is for 9 = 0. 



to the rest of the system, averaged over the whole of the system. For pure states, 
linear entropy can serve as a good measure for the entanglement of a single particle 
with the rest of the system. It is given by 

E=j^[l- Trp2] (5.14) 
for a (i-dimensional particle Hilbert space. The M-W measure follows : 



E 



d 



d 



(5.15) 



1=1 



where L is the system size. In a multipartite quantum walk the dimension of each 
lattice point, after projection over and one of the coin (particle) state, is 2^ where 
M is the number of particles. So the expression for entanglement will be 
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2M / 2 ^"^^ \ 

V J = -('H-f) / 

where t is the number of steps and pj is the reduced density matrix of j*'* lattice 
point. Pj can be written as 

p. = Epil^)(^l (5.17) 

where \k) is one of the 2*'^ possible states available for a lattice point and pi. can 
be calculated once we have the probability distribution of individual particle on 
that lattice. 

Since we have M distinguishable particles, we have 2*^ configuration depending 
upon whether a given particle is present in the lattice point or not after freezing 
the state of the particle. That forms the basis for a single lattice point Hilbert 
space. Now we can calculate p^, the probability of one of the possible \k) state 
of particle in the j— th lattice point as: let us say aj'*'' is the probability of i— th 
particle to be or not to be in the in the j— th lattice point depending on If /j 
is 1 then it gives us the probability of the particle to be in the lattice point. If li 
is —1 then a^j'^ is the probability of the particle not to be in the lattice point so 
a^~^^ = 1 — aj^^. And hence we can write 



Thus the spatial entanglement can be conveniently calculated once we have the 
probability distribution of individual particle on that lattice. Since quantum walk 
is a controlled evolution, one can obtain a probability distribution of each particle 
over the entire lattice positions. In fact one can easily control the probability 
distribution by varying the quantum coin parameters during the quantum walk 
process and hence the entanglement. 

Figure 5.4 is the phase diagram of the spatial entanglement using the many particle 
quantum walk. This figure gives us the information about the entanglement in 
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the space for different number of particles with increasing number of steps of the 
quantum walk. 

Here we have chosen Hadamard operation -Bo,7r/4,o and :^(|0) + z|l)) as quantum 
coin operation and initial state of the particles respectively for the walk evolution. 
To see the variation of entanglement for a fixed number of particles with increase 
in steps, we can pick a line parallel to y axis i.e, for fixed number of particles and 
see the entanglement varying with number of steps. 

We see that the entanglement for an instance, after first iteration goes to maxi- 
mum from zero and with the further increase in the number of steps, the number 
of lattice positions exceeds the number of particles in the system resulting in the 
decrease of the spatial entanglement. The decrease in entanglement before the 
number of steps is equal to number of particles should be noted. This is be- 
cause, for Hadamard walk the spread of the probability distribution after t steps 
is between ^ and ^ and not between —t and t [61]. To see the variation of en- 
tanglement with parameter 6 more clearly, in Figure 5.5 the spatial entanglement 
for 20 particles after 20 steps of quantum walk using different values of 6 in the 
quantum coin operation B^ g o is presented. Variation of entanglement with 6 for 
two case, all particles in equal superposition state :^(|0) + (unbiased quan- 
tum walk) and in one of the basis state |0) or |1) (biased quantum walk) as initial 
state is presented. The effect of biasing the quantum walk by different magnitude 
can also be reproduced using the other two parameters ^, ( in the coin operation 

In Figure 2.3, we can see that for a single particle quantum walk, with increase 
in 6 from to tt/A in the coin operation, along with decrease in the spread, the 
maximum point (two peaks) in the probability distribution in the position space 
also decrease. Whereas, with increase in 6 from n/A to -7r/2, the maximum point 
in the probability distribution in position space increases along with decrease in 
the spread. That is, we can see that the minima of the maximum point in the 
probability distribution is for 6 = tt/A. Therefore, this attributes for the minimal 
in the spatial entanglement at ^ = 7r/4 in Figure 5.5. Note that we have ignored the 
extreme values of ^ = and ir /2 which evolves without any quantum interference 
effect. 
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Figure 5.6: Evolution of spatial entanglement for a system with different 
number of particles in a closed chain. With increase in the number of steps, 
the entanglement value remains close to asymptotic value with some peaks in 
between. The peaks can be accounted for the cross over of the left and right 
propagating amplitudes of the internal state of the particle during quantum 
walk. The peaks are more for chain with smaller number of particles. Increase in 
the number of particles in the system results in the decrease in the entanglement 
value. The distribution is obtained by using "^(|0) + as the initial states 
of all the particles and Hadamard operation -Bo,45°,o as quantum coin operation. 




— 10 particles 
— 14 particles 
— 20 particles 



For fixed number of steps of quantum walk, the spatial entanglement first decreases 
and then it starts reviving as we increase the number of particles. 

Closed chain : 

Since most of the physical system that will be considered for implementation will 
be of definite dimension we extend our calculations to one of the simple example 
of closed geometry, n— cycle. 

When we consider a many particle system in a closed chain, with number of 
lattice position equal to number of particles M , the quantum walk process does 
not expand the position Hilbert space like it does on an open chain. Therefore the 
spatial entanglement does not decrease, remains close to the asymptotic value with 
increase in the number of steps of quantum walk. Figure 5.6 shows the evolution 
of entanglement for a system with different number of particles in a closed chain. 
The peaks seen in the plot can be accounted for the cross over of the left and right 
propagating amplitudes of the internal state of the particle during the quantum 
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Figure 5.7: Value of spatial entanglement for 20 particles on a closed chain 
after 20 steps of quantum walk using different values of 9 in the quantum coin 
operation -Bo,6»,o- (a) and (b) are the distributions for particles initially in state 
^(|0) +i|l)) and state |0) (|1)) respectively. Since the system is in the closed 
chain, the quantum walk does not expand the position Hilbert space therefore 
for all values of 6 from 0° to tt/2 the entanglement value remains close to the 
asymptotic value except for a small peak at smaller values of 9. For 9 = 0° 
when number of steps equal to number of particles, the amplitudes goes round 
the chain and returns to its initial state making the entanglement and for 
9 = tt/2, for every even number of steps of quantum walk, the system returns 
to the initial state where entanglement is again 0. 

walk process. Therefore the frequency of the peaks are more for smaller number of 
particles (smaller closed chain). In general, one can also notice that the increase 
in the number of particles and number of lattice point in the closed cycle results 
in the decrease in the spatial entanglement of the system. 

In Figure 5.7, the value of spatial entanglement for 20 particles on a closed chain 
after 20 steps of quantum walk using different values of 6 in the quantum coin 
operation i?o,e,o is presented. For all values of 6 from to 7r/2 the entanglement 
value remains close to the asymptotic value except for the extreme values of 6. 
For ^ = when number of steps equal to number of particles the amplitudes goes 
round the ring and returns to its initial state making the spatial entanglement 
value = 0. For 6 = 7i/2, for every even number of steps of quantum walk, the 
system returns to the initial state where the spatial entanglement is again 0. 
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Therefore, spatial entanglement on a large lattice space can be created, controlled 
and optimized for maximum entanglement value by varying the quantum coin 
parameters and number of particles in the multi particle quantum walk. In [152] a 
scheme to implement a scalable quantum information processing using ultra cold 
molecules is proposed. Two different species of atoms held in independent optical 
lattices are entangled by translating the separate lattices to overlap. Implementing 
quantum walk on one of the species and creating spatial entanglement can possibly 
make the above scheme more robust. 
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5.4 Summary 



• We have presented the evolution of spatial entanglement in many particle 
system subjected to the quantum walk process. By considering many parti- 
cles in the one dimensional open and closed chain we have shown that the 
spatial entanglement of the system can be controlled by controlling the dy- 
namics of the quantum walk, the quantum coin parameter, initial state and 
number of steps. Developing on this approach, many particle entanglement 
can be generated and controlled in a many body system using quantum walk 
as a tool. 
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Conclusion and future perspective 



6.1 Conclusion 

In this thesis we briefly reviewed the continuous- and discrete-time quantum walk 
which are structurally identical to Schrodinger and Dirac equations, respectively. 
By simple decoupling analysis of the evolution, structural similarity of the discrete- 
time quantum walk and the Dirac equation was shown. To optimize the discrete- 
time quantum walk, the use of three parameter quantum coin operation -85,0,1^ 
from SU(2) group. Numerical data with supporting mathematical analysis, it was 
shown that parameter 9 can be used to control the variance and parameters ^ and 
C to bias and control the biasing in the walk. The use of quantum coin -B5,6i,c ^01 
walk on an n— cycle and its effect on optimizing mixing time was shown. 

Recurrence is an important phenomenon in the study of dynamics. Recurrence we 
considered was the return of unit probability amplitude at the origin during the 
dynamics. Based on the numerical data and analysis we have shown that the quan- 
tum walk evolution dominated by the interference of quantum amplitudes fails to 
satisfy complete recurrence theorem. However, fractional recurrence characterized 
by the quantum Polya number can be seen. 

Symmetries and effect of noise on quantum walk on a line and n— cycle was stud- 
ied. We considered the noise on the coin space in our study. Variants of quantum 
walks on a line was considered and showed that they are equivalent in the sense 
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that the final positional probability distribution remains the same in each variant 
(symmetry in the probability distribution). In particular, we considered variants 
obtained by the experimentally relevant operations of Z or X applied at each 
quantum walk step, with the symmetry operations given by Z and PRX. Inter- 
estingly we observed that these symmetries were preserved even in the presence of 
noise, in particular, those characterized by the phase fiip, bit fiip and generalized 
amplitude damping channels. This is important because it means that the equiv- 
alence of these variants is not affected by the presence of noise, which would be 
inevitable in actual experiments. However, the symmetry of the phase operation 
under phase noise was intuitive, considering that this noise has a Kraus represen- 
tation consisting of operations that are symmetries of the noiseless quantum walk. 
But, for the PRX symmetry under phase noise, and for any symmetry under other 
noisy channels (especially in the case of generalized amplitude damping channel), 
the connection was not obvious before the analysis was completed. When the 
studies were extended to the walk on an n— cycle, an interesting fact that we 
observed was the breakdown of the symmetries found in the walk on a line. Fur- 
ther, when the noise was introduced, above a certain noise level these symmetries 
were restored. Results to study noise model were supported by several numerical 
examples obtained by evolving the density operator in the Kraus representation 
and analytical proofs of the effect of noise on symmetries were obtained using the 
quantum trajectories approach, which we found convenient for this situation. 

In this thesis two of the applications of quantum walk was considered. Use of 
quantum walk to redistribute atoms in optical lattice and evolution of spatial en- 
tanglement in many body system. We proposed the use of the quantum walk to 
redistribute atoms, and studied its dynamics in an optical lattice and to expedite 
the process of quantum phase transition. The coherent control over the atoms 
using the coin degree of freedom during the evolution of the walk was demon- 
strated. Experimentally realizable noisy channels studied in the dynamics part of 
this thesis was used to show the additional control over the evolution and atomic 
density redistribution. Making use of the stimulated Raman transition, we have 
proposed stimulated Raman kicks to act as shift operator along with rf-pulse as 
coin operation to implement discrete-time quantum walk on Bose-Einstein conden- 
sate (BEC). Our scheme can retain the macroscopic features of the wave packet 
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or can be accordingly modified to implement the walk at an individual atom level 
which will be effective for transition from superfiuid to the Mott insulator state 
and vice versa. 

Entanglement is many body system is one of the important topic of interest. As a 
step towards using quantum walk to generate and control entanglement in many 
body system, both distinguishable and indistinguishable, we have presented the 
evolution of spatial entanglement in many distinguishable particle system. We 
considered Meyer- Wallach measure for our study. By considering many particles 
in one dimensional open and closed chain we have shown that the spatial entangle- 
ment of the system can be controlled by controlling the dynamics of the quantum 
walk, the quantum coin parameter, initial state and number of steps. 

6.2 Future perspective 

In this thesis we considered quantum walk on a line and on an n— cycle. Ex- 
tension of these studies presented to general graph and higher dimension could 
reveal other interesting features of the quantum walk dynamics. For example, 
structural similarity of discrete-time quantum walk with Dirac equation was pre- 
sented. Discrete-time quantum walk in higher dimension and its analysis following 
the decoupling approach used to walk on a line could reveal interesting structural 
similarity. 

Derivation of dependency of variance on the parameter 6 for a walk on a line 
(1 — sin(^))t^ from first principles (Fourier analysis) which has not be done in this 
thesis would be a step in a direction of deriving the dependency of variance in two- 
and higher dimensions. 

The results from symmetries and effects of noise on quantum walk can be extended 
and generalized to walk on most of the symmetric closed graphs. This symmetry- 
topology-noise interplay presented would be of relevance to quantum information 
processing systems, and have wider implications to the condensed matter systems. 

Theoretically, the evolution of the density profile with a quantum walk can be 
used in place of quantum Monte Carlo simulation [124, 125] or the time evolution 
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density matrix renormalization group (t-DMRG) to study the correlation and re- 
distribution of atoms in optical lattices [123]. We expect the quantum walk to play 
a wider role in simulating and expediting the dynamics in various physical sys- 
tems. Similarly, entanglement in many body system using quantum walk could be 
an interesting topic to explore quantum phase transition and quantum annealing 
problems in various physical systems. Presently, the spatial entanglement stud- 
ied in this thesis seems to have very little practical interest. Developing on this, 
further studies can be considered to extend the evolution of spatial entanglement 
using indistinguishable particle quantum walk. However, many particle-position 
entanglement using quantum walk will be an interesting and bigger problem to 
explore and use it to study phase transitions in quantum systems. 

In general quantum walk seems to serve as a tool to understand dynamics in various 
quantum systems which can further motivate to be used for various applications 
in quantum systems. 
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Quantum walk and Klein-Gordon 
equation 

A.l Decoupling the coupled expression 

Getting Equation (2.23) from Equations (2.22a) and (2.22b) 
From Equation (2.22b), solving for ^1/^ we get 

*l(j + l,t) = [^rU, t + 1)- cos{e)^R{j - 1, t)] . 

Therefore 

^L(i, t + l) = [*^,(j - 1, t + 2) - cos(^)^fi(j - 2, t + 1)] . 

By substituting the above for ^^(j + 1, t) and ^^(j, t + 1) in (2.22a), we get (2.23). 

A. 2 Getting the difference operator that corre- 
sponds to the differential operators 

Getting Equation (2.24) from Equation (2.23) 
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The difference Vj operator that corresponds to the differential operator d/dt is 

_ vl/0-,t + |)-v[/(j,t-|) 
h • 

By setting the small incremental time to 1 [h = 1) difference operator 

Vi = ^(j,t + 0.5)-^(j,t-0.5), 

corresponds to the difference operator d/dt. Therefore, the operator d"^ /dt^ will 
correspond to applying the difference operator in each of the above two terms, 
which yields 

^2 ^ 1 ^ [v[/(j,t + 1) - Vl>(j,t)] - [m^J^t) - Vl>(j,t - 1)] 



' h h 

(vl>(j,t + l)-2^(j,t) + vl>(j,t-l) 

when the small incremental time step /i = 1, it corresponds to d"^ /dt^. The differ- 
ence operators Vj and corresponding to corresponds to d/dj and corresponds 
to d'^ /dj'^ are also defined analogously for j keeping t constant. 
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Variation of the variance as a 



function of 6 



For a quantum walk using -Bo,6i,o as quantum coin, after t steps the probability 
distribution is spread over the interval (— t cos(6'), t cos(^)) and shrink quickly out- 
side this region. The height of the distribution with 6 is proportional to sin (6*). 
Therefore the integral of the probability distribution with in the interval can be 
approximated to, 



By approximating the probability distribution to fit the envelop of the quantum 
walk distribution. 



where K{e) = ^ cos(e)[l + cos2(2e)][l + sin(^)]. 

The position j in the interval (— t cos(6'), t cos(^)) in position space can be repre- 
sented clS db function of 6, 




(B.l) 




(B.2) 



j ^ /(0) = tcos(6') sin(0) 



(B.3) 
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where (j) range from — ^ to ^. For a walk with coin i?o,e,0! the mean of the distribu- 
tion is zero and hence the variance can be analytically obtained by extrapolating, 



tcos(6») 



tcos(6») 



P{j)rdj = / P(/(0))(/(0))Y(0)rf0. (B.4) 



1^ [l + cos^(2g)] ^/rw(' - ,o.^(,) V(tcos(e)sin(0))'x(tcos(^)cos(0))rf</(B.5) 



Vi 



coB^(e) siii^(0) 



(T^ ^t^[l + cos2(2^)]cos3(^) 



,X(e)(sm2(</,)-l) „• 2 



sin^(0)cos(0) (B.6) 



a 



tt [1 + cos2(2^)] cos=^(^) / e^('')('-'-^) dr 



tl sm{9) cos3(^) / r e^^^^^"'"^) r dr 



+ cos^(2^)]cos^(^) — 



+ cos2(2^)]cos2(^) 



2K{9) 



1 gA-(e)r2 

1 2^ 



dr 



(B.7) 



The integration of J^^ e^^^^"^^ dr can be done using Taylor expansion, 



e^^'^'dr 

1 

1 - 

l + K{e)r 

1 L 



2 ir(e)3r6 ir(e)V 



n„2n 



2! 



3! 4! 

3 T,-tQ\i 



^ i^(g) K{ef K{ef im 



5*2! 7*3! 9*4! 



+ 



+ 



K{eY 



i2n + 1) * n! 



dr 
(B.8) 
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E 



{2n + 1) *n\ 



E 



2if (^)" 4 n K{eY 
n\ (2n + l)n! 



For large n 



in 
' 2n+l 



2 in the above expression. When n is small / 



Substituting in (B.7), the variance is simplified to, 

2 tt [1 + cos2(2e)] cos3(0) tf [1 + cos2(2^)] cos^ {6) 



a 



2K{e) 



v^cos(^)[l + cos2(2^)][l + sin(^)] 



t2(l -sin(g))(l + sin(g)) 
1 + sin(^) 



(B.9) 



:i - sin(^))t2. (B.IO) 



From the solution obtained empirically through numerical integration we find that 
the variation of Be with 9 also fits (1 — sin (6')). That is, 



= Cet" ^ (1 - s\n{e))e. 



(B.ll) 
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Dipole trap for ^'^Rb atoms using 
light of different wavelengths 



Effective laser detuning A for alkalis can be calculated using, 




where Aj is detuning from Di line. 
Maximum potential depth is calculated using 



where F is natural linewidth, P is laser power, and /q is the saturation intensity 
given by. 



_ Ti^hcV he 
^° - ^A^ - 3 a57' ^^-^^ 
where F = 1/r. The numerical value of saturation intensity Jq for rubidium atoms 
is calculated to be 1.67 mW/cm^=16.7 W/m^. 
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C.l Frequency-detuning and potential depth for 
different wavelength 

Below is the table with numerical values of detuning and potential depth for dipole 
traps using laser lights of different wavelengths. 



wavelength (nm) 


Detuning A 


Potential Depth Uq W/m^ 


1064 


-25 X lo^r 


1.53 X 10-35^ 
5.5 X 10-35-^ 


850 


-6.93 X lO^r 


820 


-3.46 X lO^r 


11.04 X 10-35-^ 


800 


-0.6 X lo^r 


63.66 X 10-35^ 



The ac stark shift creates a potential, U proportional to the light intensity 



U{r,z) = Uo 



l + (^)2 



(C.4) 



r and z being the radial and axial co-ordinates, zn = ^^^^ is the Rayleigh range^ 
at wavelength A and beam waist Wq. W{z) is the beam radius as a function of 
axial position z and is given by. 



, 2 

Z 



W{z) = W^\\\\\—^ \ . (C.5) 



The overall potential due to effect of gravity is given by. 



\Jg{r, z) = -mgr - Uq 



g lV(z)2 



1 + — 

2H 



(C.6) 



^Distance at which the diameter of the spot size increases by a factor of a/2- 
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To calculate the minimum power required to trap ^ Rb atoms at distance z from 
the focal point of the beam in the axial direction of the beam, the above equation 
is differentiated, 



-mg 



4Uor 

w{zy 



(C.7) 



mg W{zY 
4^0 g^fzf 



z 

zr 



(C. 



By substituting the appropriate values for the above equation one can calculate 
the minimum power required to trap atoms at distance z from the beam focus in 
the axial direction. Below is the table with the calculated power required to trap 
^''Rb atoms using laser light with different wavelengths and beam waist. 



For light with beam waist, Wq = 50/im 



Wo = 50yum 


1064 nm 


850 nm 


820 nm 


z = 5 cm 


3.33 W 


450 mW 


210 mW 


z = 2 cm 


30 mW 


37 mW 


26.5 mW 


z =1 cm 


16 mW 


8.4 mW 


4 mW 


z = 0.5 cm 


11.6 mW 


3.8 mW 


1.5 mW 


z = cm 


9 mW 


2.7 mW 


1.4 mW 



For Wo = 100/im 



Wo 


= 100/im 


1064 nm 


850 nm 


820 nm 


z 


= 5cm 


610 mW 


110 mW 


47 mW 


z 


= 2cm 


140 mW 


31 mW 


15 mW 


z 


= 1cm 


91 mW 


23 mW 


12 mW 


z - 


= 0.5cm 


78 mW 


22 mW 


11 mW 


z 


= 0cm 


77 mW 


22 mW 


11 mW 
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For Wo = 200/im 



Wo 


= 200/im 


1064nm 


850nm 


820nm 


z 


= 5cm 


1.53 W 


200 mW 


1.1 W 


z 


= 2cm 


650 mW 


11 


iO mW 


900 mW 


z 


= 1cm 


610 mW 


11 


BO mW 


900 mW 


z - 


= 0.5cm 


600 mW 


11 


BO mW 


900 mW 


z 


= 0cm 


600 mW 


11 


BO mW 


900 mW 



The maximum photon scattering rate T^c is given by 

r -I^ 
A h- 



(C.9) 



For the laser hght of 850nm the photon scattering rate at beam waist 3.8/im and 
power 0.12 mW will be, 



r 5.5 X 10-35 f P \ fi 

r— — — — 7 = 6.25 X 10 photons/sec. 

6.93 X lO^r 1.054 X lO'^^ \W^ ' ^ ' 



(C.IO) 
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